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EXECUTIVE  SUMMARY 


The  Air  Force  Dispersion  Assessment  Model  ("ADAM")  was  developed  in 
«  1987  to  assess  the  potential  hazard  areas  arising  from  the  accidental 

release  of  several  chemicals  of  interest  to  the  U.S.  Air  Force 
( "USAF" ) .  The  models  developed  primarily  focused  on  simulating  the 
V  dispersion  in  the  atmosphere  of  vapors  and  any  entrained  liquid 

aerosols.  Dispersion  regimes  in  both  the  heavy  gas  region  and  the 
neutral  density  region  were  modeled.  Reactions  of  the  chemical  with 
ambient  moisture  were  also  modeled  and  considered  in  the  ADAM 
programs.  The  list  of  chemicals  whose  dispersion  behavior  was 
modeled  included  i)  nitrogen  tetroxide,  ii)  chlorine,  iii)  anhydrous 
ammonia,  iv)  phosgene,  v)  hydrogen  sulphide,  and  vi)  sulfur  dioxide. 
Because  of  potential  hazards  associated  with  USAF  and  industry 
handling,  storage,  and  transport  of  hydrogen  fluoride  and  liquid 
fluorine,  it  was  decided  to  add  these  two  chemicals  to  the  list  of 
chemicals  in  the  ADAM  coda.  This  report  deals  with  the  details  of 
modeling  hydrogen  fluoride  and  fluorine  dispersion  in  the  atmosphere. 


SCOPE  OF  WORK 

The  scope  of  work  described  in  thin  report  included: 

e  development/modification  of  thermodynamic  models  for  the  mixing 
of  hydrogen  fluoride  (HF) -humid  air  as  wall  as  fluorine  (F2)  - 
humid  air. 

e  integration  of  the  dispersion  models  for  HF  and  Fz  into  ADAM, 
e  evaluation  of  the  model  results  with  available  field  test  data. 


PROJECT  ACCOMPLISHMENTS 

In  this  project,  the  relevant  thermodynamic  properties  of  hydrogen 
fluorine  and  fluorine  were  gathered,  reviewed  carefully,  and  included 
in  the  ADAM  database.  In  addition,  thermodynamic  models  were 
developed  to  determine  the  final  state  of  the  released  chemical 
(either  HF  or  F2)  after  mixing  adiabatically  with  humid  air  of 
specified  mass.  These  models  were  incorporated  into  ADAM.  Source 
models  describing  the  release  and  formation  of  the  vapor  clouds  wore 
developed  and  integrated  into  ADAM.  Finally,  the  dispersion  of  HF 
and  F2  were  modeled  and  included  in  the  ADAM  routines.  Discussed 
below  are  more  details  of  the  special  behavior  properties  of  the  two 
chemicals  considered  and  analyses  performed. 


(Xi) 


HYDROGEN  FLUORIDE  (HF) 


HF  Properties 

Anhydrous  HF  boils  at  20  °C  at  atmospheric  pressure.  Release  of  this 
liquid  from  tanks  at  higher  than  20  °C  results  in  the  formation  of 
vapor  by  flash  vaporization  and  entrainment  of  liquid  droplets  into 
the  vapor.  In  chis  project,  the  detailed  physical  and  chemical 
properties  of  HF  were  gathered,  tested,  and  incorporated  into  the 
ADAM  database.  Also,  the  flashing  of  HF  has  been  modeled. 

HF  Thermodynamic  Modal 

The  vapor  in  equilibrium  with  liquid  at  20  °C  comprises  a  mixture  of 
HF  vapors  of  different  molecular  weights.  This  is  due  to  the 
polymerization  (or  association)  property  of  HF.  A  thermodynamic 
model  was  developed  to  determine  the  vapor  mass  fractions  of  the 
various  oligomers  in  the  saturated  vapor  at  20  °C  and  from  it  the 
density  of  vapor. 

Hydrogen  fluoride  vapor  not  only  polymerizes  and  dissociates  when 
mixed  with  air,  but  also  reacts  with  the  atmospheric  moisture.  The 
reaction  with  moisture  is  exothermic  while  dissociation  reaction  is 
endothermic.  Also,  the  chemical  properties  are  such  that  at  low 
vapor  pressures  and  high  temperatures  the  oligomer-to-monomer 
reaction  is  favored.  Because  of  this  reactivity  and  associative 
properties,  the  final  thermodynamic  state  of  a  mixture  of  HF  and 
humid  air  varies  considerably  depending  on  the  initial  conditions  and 
the  amount  of  air  mixed.  The  final  conditions  of  interest  to 
dispersion  are  the  cloud  density,  cloud  quality,  and  temperature. 

In  this  project,  an  elementary  thermodynamic  model  developed  by 
Schotte  of  DuPont  Inc.  was  modified  to  extend  its  range  of 
applicability  both  in  HF  initial  temperature  and  HF  concentration. 
Some  of  the  HF  property  correlations  were  modified  to  enhance  both 
the  accuracy  of  prediction  and  the  use  in  a  computer  program  for 
quick  convergence.  The  modified  thermodynamic  model  for  the  mixing 
of  anhydrous  HF  (which  may  or  may  not  contain  aerosols)  with  humid 
air  is  discussed  iu  Chapter  2.  The  temperature  range  of 
applicability  of  the  model  has  been  extended  to  an  initial  HF 
temperature  of  1000  K. 

The  thermodynamic  model  calculates  the  final  mixture  conditions  given 
HF  initial  conditions,  temperatures  and  relative  humidity  of 
atmospheric  air  and  the  mass  of  air  mixed  with  a  unit  mass  of  HF. 
These  final  conditions  include  the  mixture  density  and  temperature, 
mass  fractions  in  vapor  and  liquid  phases,  and  mass  and  mole 
fractions  of  individual  species  of  HF  (various  oligomers)  in  vapor 
and  liquid  phases.  The  model  results  have  been  compared  with  results 
from  a  model  developed  by  Clough,  et  al.  The  agreement  is  good.  No 
experimental  data  are  available  to  check  the  predictions  of  the 
thermodynamic  model  over  the  full  range  of  HF  concentrations  and 
temperatures, 
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HF  Dispersion  Model 


The  HF-Air  thermodynamic  model  has  been  integrated  into  the  ADAM 
code.  The  dispersion  results  from  ADAM  have  been  tested  against  the 
data  from  a  single  series  of  field  tests.  These  tests,  called  the 
Goldfish  Series,  were  conducted  by  an  industry  consortium  headed  by 
Amoco,  Inc.  Three  dispersion  tests  with  anhydrous  HF  releases  were 
conducted  in  a  desert  environment. 

The  ADAM  results  for  the  conditions  of  the  test  and  experimental  data 
are  compared  in  Chapter  4.  The  data  and  results  compared  include  the 
i)  source  rate,  ii)  plume  temperature,  and  iii)  variation  of 
centerline  ground  level  concentration  with  downwind  distance. 

It  is  seen  that  the  measured  and  predicted  source  rates  agree  very 
closely.  The  measured  minimum  temperatures  at  various  downwind 
locations  and  those  predicted  at  the  same  locations  agree  reasonably 
close.  Some  of  the  temperature  data  from  the  test  may  not  be 
accurate  because  their  values  are  below  what  are  predicted  by  the 
thermodynamic  model. 

The  overall  trend  in  the  predicted  concentration  variation  with 
distance  agrees  closely  with  the  test  data.  The  numerical  agreement 
in  the  predicted  vs.  measured  concentration  data  are  close  within  a 
factor  of  2  in  Tests  #1  and  #3  and  within  a  maximum  deviation  of  a 
factor  of  5  in  Test  #2 .  A  number  of  reasons  as  to  why  there  are 
discrepancies  in  the  predicted  vs.  measured  data  values  are  discussed 
in  Chapter  4.  Some  of  the  errors  may  be  due  to  significant  scatter 
in  the  concentration  data  reported  from  the  tests.  Overall,  however, 
ADAM  predictions  are  reasonable  representation  of  the  test  results 
within  the  accuracy  needed  for  hazard  area  estimation. 


FLUORINE  (Fz) 
f2  Properties 

Fluorine  is  a  highly  reactive  chemical  which  is  in  the  gaseous  state 
at  ambient  pressure  and  temperature.  It  boils  at  84.5  K  at  ambient 
pressure  and  is  highly  reactive  with  most  substances.  Fluorine  at 
normal  temperatures  reacts  with  water  (vapor  or  liquid)  forming  the 
fluorides  of  hydrogen  and  oxygen.  However,  there  is  no  evidence  of 
occurrence  of  this  reaction  when  the  fluorine  temperature  is  at  or 
near  its  boiling  point  (84.5  K) . 

Fluorine  is  normally  shipped  as  a  liquid  in  double-tanked  road  or 
rail  tankers.  The  outer  tank  contains  liquid  nitrogen  (at  77.4  K) 
which  maintains  the  F2  as  a  liquid  in  the  inner  tank.  Nitrogen 
temperature  is  maintained  at  its  boiling  point  by  allowing  a  small 
amount  of  evaporation  to  be  vented  to  the  atmosphere. 


P2  Release  scenario* 

Two  release  scenarios  have  been  considered  in  this  report,  namely, 
the  tank  burst  scenario  ("Scenario  A")  and  the  tank  breach  scenario 
(  "Scenario  B")  .  The  former  is  assumed  to  occur  when,  due  to  an 
accident,  the  outer  tank  is  breached  and  all  of  the  liquid  nitrogen 
is  lost.  The  inner  tank  heats  up  and  when  the  internal  fluorine 
pressure  reaches  the  tank  failure  pressure  all  of  the  fluorine  is 
released  instantaneously.  The  fluorine  released  flashes  forming  a 
cloud  of  vapor  and  liquid  aerosols  at  a  temperature  of  84.5  K.  Also, 
the  violence  of  release  will  entrain  air  from  the  ambient.  It  is 
assumed  that  the  amount  of  air  entrained  is  about  10  times  the  volume 
of  the  vapor  volume  produced. 

In  Scenario  B,  the  breach  occurs  simultaneously  in  both  the  outer  and 
inner  tanks  leading  to  the  release  of  both  liquid  nitrogen  and  liquid 
fluorine.  These  liquids  form  a  pool  on  the  ground  and  evaporate, 
rather  quickly,  because  of  the  very  cold  temperature  of  each  liquid. 
It  is  estimated  that  the  liquid  pool  formed  by  the  release  of  the 
entire  contents  of  a  road  tanker  of  fluorine  will  evaporate  within  13 
seconds.  The  initial  vapor  cloud  formed  will  consist  of  only  vapors 
of  fluorine  and  nitrogen  (at  77.4  K) . 


7 2  Thermodynamic  Model 

The  mixing  of  humid  air  and  cold  fluorine  vapor  containing  liquid 
aerosols  has  been  modeled.  It  is  assumed  that  there  are  no  chemical 
reactions  between  the  water  and  fluorine  because  of  the  very  cold 
temperature.  The  condensation  and  freezing  of  water  vapor  has  been 
accounted  for.  Also,  the  evaporation  of  the  fluorine  liquid  aerosols 
has  been  modeled.  In  the  case  of  mixing  of  pure  vapor,  initially 
diluted  with  nitrogen  vapor  a  similar  themodynamic  modeling  approach 
as  teha  above  is  used.  The  final  thermodynamic  state  of  the  mixture 
is  determined  and  used  in  the  subsequent  dispersion  calculations. 

It  is  seen  that  the  density  of  fluorine  vapor  at  saturated  condition 
at  ambient  pressure  is  about  5.5  kg/nr  3  which  is  about  5  times  the 
density  of  air.  The  presence  of  liquid  aerosols  and  condensed  or 
frozen  water  makes  the  cloud  even  heavier.  The  mixture  density 
decreases  monotonically  with  increased  air  dilution  in  both  the 
aerosol  case  and  in  the  initial  vapor  case. 


7 j  Dispersion  Modal  Results 

We  have  calculated  the  dispersion  behavior  of  fluorine  released  under 
scenario  A  and  scenario  B.  In  each  case,  the  release  is  assumed  to 
be  from  the  road  tanker.  It  is  further  assumed  that  both  releases 
could  be  classified  as  "instantaneous"  releases  because  in  scenario 
B  the  duration  of  evaporation  of  the  pool  formed  on  the  ground  is 
extremely  short. 
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The  short  term  exposure  (concentration)  limit  for  fluorine  is  very 
low,  and  is  1  ppm.  It  is  seen  that  the  dispersion  distances  to 
achieve  this  level  of  concentration  are  very  large,  in  tens  of 
kilometers!  This  is  because  not  only  is  the  entrainment  reduced  due 
to  the  very  heavy  (and  therefore  stratified)  cloud,  but  also  the 
degree  of  dilution  required  is  of  the  order  of  magnitude  106.  In  the 
case  of  neutral  atmospheric  condition  at  release,  the  hazard  distance 
is  calculated  to  be  greater  than  20  km,  and  in  the  case  of  stable 
condition,  the  same  distance  increases  to  60+  km!  These  results 
imply  that  during  a  significant  portion  of  the  dispersion  the  cloud 
can  be  considered  to  be  neutrally  buoyant.  As  a  matter  of  fact,  the 
downwind  hazard  distance  can  be  calculated  with  the  simple  point 
source  Gaussian  model  without  losing  much  accuracy. 


CONCLUSIONS 

Based  on  the  work  indicated  in  this  report,  it  is  concluded  that: 

1.  The  thermodynamic  model  developed  for  HF  is  robust  and  is 

applicable  over  a  wide  range  of  HF  temperatures  (290  K  to  1000 
K)  and  initial  conditions  (pure  vapor,  vapor  mixed  with  liquid 
aerosols)  and  atmospheric  humidiy  conditions  (0%  to  100% 

relative  humidity) . 

2.  While  the  reactivity  of  Fz  is  very  high  at  ordinary 

temperatures,  it  is  unlikely  to  have  any  reactions  with  water 
vapor  in  the  atmosphere  when  released  from  a  saturated  liquid 
condition. 

3.  The  integration  of  the  source  models  and  atmospheric  dispersion 
models  for  both  HF  and  F2  into  ADAM  has  been  accomplished 
successfully. 

4.  The  ADAM  predicted  concentration  variation  with  downwind 

distance  agrees  reasonably  well  with  data  from  HF  field  tests. 

5.  Disagreements  between  HF  test  data  and  predictions  have  been 
discussed  and  explained  to  the  extent  possible  (with  publicly 
available  information) . 

6.  There  are  no  experimental  dispersion  data  for  F2  with  which  to 
compare  the  ADAM  results.  Hence,  no  comparisons  have  been  made. 
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CHAPTER  1 


INTRODUCTION 


1.1  Background 

The  hazards  of  toxic  vapors  arising  from  the  accidental  releases  of 
chemicals  have  received  considerable  public  attention  because  of 
the  potential  for  affecting  large  areas  and  large  numbers  of 
people.  A  number  of  experimental  studies  and  mathematical  modeling 
analyses  have  been  undertaken  both  in  the  U.S.  and  abroad. 
Research  is  continuing  to  understand  various  aspects  of  the 
dispersion  of  chemicals  including  modeling  different  types  of 
chemical  sources,  effect  of  properties  of  chemicals  on  the 
dispersion  phenomenon,  effects  of  terrain  and  obstructions  in  the 
path  of  dispersing  clouds/plumes,  etc.  The  study  documented  here 
forms  part  of  the  continuing  research  efforts  sponsored  by  the  U.S. 
Air  Force  (USAF)  to  understand  and  model  chemical  vapor  dispersion 
so  that  both  contingency  planning  and  emergency  response  can  be 
effectively  implemented. 

Until  a  decade  ago,  most  of  the  mathematical  models  used  for  gas 
concentration  predictions  were  based  on  Gauss's  theory  of  diffusion 
of  ensemble  of  particles  (the  "Gaussian  Model") .  It  has  become 
obvious,  from  the  several  field  experiments  involving  the  release 
of  heavy  gases,  that  the  conventional  approach  of  using  the 
Gaussian  model  to  describe  the  dispersion  behavior  of  gases  and 
vapors  that  are  heavier  than  air  ("heavy  gases")  leads  to 
significantly  distorted  results  for  the  hazard  area,  downwind 
distance  extent  of  hazard  and  the  duration  of  the  hazard.  The 
vapors  of  many  chemicals  that  are  handled  by  the  USAF  are  heavier 
than  air.  Some  of  the  chemicals,  when  released  from  containment, 
form  liquid  aerosols  which  adds  to  the  negative  buoyancy  of  the 
vapor  clouds.  Some  may  even  react  with  the  ambient  moisture  or 
polymerize  upon  release  into  air  forming  new  species. 

A  heavy  gas  dispersion  model  was  recently  developed  under  the 
sponsorship  of  the  USAF  (Raj  and  Morris,  1987) .  This  model,  called 
the  Air  Force  Dispersion  Assessment  Model  (ADAM) ,  models  the 
release  of  liquid  or  gaseous  chemicals  into  the  environment,  and 
the  dispersion  of  the  vapor  clouds  formed  taking  into  consideration 
the  negative  buoyancy  caused  by  gas  density  and  the  density  of 
liquid  aerosols,  if  any.  The  model  also  takes  into  account  the 
effects  of  any  chemical  reaction  with  the  ambient  moisture.  The 
model  has  been  tested  against  several  field  test  results  and  found 
to  give  reasonable  estimates  of  the  cloud  behavior  and 
concentration  contours.  ADAM  is  expected  to  be  used  in  conjunction 
with  the  Air  Force  Toxic  Chemical  Dispersion  Model  ("AFTOX").  The 
original  version  of  ADAM  contains  the  behavior  models  for  six  (6) 
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chemicals,  nitrogen  tetroxide,  phosgene,  ammonia,  chlorine, 
hydrogen  sulphide  and  sulphur  dioxide. 

The  ADAM  code  is  versatile  in  that  the  dispersion  of  vapors  of 
chemicals  that  are  not  a  part  of  the  original  six  chemicals  can  be 
incorporated  into  the  model  provided  the  physical  and  chemical 
properties  of  the  chemicals  are  provided  in  the  ADAM  database.  In 
addition,  if  the  chemical  is  reactive  or  exhibits  other  types  of  * 

behavior,  then  it  is  necessary  to  develop  thermodynamic/reaction 
models  for  incorporation  into  ADAM.  Two  chemicals  of  interest  to 
the  USAF,  namely,  hydrogen  Fluoride  (HF)  and  fluorine  (F2)  were 
added  to  the  list  of  chemicals  in  ADAM.  HF  is  used  in  lasers  and 
in  alkylation  units  in  the  petroleum  industry.  F.  is  transported 
in  bulk  quantities  as  a  liquid  in  highway  and  rail  tankers. 

This  report  documents  the  mathematical  models  developed  to  simulate 
the  dispersion  of  HF  and  F2  in  the  atmosphere.  The  thermodynamic 
and  reactive  properties  of  the  two  chemicals  are  indicated  in  the 
report  and  these  have  been  considered  in  the  dispersion  models. 


1.2  Hydrogen  Fluoride  and  Fluorine 
1.2.1  Hydrogen  Fluoride 

The  dispersion  behavior  of  hydrogen  fluoride  in  the  atmosphere  is 
very  complex.  HF  has  the  property  of  molecular  association  and 
dissociation.  Depending  on  the  temperature  and  HF  partial 
pressure,  the  HF  forms  a  mixture  of  monomers  and  high  molecular 
weight  polymers.  In  addition,  HF  can  react  with  moisture  in  the 
atmosphere  forming  aqueous  HF.  These  properties  have  profound 
effects  on  the  density  of  the  dispersing  cloud  and  hence  on  the 
process  of  dispersion. 

Large  quantities  of  HF  are  transported  in  bulk  tanks  as  an 
anhydrous  liquid.  The  normal  boiling  point  of  HF  is  20°  c. 
Shipments  of  HF  are  made  both  by  rail  and  highway.  Rail  tank  car 
sizes  are  20  ton,  40  ton  and  80  ton,  whereas  highway  shipments  are 
generally  in  20  ton  quantities.  Occasionally,  HF  is  transported  in 
a  variety  of  cylinder  sizes  (200  lb.,  850  lb.,  1300  lb.,  etc.). 
Tanks  are  of  single  skin  construction  with  no  insulation.  Liquid 
valves  (generally  two  for  each  tank)  dip  into  the  liquid  space  and 
the  vapor  vent  valves  are  at  the  top  of  the  tank.  In  general,  the 
valves  are  all  1-inch  size  (nominal). 

Any  leak  from  a  tank  at  a  temperature  greater  than  20°  c  will  lead 
to  flash  vaporization  of  a  part  of  the  liquid.  Also,  some  liquid 
aerosol  may  be  entrained  into  the  dispersing  vapor  cloud  resulting 
in  the  formation  of  a  stratified,  negatively  buoyant  intrusion 
layer  in  the  atmosphere. 
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1.2,2 


Fluorine 


Fluorine  is  a  very  reactive  element;  it  combines  with  most  organic 
and  inorganic  materials  at  or  below  normal  temperature.  Reactions 
with  organic  materials  is  highly  exothermic;  hydrogen  containing 
compounds  react  with  fluorine  explosively.  While  this  element 
reacts  with  water  or  water  vapor  at  ambient  temperature  to  form 
oxygen  fluoride  and  hydrogen  fluoride,  it  is  uncertain  whether  such 
reactions  are  possible  at  very  low  temperatures. 

Fluorine  is  transported  in  bulk  as  a  liquid  at  very  low 
temperatures  (at  liquid  nitrogen  temperature)  and  ambient  pressure. 
F,  boils  at  about  85  K  at  normal  atmospheric  pressure.  Normally, 
about  2260  kg  (5,000  lbs)  of  fluorine  is  carried  in  each  road 
tanker.  The  chemical  is  maintained  in  a  cryogenic  liquid  state  by 
surrounding  the  fluorine  tank  with  a  liquid  nitrogen  jacket.  About 
1.3  nr  of  liquid  nitrogen  (LN2)  surrounds  the  liquid  fluorine  (LF,) 
and  is  maintained  at  77  K  temperature  by  venting  the  nitrogen  boil- 
off. 


Because  of  the  double  tank  construction  of  the  road  transport, 
accidental  releases  of  F,  can  occur  only  due  to  puncture  of  both  LN2 
and  LF2  tanks.  A  loss  or  LN2  cooling  can  result  in  the  heating  up 
of  LF2  ultimately  resulting  in  a  tank  burst  and  the  release  of  the 
entire  F2  content  explosively.  In  this  latter  release  scenario, 
fluorine  vapor  and  liquid  aerosols  are  formed  and  a  two-phase  vapor 
cloud  will  disperse  in  the  atmosphere.  This  vapor  cloud  is  very 
cold  at  release  (85  K)  and  will  condense  the  water  vapor  from  the 
ambient  into  ice.  It  is  uncertain  whether  there  will  be  any 
reaction  between  the  Fz  and  condensed  water. 


The  original  version  of  ADAM  did  not  treat  the  dispersion  of  HF  or 
F,.  The  development  of  the  appropriate  thermodynamic  and 
dispersion  models  for  both  chemicals  and  their  integration  into 
ADAM  is  described  in  this  report.  In  addition,  several 
improvements  to  the  ADAM  program  and  the  correction  of  certain 
features  was  undertaken  in  this  study  and  these  changes  are  also 
documented. 


1.3  objectives  of  the  work 

The  objective  of  the  work  was  to: 

o  develop  appropriate  thermodynamic  and  dispersion  models  to 
describe  the  dispersion  behavior  of  (each  of)  hydrogen 
fluoride  and  fluorine  chemicals  released  accidentally. 

o  integrate  the  models  into  the  ADAM  system. 


1-3 


1.4  scope  of  Work 


The  following  tasks  were  performed  in  order  to  achieve  these 

objectives: 

1.  HF-humid  air  thermodynamic  models  available  in  the  literature 
were  reviewed.  The  model  selected  was  improved  by  expanding 
its  applicability  to  a  wider  range  of  conditions. 

2.  The  HF-thermodynamic  program  was  integrated  into  ADAM. 

3.  Fluorine  release/source  models  were  developed. 

4.  Thermodynamic  model  describing  the  mixing  of  cold  fluorine 
two-phase  mixture  with  ambient  air  was  formulated  and  solved. 

5.  The  F2  release  model,  thermodynamic  model,  and  dispersion 
model  were  integrated  into  ADAM. 

6.  Routines  in  ADAM  related  to  the  determination  of  atmospheric 
stability  were  improved. 

7.  ADAM  was  modified  to  take  into  consideration  large  aerodynamic 
roughness  effects. 

8.  The  input  and  output  features  of  ADAM  were  substantially 
improved  to  make  them  very  user  friendly. 


1.5  Organisation  of  the  Report 

The  general  physical  and  chemical  properties  of  hydrogen  fluoride 
(HF),  as  they  relate  to  the  dispersion,  are  discussed  in  Chapter  2. 
Included  in  Chapter  2  are  the  details  of  the  thermodynamic  model 
for  the  hydrogen  fluoride  -  moist  air  system  and  results  from  the 
thermodynamic  model.  The  modifications  made  to  the  ADAM  code  are 
discussed  in  Chapter  3  including  the  improvements  to  the  ADAM  input 
panels.  The  results  from  ADAM  using  the  Goldfish  (HF)  field  tests 
release  data  are  compared  with  sparse  data  available  from  the  tests 
in  Chapter  4. 

The  fluorine  (Fj)  models  are  discussed  in  Chapter  5.  In  this 
chapter  are  given  the  details  of  fluorine  transport,  fluorine 
properties,  release  scenarios,  the  thermodynamic  model  for  fluorine 
mixing  with  humid  air,  etc. 

Also  discussed  are  the  results  from  the  dispersion  modal. 

Discussions,  conclusions,  and  recommendations  are  indicated  in 
Chapter  6. 


Several  appendices  are  also  included.  These  provide  additional 
details  of  the  model  derivations,  chemical  property  database 
listings  and  FORTRAN  source  codes  for  the  thermodynamic  model.  In 
Appendix  A  is  illustrated  the  calculations  of  HF  flashing.  The 
distillation  and  evaporation  of  the  two  liquid  (LP2  and  LN2)  pool  is 
described  in  Appendix  B.  The  printout  of  the  thermodynamic 
properties  of  hydrogen  Fluoride  and  fluorine  are  given  in 
Appendix  C. 
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CHAPTER  2 


HYDROGEN  FLUORIDE  -  MOI8T  AIR  MIXING 
THERMODYNAMIC  MODEL 


In  this  chapter,  we  present  a  thermodynamic  analysis  to  predict  the 
physical  state,  temperature  and  density  of  a  mixture  of  anhydrous 
hydrogen  fluoride  (HF)  and  humid  air.  First,  we  discuss  the 
physical  and  chemical  properties  of  HF  followed  by  the  presentation 
of* the  details  of  the  thermodynamic  model. 


2.1  Anhydrous  Hydrogen  Fluoride  Physical  Properties 

Hydrogen  fluoride  (HF)  is  an  extremely  reactive  chemical  which  is 
gaseous  at  ordinary  temperatures  (above  20°  C)  and  atmospheric 
pressure.  HF,  also  known  as  hydrofluoric  acid,  is  a  very  strong 
acid  which  can  severely  burn  skin  and  eyes  and  may  be  fatal  if 
swallowed.  Liquid  anhydrous  HF  and  water  solutions  above  70%  HF 
fume  copiously  and  these  fumes  can  be  dangerous  to  skin,  eyes,  and 
the  respiratory  system  (Allied,  1978) .  Liquid  and  gaseous  forms  of 
HF,  both  anhydrous  and  aqueous,  are  colorless. 

Anhydrous  HF  is  rather  unique  among  hydrogen  halides,  especially 
with  respect  to  one  property:  its  tendency  toward  association  or 
polymerization.  In  both  the  liquid  and  gaseous  states,  anhydrous 
HF  is  believed  to  exist  mostly  as  a  polymer,  though  at  high 
temperatures  and  low  pressures,  the  average  molecular  weight  of 
anhydrous  HF  gas  may  approach  that  of  the  monomer.  At  this  time, 
there  is  conflicting  evidence  as  to  the  nature  of  associated  HF 
molecules.  Some  information  suggests  open  chains  while  other 
information  suggests  cyalic  forms  or  combinations  of  the  two. 
Nevertheless,  it  is  generally  agreed  that  the  degree  of  association 
is  increased  in  the  gaseous  state  as  HF  goes  from  high  temperature 
or  low  pressure  toward  low  temperature  or  high  pressure  or  into  a 
liquid  state.  The  association  property  of  anhydrous  HF  gas  and 
other  factors  cause  it  to  behave  in  a  manner  quite  different  from 
that  of  an  ideal  gas.  Therefore,  even  the  property  measurements 
and  characterizations  are  difficult. 

The  physical  and  chemical  properties  available  in  the  literature 
are  summarized  in  Table  2.1.  Figure  2.1  shows  the  vapor  pressure- 
saturation  temperature  relationship.  The  normal  boiling  point  is 
19.54°  C  which  means  the  chemical  is  in  a  liquid  for  most 
transportation  conditions.  Figure  2.2  shows  the  variation  of 
apparent  molecular  weight  of  HF  vapor  as  a  function  of  temperature 
and  pressure.  From  the  figure,  it  can  be  seen  that  saturated  HF 
vapor  at  the  normal  boiling  point  has  an  apparent  molecular  weight 
of  68  kg/kmol  which  is  about  3.5  times  the  molecular  weight  of  the 
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TABLE  2.1 


SOME  THERMODYNAMIC  PROPERTIES  OP 
ANHYDROUS  HYDROGEN  FLUORIDE  (HF) 


t 

C* 

PROPERTY 

SI  UNITS 

FPS  UNITS 

OTHER  UNITS 

Molecular  Weight  20.00M  kg/kmol  20.0014  lb/lb  mol 

Boiling  Point 'at  292.69  K  62.4  *F  19.54  *  C 

Praoaura 

Critical  Penalty  290  £  SO  kg/m1  16.1  ♦  1.9  lb/ft*  0.29  £  0.03  g/co 

Critical  Prauura  (64.662  ±  3.378)  x  10*  N tut  941  £  1.9  Tb/ft* 

Critical  Temperature  461.15  ♦  3  K  370  +  5  *F  188  +  3  *  C 

Density  of  Liquid  873.5  kg/m1  8.124  lb/gal  0.9735  g/cc 

at  20*C 

Fraccing  Point  189.78  K  -118,07  *F  -83.37  *C 

Heat  of  Formation  of  -13,584  MJ/kg  -5840  Btu/lb  -3.245  kce'l/g 

Ideal  Oaa  at  25*C 

Heat  of  dilution  to  381  kJ/kg  — -  - 

SOX  aquaoua 

Heat  of  dilution  to  731  kJ/kg  ---•  --— 

SOX  aquaoua 


Hast  of  Fusion  at 

Malting  Point 

196.48  kJ/kg 

64.47  Btu/lb 

46.93  oel/g 

Heat  of  Vaporization 
at  NBT 

208.06  kJ/kg 

161.01  Btu/lb 

69.45  cal/g 

Specific  Hast  of  Liquid 
at  NBT 

2524.16  J/kg  K 

0  603  Btu/lb  *F 

0.603  oel/g 

Surface  Tension  at  NBT 

8.6  x  10'1  n/m 

.... 

8.6  dyna/em 

Surface  Tension  at  0*C 

10.2  x  10-'  n/m 

.... 

10,2  dyne/cm 

Refractive  Index  at  25* 
(and  5692.6  A) 

1.1574 

mm  mm 

Toxic  Concentrations  of  Vapor 

Threshold  Limit  Value 

3  ppm  (2  mg/rn’) 

Detectable  by  Small 

2  to  3  ppm 

Throat  Irritation 
Threshold 

5  ppm 

Dangerous  for  Short 

50-250  ppm 

Term  Exposure 


Sources  of  Data: 

1)  Allied  (1978) 

2)  Cham  Ind  Assn  (1976 


Pressure  PSIA 


o  I _ 1 - 1 
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FIGURE  2.1: 
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FIGURE  2.2:  APPARENT  MOLECULAR  WEIGHT  OF  THE  VAPOR 
VERSUS  TEMPERATURE 

Source:  Allied  (1978) 


FIGURE  2.3: 


SATURATED  HF  VAPOR  DENSITY 
VS.  TEMPERATURE 
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Source: 


Allied  (1978) 


monomer!  Figure  2.3  shows  the  variation  of  saturated  vapor  density 
with  temperature.  Saturated  vapor  at  the  normal  boiling  point  has 
a  density  of  3  kg/m3  (about  2.5  times  the  density  of  air  at  20°  c)  . 

2.2  Physical  Processes  in  KF-Humid  Air  Mixing 

In  this  section,  we  discuss  the  physical  processes  occurring  during 
the  mixing  of  anhydrous  HF  at  specified  conditions  with  humid  air 
of  a  given  mass  at  specified  temperature  and  relative  humidity. 
The  mixing  is  assumed  to  occur  at  constant  pressure. 


2.2.1  Description  of  the  Problem 

For  an  accidental  release  of  HF,  the  source  of  the  escaping 
chemical  is  a  storage  vessel  containing  gaseous  and  liquid  HF. 
There  may  be  an  inert  gas  such  as  nitrogen  in  the  HF  source  tank, 
but  this  inert  is  neglected  in  the  calculations  as  it  would  only 
appear  as  a  transient  component  in  any  large  HF  loss.  The  released 
HF  may  be  either  a  vapor  or  a  mixture  of  vapor  and  liquid.  We 
assume  any  liquid  lost  from  the  tank  is  carried  as  aerosol  droplets 
and  is  well  distributed  in  the  HF  vapor  phase.  Since,  at  one  bar, 
HF  boils  at  about  19.5°C  (292.7  K) ,  if  any  liquid  were  present  in 
the  leaking  HF  and  the  ambient  temperature  is  greater  than  19.5°F, 
the  temperature  is  set  at  19.5t’c  due  to  flash  expansion  cooling. 
For  releases  containing  only  HF  vapor,  the  temperature  must  be 
19 . 5°C  or  greater. 

Thus  the  initial  input  to  the  problem  provides  the  specification 
of! 

o  total  mass  flow  of  HF 

o  fraction  of  the  flow  that  is  liquid 

o  temperature  of  the  HF.  If  liquid  is  present  and  the 

ambient  temperature  is  19.5°C  or  greater,  the  temperature 
is  set  to  19 . 5°C 

The  HF  accident  is  assumed  to  occur  into  a  one  bar  environment. 

Upon  release  of  the  HF,  ambient  air  is  entrained.  The  conditions 
of  the  air  that  must  be  specified  are: 

o  mass  ratio  (total  air/total  HF) 

o  air  temperature 

o  relative  humidity  of  the  a.'.r 

The  thermodynamic  analysis  assumes  that  the  only  independent 
variable  is  the  mass  ratio,  J ,e.,  for  all  possible  end  states,  each 
increment  of  the  air  mixed  in  has  the  same  temperature  and  relative 
humidity. 


2.2.2  Hydrogen  Fluoride  Polymerisation 
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One  of  tha  complications  when  dealing  with  hydrogen  fluoride  (with 
or  without  air  dilution)  is  that  the  actual  "chemical"  state  of  the 
HF  is  not  that  of  a  simple  species  but,  due  to  the  interactions 
between  HF  molecules,  the  true  state  is  a  mixture  of  oligomers. 
The  predominant  components  of  this  mixture,  besides  the  monomer, 
are  the  dimer  (HF)2,  the  hexamer  (HF)e,  and  the  octamer  (HF) e,  i.e., 

2HF  ♦♦  (HF)2 
6HF  «  (HF)6 
8HF  ♦*  (HF), 

where  the  double  arrows  indicate  that  a  chemical  equilibrium  is 
attained. 

The  presence  of  such  oligomers  in  the  HF  affects  the  volumetric 
properties  of  the  vapor,  i.e.,  as  more  of  the  higher  oligomers  are 
formed,  the  volume  occupied  by  a  given  mass  of  HF  decreases.  The 
oligomer  concentration  varies  with  pressure  (and  dilution)  as  well 
as  with  temperature.  As  the  pressure  is  lowered  -  or  as  the 
fraction  of  air  increases  so  as  to  decrease  the  partial  pressure  of 
HF  -  the  distribution  of  oligomers  shifts  toward  the  lighter 
entities.  As  a  limit,  at  very  low  HF  pressures,  only  HF  monomer 
remains . 

Similarly,  as  the  temperature  increases,  the  distribution  changes 
to  more  of  the  smaller  oligomers.  Thus  one  can  infer  that  the 
breaking  of  oligomers  is  an  endothermic  process  whereas  the 
formation  is  exothermic. 

Knowing  thi^,  one  can  begin  to  visualize  the  interplay  between  air 
dilution  and  oligomer  distribution.  For  example,  if  dry  air  were 
mixed  with  HF  vapor  (no  liquid)  as  a  simple  case,  the  dilution 
process  would  decrease  the  fraction  of  higher  oligomers  and,  in  so 
doing,  endothermic  reactions  would  cause  a  drop  in  the  mixture 
temperature  until  a  new  equilibrium  state  was  achieved. 

Therefore,  any  thermodynamic  model  must  have  the  ability  to  track 
the  HF  oligomer  distribution  to  maintain  chemical  equilibrium  at 
the  system  temperature  and  HF  partial  pressure. 


2.2.3  Aqueous  HF  Aerosols 

Water  and  HF  are  miscible  in  all  proportions  and  form  aqueous 
non-ideal  solutions  that  show  negative  deviations  from  Raoults's 
law.  What  results  physically  are  solutions  in  which  the  partial 
pressures  of  both  water  and  HF  are  depressed  below  what  might  be 
expected  from  Raoult’s  law.  (Raoult's  law  states  that  the 
component  partial  pressures  are  equal  to  the  liquid  mole  fraction 
multiplied  by  the  pure  component  vapor  pressure.)  The  consequences 
of  such  behavior  are  that  aqueous  HF  aerosols  are  readily  formed 
when  humid  air  encounters  HF  vapor.  This  process  removes  water  and 
HF  from  the  vapor  phase,  liberates  energy  due  to  the  phase  change 
of  condensation,  and  increases  the  density  of  the  resulting 
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vapor-aerosol  mixture.  The  concentration  of  water  and  of  HF  in  the 
aerosol  as  well  as  the  extent  of  condensation  can  vary  over  a  wide 
range.  As  in  any  thermodynamic  analysis,  one  must  allow  for  euch 
events  in  energy,  mass,  and  volume  balances  to  be  certain  that  the 
mixture  is  in  vapor-liquid  equilibrium. 


2.2.4  Adiabatic  vs.  Mon-adiabatic  Mixing 

Most  models  of  the  mixing  of  air  with  a  released  chemical  assume 
there  is  a  negligible  loss  or  gain  of  energy  from  the  environment. 
Assuming  this  is  true  since  the  pressure  is  constant  during  the 
entrainment  of  air,  one  would  base  the  thermodynamic  analysis  on 
the  conservation  of  enthalpy  of  the  species  involved,  i.e.: 

H  (HF)  +  H  (wet  air)  «  H  (mixture)  ■  constant  (2-1) 

where  H  (HF)  and  H  (wet  air)  represent  the  total  enthalpies  of  the 
HF  and  wet  air  streams  before  mixing  whereas  H  (mixture)  would  be 
the  total  enthalpy  of  the  final  mixture.  Equation  (2-1)  will  form 
the  basis  of  the  thermodynamic  analysis  to  follow. 

However,  if  one  desired  to  include  energy  losses  or  gains  from  the 
environment  during  mixing,  Equation  (2-1)  would  bo  modified  to 

H  (HF)  +  a  (wet  air)  +  Q  “  H  (mixture)  -  constant  (2-2) 

where  Q  is  positive  if  energy  is  added  to  the  cloud  end  negative  if 
energy  is  lost.  The  value  of  Q  is  associated  in  some  manner  with 
the  mixing  ratio. 

we  note  that  while  enthalpies  are  state  functions  and  are 
independent  of  the  path  between  the  initial  and  final  states,  the 
parameter  Q  is  clearly  related  to  the  history  of  the  cloud  from  the 
initial  states  of  anhydrous  HF  and  wet  air  to  the  final  mixed 
state.  The  "path"  of  the  cloud  with  mixing  ratios  from  0  to  some 
final  value  affects  the  magnitude  of  Q.  Normally,  however,  this 
path  specific  property  is  ignored  and  the  value  of  Q  employed  is 
based  on  "average"  cloud  properties. 


2.2.5  Summary  of  the  Principal  Assumptions 

Several  assumptions  have  been  noted  or  inferred  in  the  discussion 
to  this  point.  We  summarize  these  below. 

o  the  hydrogen  fluoride  source  is  chemically  pure, 

anhydrous  HF;  vapor  or  vapor  plus  liquid  may  be  present. 

o  the  HF  is  mixed  with  constant  temperature,  constant 

relative  humidity  air  to  form  a  well-mixed  cloud  that  is 
in  chemical  equilibrium  with  respect  to  the  HF  oligomers 
and  in  phase  equilibrium  with  respect  to  the  HF-H?0 
aerosols.  All  kinetic  steps  are  instantaneous  in  nature 
and  no  supersatur ition  states  are  allowed. 
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the  entire  mixing  process  is  isobaric.  If  adiabatic,  the 
total  enthalpy  of  the  system  is  conserved  [Equation  ( 2 - 
1)  ] .  If  energy  is  allowed  to  enter  from  the  environment, 
the  enthalpy  conservation  relation  is  modified  [Equation 
(2-2)]. 


2.3  Review  of  Existing  Thermodynamic  Models 

Two  thermodynamic  models  have  been  developed  to  treat  the  mixing  of 
hydrogen  fluoride  and  wet  air.  These  will  be  compared  and 
contrasted  below.  However,  as  they  differ  primarily  in  the  nature 
of  the  property  correlations  employed  and  in  the  computer  logic,  it 
is  convenient  to  introduce  the  general  program  outline  separately 
so  as  to  be  able  to  emphasize  the  logic  of  the  approaches.  Later, 
when  examining  specific  treatments,  the  differences  can  be  more 
readily  appreciated. 

Recall  that  in  Section  2 . 1  we  indicated  the  general  problem  was  to 
calculate  the  physical  composition  of  a  mixture  of  HF  and  humid  air 
based  on  utilizing  as  input 

e  HF  (flow,  fraction  liquid,  temperature) 

and 

e  Air  (flow,  temperature,  relative  humidity) 

To  achieve  a  solution,  we  employ  Equation  (2-1)  or  Equation  (2-2). 

We  first  must  decide  if  the  final  state  is  all  vapor  or 
whether  there  are  aqueous  HF  aerosols  present.  To  accomplish  this, 
we  assume  that  initially  there  are  na  aerosole  present.  If  this 
assumption  is  proved  false,  then  we  begin  again  but  with  the 
knowledge  that  aerosols  are  present. 

Beginning  with  the  assumption  of  no  aerosols,  a  final  temperature 
is  assumed.  From  the  mixing  ratio  and  input  air  relative  humidity, 
the  final  partial  pressures  of  HF  and  H;0  are  computed.  Note  that 
at  this  stage  all  components  are  assumed  to  be  in  gaseous  phase. 
Using  the  assumed  temperature  and  effective  HF  partial  pressure 
after  dilution,  the  oligomeric  distribution  in  the  HF  is 
calculated.  Comparing  the  component  distribution  with  that  in  the 
initial  HF,  the  degree  of  the  oligomer  reactions  is  found.  An 
enthalpy  balance  1b  carried  out  between  the  initial  and  final 
states  employing  heats  of  reaction  and  heat  capacities  of  the 
various  components.  Closure  of  the  enthalpy  balance  indicates  that 
the  correct  final  temperaturo  was  assumed.  Non-closure  requires  a 
new  choice  of  final  temperature  and  the  computations  are  repeated 
until  there  is  closure. 

Having  a  solution  for  the  single  phase  (vapor)  case,  one  now  tests 
to  ascertain  if  a  liquid  phase  could  exist.  In  all  cases,  vapor- 
liquid  equilibrium  correlations  for  HF  and  water  may  be  expressed, 
in  general,  as 
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(2-3) 

(2-4) 


r* 


Phf  ■ 

Pu  “  f2(T'XHf) 

where  pHf  is  the  effective  partial  pressure  of  hydrogen  fluoride  in 
the  vapor  phase  based  on  the  monomer  molecular  weight,  p  is  the 
partial  pressure  of  water  vapor  in  the  gaseous  phase,  T  is  the 
temperature,  and  xHf  is  the  mole  fraction  hydrogen  fluoride  in  the 
liquid  phase  (again  with  the  monomer  molecular  weight) .  The  form 
of  the  functions  f,  and  f2  in  equations  (2-3)  and  (2-4)  have  been 
developed  empirically  from  experimental  data.  In  the  present 
situation,  values  of  T,  pHF  and  pH  are  available  (with  the 
assumption  of  nfi  liquid  present) . 


Testing  for  a  liquid  phase  can  be  accomplished  in  various  ways.  As 
an  illustration,  suppose  one  employed  Equation  (2-3)  to  compute  xH,. 
knowing  T  and  p^.  Then,  using  this  value  and  the  known  value  of  T, 
with  Equation  (2-4),  one  determines  an  expected  value  of  the  water 
vapor  partial  pressure.  Let  this  calculated  water  vapor  pressure  be 


Pm*  ‘ 


The  liquid  test  is  then 


if  pu*  >  pw,  a  liquid  phase  exists 

if  pw*  <  pH,  no  liquid  phase  exists 


If  no  liquid  phase  exists,  the  problem  is  solved  and  the  computed 
temperature  and  vapor  phase  compositions  are  employed  to  estimate 
the  cloud  density.  However,  if  a  liquid  phase  exists,  the 
calculations  in  Step  l  are  invalid  and  one  proceeds  to  Step  2. 


Step  2 »  we  now  know  that  HF-water  aerosols  are  present  in  the  HF- 
air  cloud,  but  we  do  not  know  their  extent  or  composition.  Nor  do 
we  know  the  mixture  temperature.  The  computational  procedure 
becomes  more  complex,  but,  in  essence,  a  double  trial  and  error 
technique  is  initiated  wherein  both  the  final  temperature  and  final 
HF  composition  in  the  aerosol  are  assumed.  Then,  employing  the 
enthalpy  balance  (Equation  2-1  or  2-2),  the  phase  equilibrium 
relations  (Equations  2-3  and  2-4),  the  chemical  equilibrium 
relationships  for  the  various  HF  oligomers  in  the  vapor,  along  with 
mass  balances,  one  can  arrive  at  a  solution. 


This  computational  procedure  has  been  usad  by  Clough  and  his 
colleagues  at  the  Health  and  Safety  Executive  in  Great  Britain 
(Clough  et  al.,  l9S7a,b)  and  by  Schotte  (1987,  1988)  at  duPont  in 
the  United  States. 


Models  developed  by  both  Clough  et  al.  and  by  Schotte  are 
acceptable.  The  original  Schotte  model  (1987)  was  limited  to  HF 
concentrations  below  20  mole  percent  HF,  but  the  range  was  later 
expanded  to  cover  0  to  100%  HF.  Also,  the  Schotte  model  lacks  a 
satisfactory  single  phase  methodology,  but  the  modal  was  originally 
developed  to  treat  cases  where  HF  aerosols  would  be  present.  As 
noted  later,  the  TMS  version  of  the  Schotte  model  does  have  single 
phase  capabilities.  Both  the  Clough  et  al.  and  Schotte  models  have 
temperature  limitations  in  that  the  vapor-liquid  equilibrium  data 
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upon  which  the  correlations  were  based  did  not  exceed  about  60°C. 

In  treating  the  vapor-equilibrium  of  HF  and  water,  as  noted 
earlier,  the  system  is  highly  non-ideal  and  exhibits  negative 
deviations  from  ideal-solution  behavior.  A  number  of  investigators 
have  studied  this  system  (Munter  et  al.,  1947,  1949;  Vieweg,  1963; 
Tyner,  1949;  Brosheer  et  al.,  1947;  and  Johnson  et  al.,  1973). 
Schotte  (1987)  correlated  the  experimental  data  in  equation  form  to 
relate  the  partial  pressures  of  both  water  and  HF  to  temperature 
and  mole  (or  weight)  fraction  HF  in  the  liquid.  The  partial 
pressures  of  HF  were  based  on  an  effective  mole  fraction  HF  as 
shown  later  in  Eq.  (2-14).  Clough  et  al.  (1987a, b)  have  treated 
aqueous  HF  solutions  by  standard  phase  equilibrium  relationships 
and  employed  a  three  parameter  correlation  for  the  activity 
coefficients  as  a  function  of  temperature  (Wheatly,  1985) .  The 
parameters  were  assumed  composition  dependent.  Values  of  the 
coefficients  are  tabulated  in  their  paper.  Comparing  calculated 
vs.  experimental  partial  pressures  of  HF  and  water  over  a  wide 
temperature  range,  the  relations  given  by  Schotte  were  shown  to  be 
significantly  more  accurate. 

In  the  enthalpy  balances,  Schotte  (1987,  1988)  computed  the 
specific  enthalpy  of  an  aqueous  HF  solution  relative  to  the 
enthalpies  of  pure  vapor  HF  and  HzO  at  a  reference  temperature  of 
25°C.  The  appropriate  quantities  of  HF  and  Hzo  were  first  condensed 
and  then  mixed.  The  enthalpy  of  mixing  was  obtained  from  data  by 
Johnson  et  al.  (1973).  Then  the  solution  temperature  was  raised 
(lowered)  from  25°C  to  the  desired  value  using  solution  heat 
capacities  reported  by  Thorvaldson  and  Bailey  (1946) ,  Kozhevikov  et 
al.  (1982),  and  Franck  and  Spalthoff  (1957). 

In  contrast,  Clough  et  al.  (1987a, b)  employed  their  activity 
coefficient  correlation  and,  by  differentiation  with  respect  to 
reciprocal  temperature,  were  able  to  obtain  partial  molar 
enthalpies  and,  also,  solution  enthalpies  as  a  function  of 
composition  and  temperature. 

When  comparing  calculated  enthalpies  of  mixing  (Schotte,  1988)  vs. 
the  experimental  results  of  Johnson  et  al.  (1973),  Schotte's 
method  was  shown  to  be  more  accurate.  An  illustrated  comparison  is 
shown  in  the  table  below. 
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Enthalpy  of  Mixing,  J/mole  HF 

XHF 

Experimental 

Calculated 

Johnson  et  al.(l973) 

Clough  et  al. 

Schotte 

(1987a, b) 

(1988) 

0.0909 

-18,500 

-24,700 

-18,690 

0.2000 

-17,900 

-24,300 

-18,100 

0.3333 

-16,670 

-22,900 

-16,750 

0.4000 

-15,840 

-21,760 

-15,800 

0.5000 

-14,210 

-19,600 

-14,100 

As  noted  earlier,  HF  in  the  vapor  phase  is  comprised  of  oligomers. 
Clough  at  al.  ( 1937a, b)  only  considered  the  monomer,  dimer,  and 
hexamer  forma  whereas  Sohotte  (1937,  1988)  expanded  the  list  to 
include  octamars  and  an  (HF*H2o)  complex  suggested  by  Thomas 
(1975).  Chemical  equilibrium  constants  for  these  reactions  are 
given  by  Sohotte  (1987)  and  by  Clough  et  al.  (1987).  In  essence, 
these  values  were  derived  from  vapor  density  data  for  HF  assuming 
it  to  bs  composed  of  an  ideal-gas  mixture  of  HF,  (HF)2,  (HF)4,  and 
,  possibly,  (HF) B  and  (HF«H20) .  Other  studies  include  those  of 
Rushmere  (1954),  Maclean  et  al.  (1962),  vanderzee  and  Rodenburg 
(1970),  Jarry  and  Davis  (1953),  and  Armitaga  et  al.  (1963).  There 
is  good  agrssment  between  the  two  models  for  the  chemical 
equilibrium  constants  and  for  the  enthalpies  associated  with  the 
reactions. 

in  ths  two  models,  other  physical  properties  are  necessary.  Vapor 
heat  capacities  (or  enthalpies)  are  required  for  all  components  as 
a  function  of  temperature.  The  vapor  phase  is  normally  considered 
an  ideal-gas  mixture  so  vapor  enthalpies  of  mixing  or 
compressibility  factor  deviations  are  not  required.  The  HF  is 
normally  treated  as  the  monomer  and  polymerization  effects  are 
considered  separately.  For  both  HF  and  water,  vapor  pressures, 
liquid  densities,  and  phase  change  enthalpies  are  also  needed. 
Thera  are  many  sources  of  reliable  data  for  the  required  properties 
in  this  group,  e.g.,  Reid  et  al.,  (1987),  Shaft  et  al.,  (1973), 
Allied  Chem.  Corp.  (1978),  Jarry  and  Davis  (1953),  Yabroff  et  al., 
(1964),  and  Vanderzee  and  Rodenburg,  (1970,  1971).  In  comparing 
the  two  models,  there  is  reasonable  agreement  between  the  physical 
properties  employed  although  the  heat  capacity  of  water  vapor  (at 
constant  volume)  as  used  by  Clough  appears  to  bs  far  too  small 
(3.04  J/mole  K) .  In  running  any  comparisons,  this  heat  capacity 
was  increased  to  25.9  J/mole  K. 

Sohotte  (1988)  has  made  a  few  comparisons  between  the  results  of 
his  model  and  that  of  Clough  et  al.  ( 1987a, b).  In  Figure  2.4, 
calculated  cloud  temperatures  are  shown  for  both  models  for  a  case 
wherein  HF  vapor  (no  liquid)  at  19.5°c  was  mixed  with  air  at  20“C 
and  50%  relative  humidity.  Ths  agreement  between  the  two  models  is 
excellent.  In  Figures  2.5  to  2.7,  for  the  same  case,  comparisons 
are  made  for  the  fraction  of  HF  which  condensed  into  the  aerosols, 
the  computed  aerosol  compositions  and  for  the  cloud  densities.  In 
all  cases,  the  results  from  both  models  are  in  good  agreement. 
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Figures  2.8  through  2.11  show  similar  results  for  a  case  in  which 
HF  containing  50%  liquid  is  mixed  with  air  at  20°c  and  95%  relative 
humidity.  As  before,  there  is  surprisingly  good  agreement  between 
the  two  models. 

In  conclusion,  it  would  appear  that  either  model  could  have  been 
selected  to  serve  as  the  basis  to  develop  the  TMS  version  of  the 
HF-wet  air  mixing  model.  However,  because  the  Schotte  case  seems 
to  provide  more  accurate  physical  property  correlations,  we  chose 
this  model  for  further  development. 


2. 3.1  TMS  Modification  to  the  Schotte  Model 

The  Schotte  model  is  discussed  below.  The  TMS  version  of  this 
model  differs  in  four  areas. 

o  addition  of  a  separate  subroutine  to  test  whether  only 
the  vapor  phase  exists  and  to  calculate  the  final 
temperature,  vapor  composition,  and  cloud  density  of  the 
vapor  only  mixture. 

o  addition  of  a  subroutine  to  model  the  mixing  of  high 
temperature  HF  gas  (1000  K) . 

o  providing  an  option  to  have  an  energy  gain  or  loss  during 
the  mixing  of  the  HF  with  humid  air. 

o  changes  were  made  to  incorporate  the  model  into  ADAM. 


2.3.2  schotte  (1988)  Correlations 

The  modifications  to  the  Schotte  model  are  presented  after  the 
basic  model  structure  and  correlations  are  described.  It  should  be 
noted  that  the  Schotte  model  as  presented  is  significantly 
different  from  tho  one  originally  published  (Sohotte,  1987)  and  is 
contained  in  an  internal  duPont  document  that  was  made  available  to 
TMS  (Schotte,  1988)  . 


Partial  Pressures  of  H7  and  Water  over  Aqueous  Solutions 

The  following  equations  are  correlations  used  to  calculate  partial 
pressures: 

in  pHr  -  AH,  +  aHf/T  (2-5) 

£n  (pu/(l-xHf)  ]  -  Ay  +  By/T  (2-6) 
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where  pKF  and  pH  are  partial  pressures  in  mm  Hg,  T  is  in  Kelvin,  and 
xH,  is  the  mole  fraction  of  HF  in  the  liquid.  xHF  is  based  on  the 
monomer  molecular  weight  (20.01).  The  valid  temperature  range  of 
the  correlations  is  0  to  60°c. 


The  coefficients  in  Equations  2-5  and  2-6  are  determined  using  the 
following  correlations: 


Aue  -  16.9181  4-  21.7958X,,  -  52.3860X  *  +  82.4252XHf5 
HF  -  106.184xHf*  +54.H^291XH,5 

Bu.  -  -5902.78  -  586 . 903XM.  +  1340.82X.,2  +  6822. 09xM.3 
Mf  4  2113. 93XH,4  -  68f8.84xH,s  W 


(*“7) 

(2-8) 


For  XH,  <  0.4738 

AM  -  21.1017  -  3.22961xHr  4  7.90730xHF*  (2-9) 

Bw  -  -5387.02  4  1483. 60xh,  -  4818. 83xH,2  (2-10) 


For  XHF  >  0.4738 

Ay  -  -0.0628905  +  43.2439xHp  +  4.26882XHF2  (2-11) 

By  -  2639.16  -  17849. 4XHf  4  192.177xHF2  (2-12) 


Vapor  Phase  Hydrogen  Fluoride  Composition 

The  hydrogen  fluoride  in  the  vapor  phase  is  assumed  to  be  composed 
of  the  HF  monomer  and  several  oligomers  as  well  as  in  a  complex 
with  water.  We  summarize  the  species  below  with  the  designation 
used  in  the  equations  to  follow. 


Species 

Pressure 

Mole  Fraction 

Molecular  Weight 

Partial 

HF 

yi 

m,  »  20.01 
m2  -  2m, 

P,  -  Yif 

(HF)2 

y2 

p2  -  y2P 

(HF)* 

y4 

m6  ■  6m, 

P4  -  y6J 

(HF)e 

yB 

m8  -  Bm, 

Pb  “  yBpn 

hf*h2o 

y* 

mc  “  38.03 

Pe  “  VcP 

The  effective  mole  fraction  hydrogen  fluoride  in  the  vapor,  yHF,  is 


y„ r  -  Yi  +  y2  +  y6  +  ya  +  ye  (2“13> 

and  the  effective  partial  pressure  of  the  hydrogen  fluoride  (as  in 
Equation  2-5)  is 

P„f  "  yMfP  -  Pi  +  P2  +  Pb  +  Pa  +  Pc  p-m> 
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When  chemical  equilibrium  effects  are  considered,  an  an  example, 


2HF  *♦  (HF)2 


then 


K2  -  f2/f,*  (2-15) 

where  f2  and  f,  represent  the  fugacities  of  the  dimer  and  monomer 
HF.  While  Schotte  employed  fugacitiea  in  hie  treatment,  he 
assumed,  for  all  HF  epecies,  > 

f  j  -  yjd,P 

d1  is  the  fugacity  coefficient  of  the  monomer.  The  numerical 
values  of  0,  ranged  only  from  about  0.98  to  1.00  at  the  low 
pressures  found  in  atmospheric  mixing.  Thus,  while  0,  is  retained 
in  the  computer  program  for  purposes  of  clarity,  we  have  set  0,  - 
1. 


Then,  Equation  2-15  becomes 


k2  -  yj/y^P 

(2-17) 

or 

y,  - 

(2-13) 

Similarly, 

y4  -  yitpIK« 

(2-19) 

y#  -  yi“p7K8 

(2-20) 

For  the  complex, 

HF  +  H20  -  HF 

•  h2o 

Ke  -  fc/(f,fM) 

(yep)/((y1p)  (pw-ycp)  ]  yc/(yiPH) 

(2-21) 

ye  -  yiPA 

(2-22) 

Substituting  Equations  2-18  through  2-22  into  Equations  2-13  and 
2-14, 


Pur  “  Vi*  +  yi*P2K2  +  Yi6P6K6  +  y,"P8KB  +  y^PA  (2'  23) 

Thus,  knowing  pHf  and  pM,  one  can  obtain  a  value  of  y,  and, 
therefore,  y2,  y6,  y8,  and  ye. 
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The  equations  for  the  chemical  equilibrium  constants  are: 


In  K2  -  6,429. 729/T  -  24.14682  (2-24a) 
£n  K4  -  21,101 . 965/T  -  69.73267  (2-24b) 
£n  Ka  »  25,225. 720/T  -  83.47307  (2-25a) 
£n  Ke  -  3,154/T  -  11.425  (2-25b) 


T  is  in  Kelvin  and  K,  has  the  unit  atm'1;  K4,  K8,  and  Ke  have  the 
units  atm’5,  atm"7,  ana  atm-1. 

Material  Balanoes 

Let  M.  be  the  moles  of  air  added.  With  this  air  there  are  Mu  moles 
of  water  where 

Mw  -  pw°MA/(P-pM0)  (2-26) 

and  p°  is  the  partial  pressure  of  the  water  in  the  air  that  is 
entrained. 

With 


mh.  -  moles  of  HF  added  based  on  the  monomer, 

mhp  “  mhpl  +  mhpv  (2-27) 

with  Mhh  and  MHFV  as  the  moles  of  HF  (as  the  monomer)  in  the 
escaping  hydrogen  fluoride  vapor  and  liquid  fractions. 

After  mixing  and  attaining  an  equilibrium  state,  the  hydrogen 
fluoride  moss  balance  is  (xNF  the  mole  fraction  HF  in  the  aerosol 
(as  the  monomer)  and  L  as  the  moles  of  liquid  aerosol) : 

Mh,l  +  MHFV  -  xHFL  +  (yt  +  2y2  +  ey6  +  sy8  +  ye)MA/yA  (2-28) 

where  MA/yA  are  the  moles  of  vapor  since  yA  is  the  mole  fraction 
air.  For  water, 

Mw  -  <l-xHF)L  +  (yw+y0)MA/yA  (2-29) 

If  no  aerosol  is  formed,  L  -  0. 


2.3.3  solution  to  the  Equations  of  the  Model 

It  was  noted  earlier  that  one  must  first  determine  if  an  aerosol  is 
present.  Depending  upon  the  outcome  of  that  test,  different 
enthalpy  balances  are  employed  to  obtain  closure  of  the  enthalpy 
conservation  equation  (Equation  2-1  or  2-2).  If  no  aerosol  is 
formed,  a  single  trial  and  error  procedure  is  used  with  the  final 
temperature  as  the  search  variable.  If  aerosol  is  formed,  the 
procedure  is  more  complex  as  a  double  iteration  is  required  with 
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temperature  and  aerosol  composition  as  the  search  variables.  In 
both  cases,  enthalpies  of  reactions  related  to  variations  in  the 
oligomer  population  must  be  employed.  These  enthalpies  may  be 
obtained  from  the  leading  term  in  Eqs.  (2-22)  through  (2-25)  by 
multiplying  by  (-RT) .  with  R  -  8.314  x  10'3  kJ/mole  K, 

AHj  -  -53.46  kJ/mole  HF 

AH6  -  -175.44  kJ/mole  HF 

AH#  -  -209.73  kJ/mole  HF 

AHe  -  -26.22  kJ/mole  HF 

In  tracking  the  enthalpies  of  the  various  HF  species  for  use  in 
Equation  2-1  or  2-2,  one  employs  the  monomer  as  the  basis  since  all 
polymerization  equilibrium  constants  and  reaction  enthalpies  are 
relative  to  the  monomer.  For  example,  suppose  the  initial  HF 
stream  were  all  vapor  at  soma  Tt  and  P(.  One  first  computes  the 
moles  of  monomer  (N,)  ,  dimer  (N2) ,  ...  in  the  initial  stream.  Then, 
this  stream  is  assigned  a  "chemical"  enthalpy  as 

NjAHj  +  N6AH6  +  NgAHj  +  NeAHc 

which  would  be  the  true  enthalpy  of  the  monomer  if  all  the 
oligomers  were  broken  down  to  the  monomer  to  form  N1  +  2N,  +  6N&  + 
8NB  +  Ne  moles  of  monomer.  Note  that  this  stream  would  then  have 
a  negative  initial  enthalpy  if  the  heat  of  formation  of  HF  monomer 
is  given  a  value  of  zero  at  T(  since  the  breaking  of  oligomers  to 
form  monomer  is  an  endothermic  process. 

For  the  hydrogen  fluoride  enthalpy  in  the  final  state  at  T(,  assume 
that  no  aerosol  has  formed.  There  is  then  a  different  distribution 
of  HF  oligomers  (N,  • ,  N, ' ,  N4',  N#' ,  Nc').  To  determine  the  enthalpy 
of  this  stream,  the  oligomers  are  decomposed  to  the  monomer  as 
before  and  the  total  moles  of  monomer  then  cooled  to  T,.  (Since  no 
aerosol  was  assumed,  the  total  moles  of  HF  monomer  would  equal 
those  initially) .  The  enthalpy  of  this  final  HF  is: 


Nj'AH2  +  N6'AH6  +  N8'AHe  +  Ne'AHe  +  (N, '  +2N?'  +6N*'  +8N# 


'+V)  k 


f 

(HF) dT 


where  Cp(HF)  is  the  isobaric  heat  capacity  of  the  monomer  HF. 

If  aerosols  are  formed,  the  final  vapor  KF  enthalpy  is  computed  as 
above,  but  some  HF  then  exists  in  the  aerosol  and  this  HF  would 
have  an  enthalpy  equal  to  vapor  HF  (monomer)  at  Tf  less  the  molar 
enthalpy  of  condensation  and  the  enthalpy  of  mixing  the  HF  with 
water. 
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a. 3. 4 


Cloud  Density 


Computations  in  tha  computer  program  are  made  using  moles  sinoa 
this  is  more  convenient  to  employ  in  tha  chemical  and  phase 
equilibrium  balances.  For  tha  final  state,  one  has  a  vector  of 
mole  fractions  as  well  as  the  moles  of  aerosol  (L)  and  the  mole 
fraction  HF  in  the  aerosol  (xH,,)  .  Thus,  the  mass  of  the  cloud  per 
mole  of  vapor  phase  is 

W  -  20.01(y1+2yJ+6y4+8y,)+18.02yw+38.03yB+28.97yA  +  (LM'yA/MA)  (2-30) 

where  the  mean  molecular  weight  of  the  aerosol  is 

M'  -  20 .  OlXgp  +  18.02  (1-Xh,,)  (2-31) 

The  cloud  volume  per  mole  of  vapor  is 

V  -  (ZRT/P)  +  LM‘yA/MApL  (2-32) 

Normally,  the  compressibility  factor  Z  is  set  equal  to  1.0.  The 
second  term  in  Equation  2-32  is  small  compared  to  the  first  so  only 
an  approximation  to  pL  is  necessary.  Schotts  set  it  equal  to  1200 
kg/m1,  independent  of  composition  and  temperature.  Finally,  the 
cloud  density  is 

pc  -  W/V  (2-33) 


2 . 3 . 3  Computer  Program 

A  subroutine  called  HFTHRM  was  written,  incorporating  the  procedure 
outlined  in  the  above  seotions,  to  calculate  tha  final 
thermodynamic  conditions  of  a  mixture  of  anhydrous  HF  and  humid 
air.  The  FORTRAN  subroutine  call  statement  and  tha  definition  of 
tha  parameters  are  shown  in  Table  2.2.  The  routine  was  written  to 
be  compatible  with  the  overall  structure  and  calling  conventions 
used  in  ADAM. 


2.4  Results  and  Discussions 

The  HF-humid  air  thermodynamic  model,  HFTHRM,  was  exercised  for 
several  test  cases.  The  results  are  presented  in  a  series  of 
figures,  with  variation  in  the  values  of  selected  parameters.  The 
principal  parameter  varied  is  the  dilution  ratio  (i.e.,  the  ratio 
of  mass  of  air  mixed  with  a  unit  mass  of  hydrogen  fluoride) .  The 
initial  conditions  of  HF  used  include  the  following: 

i)  saturated  HF  vapor  at  ambient  pressure; 

ii)  saturated  HF  vapor  +  liquid  aerosols  at  ambient  prsssure; 

iii)  hiqh  temperature  HF  vapor  at  ambient  pressure. 

Table  2.3  shows  the  parameters  varied  and  the  abscissa  and 
ordinateB  of  the  various  figures. 
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SABLE  2.2 


FORTRAN  SUBROUTINE  FOR  HF-BUMID  AIK  MIZINO 
THERMODYNAMIC  MODEL 

C _ _ 

SUBROUTINE  HFTHRM(MHF,  FL,  THF,  MAIR,  TAIR,  RH,  Q, 

&  TMIX,  RHOMIX,  VOLMIX,  NSPECS,  SPLIST,  CSOL,  CL1Q,  CVAP) 

C _ 

C  This  is  the  TMS  subroutine  for  calculating  thermodynamic 
C  properties  of  the  mixture  resulting  from  the  isobaric  mixing  of 

C  hydrogen  fluoride  and  humid  air.  The  initial  condition  of  HF  can 

C  be  a  saturated  vapor  or  a  saturated  vapor  containing  liquid  aerosols. 
C 

C  The  program  is  modeled  after  the  work  of  william  Schotte, 

C  Ind.  A  Eng.  Cham.  Res.,  76,  p.  300  (1987). 

c 


c 

»** INPUT 

PARAMETERS*** 

**UNITS»* 

C 

c 

MHF 

■ 

Total  masa  of  HF  stream,  vapor  t  liquid 

kg 

(or 

kg/s) 

c 

FL 

m 

Mass  fraction  liquid 

c 

THF 

m 

Temperature  of  HF  etream 

K 

c 

MAIR 

m 

Mass  of  moist  air  stream 

ky 

(or 

kg/s) 

c 

TAIR 

■ 

Dry  bulb  temperature  of  the  air  etream 

K 

c 

RH 

m 

Relative  Humidity  of  the  air  stream 

% 

c 

Q 

m 

Excess  heat  (+  added,  -  extracted)  from 

ay stem  J 

(or 

J/B) 

c 

c  ***OUTPUT  PARAMETERS *  *  * 

C 

n  TMIX  ■  Final  temperature  of  the  mixture  K 

C  RHOMIX  *>  Density  of  the  mixture  kg/m*3  (or  /a) 

C  VOLMIX  ■  Final  volume  of  the  mixtura  ni~3  (or  /s) 

C  NSPECS  -  Number  of  chemical  species  in  the  mixture  *  7 

C  SPLIST  »  List  of  specie  names  in  3  letter  codes 

C  CLIQ(n)  ■  A  vector  of  liquid  mass  fractions  for  n  s  nspecs 

C  CSOL ( n)  »  A  vector  of  solid  mass  fractions  for  n  s  nspecs 

C  CVAP(n)  ■  A  vector  of  vapor  mass  fractions  for  n  s  nspecs 

C  CVAP (nspec+1 )  m  Total  mass  of  all  epeciss  in  vapor  phase  (kg) 

C  CVAP (nspec+2 )  -  Total  moles  of  all  species  in  vapor  phase  (Kmol) 

C  CSOL(nspec+l)  ■  Total  mass  of  all  species  in  Bolid  phase  (kg) 

C  CSOL (nopec+2 )  *  Total  moles  of  all  species  in  solid  phase  (Kmol) 

C  CLIQ(nspec+l)  x  Total  mass  of  all  species  in  liquid  phase  (kg) 

C  CLIQ ( nspec+2 )  *  Total  moles  of  all  species  in  liquid  phase  (Kmol) 

C _ 

C  ***  DEFINITIONS  OF  OTHER  PARAMETERS  *** 

C 

C  SPLIST ( 1 )  -  'H20'J  SPLIST(2 )  -  'AIR'j  SPLIST(3)  -  'HF1' 

C  SPLIST(4)  x  1 HF2 ' ;  SPLIST( 5 )  x  'HF6'j  SPLIST(6)  x  >HF8' 

C  SPLIST ( 7 )  -  ' WHF ’ 
c 

c  CL1Q(1)  through  CLIQ(7)*  mass  fraction  of  species  in  the  liquid  phase 

c  corresponding  to  the  specie  number, 

c 

c  CLIQ(S)  «  Total  mass  of  Liquid  phase  (kg) 

c 

c  CVAP ( 1 )  through  CVAP(7)*  mass  fraction  of  species  in  the  vapor  phase 
c  corresponding  to  the  upecie:  number, 

c 

CVAP(8)  x  Total  mass  of  vapor  phase  (kg) 
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TABLE  2.1 


INDEX  TO  HP  AND  AIR  PARAMETERS  VARIED  IN  DIFFERENT  FIGURES 


Figure 

# 


■ 


■ 


Abscissa 

(x) 


Dilution 

Ratio 


. Jtlon 
Ratio 


Dilution 

Ratio 


Dilution 

Ratio 


Dilution 

Ratio 


Dilution 

Rutlo 


Dilution 

Ratio 


Dilution 
la  do 


Dilution 

Ratio 


Dilution 

Ratio 


Initial  Conditions  of 
HF  and  Air 


Cloud  Tempo retura 

19. 5 VC 

•  Mata  Fructlon  of  HF 
Condensed  Into 
Aorotol  Phntt 

19.  SVC 

Matt  Concnntratlon 
of  IIF  in  the  Aqueous 
Solution  (Aerosol) 

19.5VC 

Density 

19.  SVC 

Temperature 

19. 5 VC 

Mate  Fraction  of  IIF 
Condensed  Into 
Aerosol  Phase 

19 .  SVC 

Mass  Concsntratlon 
of  HF  In  the  Aqueous 
Solution  (Aerosol) 

19.  SVC 

Cloud  Density 

19. 5 VC 

Temperature 

1000  K 

Density 

1000  K 

Condition 

Temp 

»  RJI 

Sat. Vapor 

20V 

50 

Sat. Vapor 

20V 

SO 

Sat. Vapor 

20V 

50 

Sat. Vapor 


&  Liquid 


A  Liquid 


A  Liquid 


A  Liquid 


Vapor  at 
Amb.Pr. 


Vapor  at 
Ainb .  Pr . 


20°C  91 


20V  95 


20V  95 


I 


20*C  50 


20*C  50 


Figure  2.4  shows  the  variation  of  HF-air  mixture  temperature  as  a 
function  of  dilution  ratio  for  the  case  of  initial  HF  being  a 
saturated  vapor  at  ambient  pressure.  The  predictions  by  the  HFTHRM 
program  (modified  Schotte  model)  and  by  the  model  of  Clough,  et 
al.f  are  shown.  It  is  seen  that  the  predictions  by  both  models  are 
very  close. 

Note  that  for  dilution  ratios  less  than  about  10  the  temperature  of 
the  mixture  decreases  (substantially)  even  though  the  initial 
temperature  of  both  air  and  vapor  are  20°  C.  This  is  because  the 
saturated  HF  vapor  at  ambient  pressure  is  a  mixture  of  several 
associated  molecules  ("oligomers")  and  as  dilution  takes  place,  the 
molecules  dissociate  towards  a  monomer.  The  dissociation  reaction 
is  endothermic  resulting  in  the  lowering  of  temperature.  Above  a 
dilution  ratio  of  10,  almost  all  of  the  oligomers  have  dissociated 
and  additional  dilution  increases  the  temperature  of  the  mixture, 
ultimately  approaching  the  air  temperature  at  large  dilutions. 

Figure  2.5  shows  the  variation  of  the  mass  fraction  of  HF  that  is 
condensed  (into  the  liquid  aerosol)  with  variation  in  the  dilution 
ratio.  Condensation  occurs  because  of  lowering  of  temperature  as 
well  as  due  to  the  lowering  of  HF  partial  pressure  in  the  mixture. 
However,  at  large  dilution  ratios  the  condensed  HF  evaporates  and 
the  final  mixture  consists  essentially  of  HF  monomer  vapor  and  air. 

The  strength  of  the  HF  acid  (aqueous  in  the  form  of  aerosol)  with 
dilution  is  shown  in  Figure  2.6.  An  increase  in  the  dilution 
results  in  condensation  of  water  vapor  (due  to  lowering  of  mixture 
temperature)  and  dilutes  the  aqueous  solution.  Hence,  the  mass 
concentration  of  HF  in  the  liquid  aerosol  decreases.  At  very  high 
dilution  ratios  (>  10),  the  partial  pressure  of  HF  in  the  vapor 
phase  is  low  and  because  of  equilibrium  the  concentration  of  HF  in 
the  liquid  phase  also  is  low,  even  though  a  part  of  the  water  in 
the  aqueous  phase  may  start  to  re-evaporate . 

The  variation  of  cloud  density  with  dilution  ratio  is  shown  in 
Figure  2.7.  It  is  seen  that  for  dilution  ratio  below  50  for  air  at 
normal  ambient  temperature,  the  cloud  density  decreases 
continuously  as  dilution  increases  for  all  relative  humidities. 
For  air  at  or  below  50%  relative  humidity,  the  density  of  the 
mixture  is  always  greater  than  that  of  ambient  air  at  20°  C. 
However,  for  relative  humidities  greater  than  50%,  thero  is  a 
region  beyond  dilution  ratio  >  50  in  which  the  cloud  is  lighter 
than  air.  ThiB  is  because  of  the  residual  heat  of  reaction  between 
HF  and  high  humidity  in  the  air  at  higher  relative  humidities. 

The  variations  of  cloud  parameters  with  dilution  when  the  initial 
HF  contains  anhydrous  liquid  aerosol  particles  are  indicated  in 
Figures  2.8  through  2.11.  It  is  seen  that  the  behavior  of  all  of 
the  mixture  parameters  are  similar  to  that  described  earlier  except 
that  the  liquid  fractions  are  higher.  The  density  of  the  cloud  is 
higher  at  low  dilutions  (<  10)  compared  to  the  case  where  the 
initial  HF  was  all  in  vapor  form.  However,  the  differences  in 
cloud  density  for  dilutions  >  10  are  negligible. 
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FIGURE  2.4:  Variation  of  HF-Hunid  Air  Mixture  Temperature 
Dilution  Ratio 
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FIGURE  2.11:  Variation  of  Mixture  Density  with  Dilution  Ratio 


The  effect  of  high  temperature  HF  release  are  indicated  in  Figures 
2.12  and  2.13.  Figure  2.12  shows  the  variation  of  mixture 
temperature  with  dilution  for  a  1000  K  HF  release  at  ambient 
pressure.  At  this  high  temperature,  almost  all  of  the  HF  released 
is  in  the  form  of  monomer  and  the  dilution  further  promotes  the 
monomer  phase.  The  mixture  temperature  decreases  continuously. 
However,  it  is  noticed  that  the  temperature  decrease  is  not  the 
same  as  would  be  the  ease  when  two  ideal  gases  at  different 
temperatures  are  mixed.  This  is  because  higher  molecular  weight 
oligomers  are  formed  during  dilution  of  high  temperature  HF.  The 
association  reaction  is  exothermic,  thus  the  rate  of  temperature 
drop  is  lessened  compared  to  that  in  the  case  of  ideal  gases. 
However,  at  large  dilution  ratios  the  partial  pressure  of  HF 
decreases  rapidly  resulting  in  the  promotion  of  dissociation 
eactions  to  monomer.  Hence,  the  curvature  of  the  temperature  vs. 
dilution  flattens. 


The  density  variation  ie  shown  in  Figure  2.13.  The  initial  density 
of  the  high  temperature  HF  monomer  is  quite  low  (compared  to 
ambient  air  density).  In  fact,  at  1000  X,  HF  acts  very  much  like 
an  ideal  gas.  Dilution  of  this  with  ambient  temperature  air 
results  in  a  monotonic  increase  in  mixture  density.  Note,  however, 
the  ohangee  in  the  curvature  of  the  density  vs.  dilution  curve. 
This  is  due  to  the  association  and  dissouation  reactions  as 
explained  in  the  previous  paragraph. 
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Dilution  Ratio  =  (Kg  Air/Kg  HF) 

FIRJRE  2.13:  VARIATION  OF  MIXTURE  DENSITY  KITH  AIR  DILUTION  FOR  HIGH 
•EKPERATURE  HF  RELEASE 
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CHAPTER  3 


REVISIONS,  MODIFICATIONS,  ENHANCEMENTS 
AND  CORRECTIONS  TO  ADAM  CODE 

Several  modifications  were  made  to  routiner  in  ADAM  to  make  the 

program  more  robust  and  applicable  over  a  greater  range  of 

oarameters.  In  this  chapter,  we  discuss  the  theoretical  and 

physical  bases  of  these  modifications,  and  the  changes  in  detail. 

The  principal  changes  made  include: 

1.  Calculation  of  atmospheric  wind  velocity  profile  for  flow  over 
very  large  aerodynamic  roughnesses  and  the  impact  on 
dispersion  of  a  heavy  gas  Cxoud  whose  vertical  depth  is 
smaller  than  the  mean  height  of  the  roughness  elements. 

2.  Determination  of  atmoBpheL '  wind  friction  velocity  under 
specified  atmospheric  condi\  .  s. 

3.  Modification  of  algorithms  to  calculate  certain  diffusion 
property  values. 

4.  Enhancement  of  the  calculation  procedure  to  determine  the 
atmospheric  stability  value. 

5.  corrections  of  the  coda  to  remove  certain  incorrect  equations 
or  parameter  value  estimations. 

Each  of  these  is  discussed  below. 


3.1  Correction  of  Large  Aerodynamic  Roughness  In  The  Dispersion 
Path 

The  roughness  of  the  ground  over  which  wind  blows  influences  the 
distribution  of  time  averaged  wind  profile  variation  with  height. 
In  the  report  by  Raj  and  Morris  (1987),  equations  were  presented 
for  average  wind  velocity  distribution  with  height  under  different 
types  of  atmospheric  stabilities.  These  equations  contain  a  length 
scale  z  representing  the  "aerodynamic  roughness"  of  the  path  over 
which  the  wind  is  flowing.  The  correlations  presented  in  the 
referenced  report  are  accurate  only  for  the  case  where  the 
magnitude  of  the  aerodynamic  roughness  (z0)  is  small  compared  to 
the  vertical  depth  of  the  dispersing  vapor  cloud. 

Occasionally,  the  aerodynamic  roughness  of  the  path  may  be  large 
compared  to  the  depth  of  the  cloud  (example:  a  large  forest  of 
tall  trees  or  tall  buildings) .  The  dispersion  of  a  cloud  under 
these  conditions  is  very  complex  and  will  depend  on  not  only  the 
height  of  the  roughness  elements,  but  also  the  fractional 
horizontal  area  occupied  by  the  solid  elements,  the  mean  distance 
between  the  solid  elements  as  a  fraction  of  the  height  of  solid 
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elements,  etc.  The  cloud  may  be  aplit  up,  caught  up  in  the  wake 
cavities  of  the  tall  elements,  be  mixed  with  air  more  due  to 
higher-than-ambient  turbulence,  etc.  However,  the  mean  velocity  of 
the  cloud  within  the  tall  elements  will  be  lower  because  the  wind 
speed  within  the  interstices  of  the  elements  will  be  lower  than  in 
the  upstream  of  the  high  aerodynamic  roughness  region.  These 
phenomena  are  extremely  difficult  to  model.  Only  an  approximation 
is  Included  in  the  modified  version  of  ADAM. 

In  the  previous  version  of  ADAM,  the  occurrence  of  a  large 
aerodynamic  roughness  region  (simulation  in  the  code)  would  have 
resulted  in  the  calculation  of  negative  wind  speeds  (for  z  <  z0)  - 
a  physically  unacceptable  solution.  We  have  modified  the  wind 
velocity  profile  (described  below)  to  take  into  account  large 
aerodynamic  roughness  heights  and  atmospheric  stabilities. 


3.1.1  Assumptions 
We  assume  that 

1.  The  mean  wind  opoed  varies  linearly  with  height  above  ground 
in  the  region  of  large  aerodynamic  roughness,  up  to  a 
"critical  height"  (zc). 

2.  The  magnitude  of  the  critical  height  is  dependent  on  the 
magnitude  of  the  aerodynamic  roughness  z0  and  not  on  the 
atmospheric  stability. 

3.  The  wind  velocity  distribution  above  the  critical  height  is 
the  same  as  in  undisturbed  flow. 

The  value  of  this  critical  height  (ze)  is  determined  by  forcing  the 
continuity  of  both  the  wind  velocity  value  as  well  as  its  slope  at 
z  -  z  from  the  velocity  profiles  above  z  -  z  and  below  z  -  z  . 
This  is  illustrated  schematically  in  Figure  3.1  and  discussed  below 
in  more  detail.  All  three  typea  of  atmospheric  stabilities  are 
discussed. 


a)  Neutral  Stability  Atmosphere  (Monin-Obukov  Length.  L  -  0) 

The  wind  speed  distribution  with  height  ie  assumed  to  be  given  by: 


U„(z) 


(u./k)  ln(z/z0) 
(u./k)  (z/ze) 


for  z  >»  zc  (3.1a) 

for  z  <»  zc  (3.1b) 


where  ze  is  a  yet-to-be  determined  height  whose  value  depends  on 
the  value  of  the  aerodynamic  roughness  height  z  .  It  is  noted  that 
the  form  of  Equation  3.1a  is  identical  to  the  distribution  under 
neutral  stability  when  z0  is  small.  Equation  3.1b  represents  the 
linear  distribution  assumption  that  we  have  made. 
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Wind  Speed  [m/s] 

FIGURE  3.1:  Wind  Velocity  Profile 


To  calculate  the  value  of  zc,  we  Impose  the  following  conditions: 
1.  Uw(z)  i®  continuous  at  z  »  zc. 


2.  d  Uu(z)  is  also  continuous  at  z  -  z_. 

dz 

It  can  be  shown  that  with  the  above  conditions  and  Equations  3.1a 
and  3.1b,  we  get 


ze  -  s  (3.2) 

whore  e  is  the  base  of  the  natural  logarithm. 

The  average  wind  speed  over  a  specified  height  z  is  given  by 

z 

uTTZT  -  (l/Z)  }  U„(B)  dZ  (3.3) 

Z«0 

Substituting  Equations  3.1a  and  3,1b  in  Equation  3.3  and  noting 
Equation  3.2,  we  get 

_ _  <u*/k)  [0.5  (*/■„)]  -  0.5  U  (z)  for  z  <-  z  (3.4a) 

- 

(u./k)  [0.5*  (z/zc)+ln(z/z0)  -  1]  for  z  >  ze  (3.4b) 

(b)  Stable  Atmosphere  (L  >  0) 

The  wind  speed  distribution  in  a  stable  atmosphere  again  is  assumed 
to  be  given  by: 


U„<z)- 
where , 
a 

where, 

L 


(u./k)  [ln(z/z0)  +  5.2  a) 

for  z  >  z0 

(3.5a) 

UN(ze)  (z/ze) 

for  z  £  ze 

(3.5b) 

f  (*/b) 

m  S 

for  z  <  L 

(3.6a) 

U 

for  z  >■  L 

(3.6b) 

Monin-Obukhov  turbulence  length  scale  in  the  atmosphere 


where  Uw(ze)  is  the  wind  speed  at  a  height  z  .  The  values  of  zc  and 
Uy(ze)  are  found  by  imposing  the  continuity  conditions  for  both 
velocity  function  and  its  derivative  with  height  at  z  -  z  .  There 
are  two  cases  to  consider;  namely,  zc  <  L  and  zc  >  L.  Applying  the 
continuity  conditions  and  using  the  velocity  distribution  Equations 
3.5a  and  3.5b,  it  can  be  shown  that  in  both  cases  z  «  ez0  (as 
indicated  in  Equation  3.2). 
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(3.7a) 


Also,  the  wind  speed  Uc  at  z  - 


*  Uc 


u. 


1  +  5.2  «c 


ze  ia  given  by: 


where, 

'  (*e/L)  for  ze  <  L 

«e 

„  1  for  ze  >-  L 

The  average  wind  speed  ia  calculated  using  Equations  3.3 
3.5b  with  conditions  of  Equations  3.2  and  3.6. 

i.e. , 

L3L  .  — L-  f  -JSJUil _ a, 

U.  z  J  u. 

Z«0 


The  above  equation  reduces  to 
Case  1  (z  <  ze) 

0w(z)  -  0.5  (z/ze)  ue  -  0.5  UM(Z) 

Case  2  (b  >  bg) 

Equation  (3.8)  can  be  written  as 


zc  z 


(3.7b) 

(3.7C) 
,  3.5a  and 

(3.8) 

(3.9) 

(3.10) 
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With  ze  -  ezf)  and  substituting  Equations  3.5b  and  3.5a  respectively 
in  the  integrals  on  the  RHS  of  the  above  equation  and  simplifying, 
we  get 


k  u 
u* 


k  Uc 

0.5  (ze/z)  -  +  [ln(z/z0)  -  l]  + 

u* 


(3.11) 


Again  depending  on  the  relative  values  of  a,  L  and  z  ,  there  are 
three  oases  to  consider.  The  results  are  summarized  below. 


^  U“if) 

u* 


k  ue 

0.5  (zc/z)  - -  "■■■  +  [{n  (z\z0)  -  1]  +  5 . 2x 

u. 


where  x 


x 


X 


[  1  -  zc/z]  for  z  >  zc  >  L  (3.12a) 

[  (z/L)2  -  (zc/L)2  ] 

- — - - -  for  L  >  z  >  z  (3.12b) 

2  (z/L) 


[  2  (Z/L)  -  (Zc/L)2  -  1] 
2  (Z/L) 


for  z  >  L  >  Z0 


(3.12c) 


(o)  Unstable  Atmosphere  (L  <  0) 


The  veloaity  distribution  in  an  unstable  atmosphere  for  small 
values  of  z0  is  given  by 


k  Uw(z)  -1  -1  (x  -  l)  (x0  +  1) 

-  -  2  [tan  (x)  -  tan  (x0)]  +  In  -  (3.13a) 

u*  (x  +  1)  (x0  -  1) 

where, 

x  -  [  1  -  15  z/L]0-25  and  x0  -  [1  -  15  z0/L]°‘25  (3.13b) 


We  assume  that  within  reasonable  limits  for  z0,  the  Equations  in 
3.13  will  apply  to  the  region  above  a  value  of  z  -  ze,  when  a 
surface  with  large  values  of  aerodynamic  roughness  is  encountered. 
We  assume  that  in  the  case  of  unstable  weather,  also  we  have  the 


following  relationship: 

*e  "  e  z0 

(3.14) 

xe  -  [1  -  15  Ze/L]°'2s 

(3.15) 
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U«  - 


U.  ,  .  (xe  “  1)  (x0  +  1) 

—  2  [tan-1  x.  -  tan*1  x0]  +  tn  -  — — —  — — 

k  ( Xq  -  l)  (xe  +  l) 


and  the  velocity  distribution  is 

(z/zj  for  z  <-  z. 


‘  {-i 


(3.16) 


(3.17a) 


£*-  2  [tan-1x  -  tan-’x,,]  +  in  [*a  $ 

for  z  >  zB  (3.17b) 

Tha  average  velooity  ovar  a  height  z  can  be  determined  by  the 
equation: 


Tuz> 


(i/z) 


z 

/  “■ 

J  Z«0 


(z)  dz 


(3.18) 


Since  it  is  assumed  that  the  velocity  distribution  is  linear  for  z 
<  ze,  we  can  show  that 


UN(z)  -  0.5  (z/zc)  Uc  for  z  <-  ze  (3.19a) 

However,  for  z  >  ze,  the  integration  in  Equation  3.18  has  to  be 
performed  in  two  regions.  In  the  region  z  >  z#,  the  velooity 
distribution  is  given  by  Equation  3.17b  -  a  very  complicated 
expression  for  integration.  We,  therefore,  make  a  simplifying 
assumption  that  for  obtaining  mean  velocity  above  z  *»  ze,  the 
following  approximate  form  of  the  velocity  distribution  is  assumed 


UH(z)  “  Uc  (z/zc)<1/7)  for  z  >-  zc  (3.20) 


Using  Equations  3.20  and  3.17a  in  Equation  3.18,  we  gat 


Uw(z)  -  Uc  (zc/z) 


[  7  (z/zc)8'7  -  3  ] 


for  z  >  ze  (3.19b) 


or 


UM(Z)  -  (7/8)  UM(Z)  -  (3/8)  (ze/z)  uc  for  Z  >  Ze  (3.19C) 
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It  should  be  noted  that  the  value  of  L,  the  Mon  in-Obukhov 
turbulence  length  scale  in  the  atmosphere  and  the  friction  velocity 
vi*  are  calculated  using  the  z0  value.  Therefore,  the  use  of  the 
above  equations  may  not  be  appropriate  if  the  meteorological 
information  is  taken  at  one  place  and  the  cloud  dispersion  is  at  a 
different  place  where  the  aerodynamic  roughness  is  substantially 
different. 


3.2  Calculation  of  Atmospheric  Turbulence  Friction  Velocity 

The  program  code  in  ADAM  provides  the  user  the  flexibility  of 
either  inputting  the  value  of  the  stability  of  the  atmosphere  or 
having  the  program  calculate  the  stability  based  on  the  location  of 
the  place,  time  of  day,  and  other  local  conditions,  in  the  case 
the  stability  value  is  input,  it  is  necessary  to  calculate  the 
characteristics  of  atmospheric  turbulence.  This  section  gives  the 
details  of  the  calculation  of  wind  friction  velocity  and  Monin- 
Obukhov  turbulence  length  scale. 


3.2.1  Assumptions 

It  is  assumed  in  these  calculations  that 

a)  the  atmospheric  stability  value  (SP)  and  the  aerodynamic 
roughness  z0  are  given 

b)  the  mean  wind  speed  Uw  is  specified  at  a  height  zM 

c)  zM  »  ez0 

where  "e"  is  the  base  of  the  natural  logarithm. 

The  objective  of  the  calculations  is  to  determines  the  values  of  L, 
the  Monin-Obukhov  length  and  u,  the  wind  friction  velocity. 

3.2.2  Analysis 

Following  Kunkel  (1986),  we  write  the  equation  for  the  atmospheric 


stability  as 

SP  -  A  +  B  log,0  (100  Zq) 

(3.21) 

where,  A  ■  3.5  +  21.67/L 

(3.22a) 

B  -  0.48 

for  abs(l/L)  >  0.015 

(3.22b) 

B  -  4  3.63  (abS  (1/L)1,08 

for  abs  (1/L)  <  0.015 

(3.22C) 

B  =  -B  for 

L  <  0 

(3.23) 

3-B 


where  SP  is  the  given  value  of  atmospheric  stability,  z.  the 
aerodynamic  roughneos,  and  L,  the  Monin-Obukhov  atmospheric 
turbulence  length  scale.  Given  the  value  of  SP  and  z0,  the  value 
of  L  is  determined  by  an  iterative  search  procedure. 

The  wind  friction  velocity  u*  is  determined  knowing  the  L  value  and 
wind  velocity  at  a  specified  height.  This  procedure  is  indicated 
below 


Stable  Atmosphere  (L  >  0) 


X..U»U1 

u. 


■  in  (z/z0)  +  5.2  # 
z/L  for  z  <  L 


1  for  z  2  L 

Neutral  Atmosphere  [abs(L)  -<*) 

k  U(z) 


u. 


«n  (z/z0) 


Unstable  Atmosphere  (L  <  0) 

rk  u  (■)!  I  (X  -  1)  (x0  +  1) 

-  2  [tan'1  x  -tan*1  x0]  +  fn 


u. 


(x0  -  1)  (X  +  1) 


(3.24a) 

(3.24b) 

(3.24c) 


(3.25) 


(3.26a) 


whore 

X  -  (1-15  z/L)0-”  (3.26b) 

x0  -  (1  -  15  20/L)°-”  (3.26C) 


The  value  of  u*  ie  determined  (by  a  process  of  iteration,  if 
necessary)  using  one  of  the  above  velocity  distributions 
appropriate  for  the  stability  value  SP.  These  calculations  have 
been  coded  into  a  subroutine  called  FRICTVEL. 
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3.3  correlation  For  Molaoular  Diffusion  coefficient 

In  calculating  the  value  of  mass  transfer  coefficient  to  estimate 
the  rate  of  evaporation  from  a  pool  of  liquid  on  the  ground,  it  is 
necessary  to  knew  the  value  of  the  molecular  diffusion  coefficient 
for  the  chemical  vapor  in  air  [see  Equation  2.5.9  in  the  report  by 
Raj  and  Morris  (1987)].  The  calculation  procedure  for  determining 
the  diffusion  coefficient  has  been  modified  to  take  into  account 
the  more  recent  correlations. 

The  correlation  used  is  the  Fuller,  et  al.,  correlation  (Reid, 
Prausnitz  and  Poling,  1987)  and  is  given  by 


1.43  x  10‘n  T1,75 


Dab 


P  MAI10J  [  Al/3  +  B,/J  ]J 


(3.25) 


where 

Dab  -  Diffusion  coefficient  for  Species  A  in  Species  B  (m5/*) 


T  -  Temperature  (X) 

P  ■  Pressure  (atm) 

Mah  -  Harmonic  mean  value  of  molecular  weight  of  Speaias  A 
and  Species  B  »  2/[l/MA  +  1/MB) 

A  «  Sum  of  atomic  diffusion  volumes  of  the  atomic  component  of 
Specie  A 

B  -  Sum  of  atomic  diffusion  volumes  of  the  atomic  component  of 
Specie  B 


Table  3.1  shows  atomic  diffusion  volumes  for  various  elements. 
Also  indicated  in  the  table  are  the  volume  contraction  for 
different  types  of  molecules  (aromatic  ring,  Heterocyclic  ring, 
etc. )  . 


3.4  Other  Corrections  To  and  Enhancements  in  ADAM 
3.4.1  corrections 

Several  corrections  were  made  to  ADAM  to  briny  the  code  to  coincide 
with  theoretical  analysis.  These  include: 
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i)  correcting  an  error  in  the  program  code  simulating  the 
mixing  of  nitrogen  tetroxide  and  air  ("NOXTHRM").  The 
error  in  the  determination  of  the  equilibrium  constant 
value  in  the  subroutine  EQUBCQNS  was  corrected. 

ii)  modifying  an  incorrect  exponent  value  and  the  density 
calculation  in  the  compressed  gas  release  subroutine 
GASREL. 


TABLE  3.1 

ATONIC  Dirruiiou  VOUMtl 

Atomic  and  Structural  Mffuaton  Voltm 

Ineramanti 

c 

- 1179 

7 

16.7 

H 

2.31 

Cl 

21.0 

0 

6.11 

Ir 

21.9 

N 

A. 36 

1 

29,6 

Aromatic  Ring 

•16.3 

1 

22.9 

murooyctlc  Ring 

•18.3 

Dlffutlon  Vo  Luma  ol  limpla  Kalaeulaa 

Ha 

zT57 

CO 

Ha 

3.98 

co2 

M 

26.9 

1 

16.2 

33.9 

Kr 

26.3 

NR, 

H.O 

20.7 

Hi 

37.7 

13.1 

He 

6.12 

"6 

71.3 

oj 

6.66 

Cl} 

36.6 

N! 

18.3 

Br! 

69.0 

A?r 

16.3 

19.7 

4 

61.8 

jgourctn  Raid,  Pruunnlti  1  Holing  <19671 

3.4.2  Enhancements 

i)  Leak  From  Pipes 

ADAM  was  improved  to  give  proper  liquid  release  rates  for  a  pipe 
leak  situation  when  the  pipe  pressure  was  ambient  but  the  hole  was 
on  the  pipe  wall  wetted  by  the  liquid.  The  release  rate  is 
conservatively  estimated  by  assuming  the  hole  to  be  at  the  lowest 
position  on  the  pipe  wall  (i.e,,  largest  hydrostatic  head  for 
liquid  leak) . 

ii)  Atmospheric  Parameters 

The  atmospheric  parameter  calculation  subroutine  has  been  modified 
to  taka  into  aacount  the  possible  differences  in  the  value  of  the 
aerodynamic  roughnoes  at  the  meteorological  tower  site  and  at  the 
place  of  vapor  dispersion.  The  atmospheric  stability  and  other 
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parameters  are  calculated  on  the  basis  of  z0MT,  the  aerodynamic 
roughness  at  wind  measuring  site  (or  the  meteorological  tower) . 
The  Gaussian  dispersion  parameters  o  and  at  are  calculated  using 
the  value  z0,  the  aerodynamic  roughness  at  the  dispersion  location. 
Also,  the  correlations  with  which  the  values  of  oy  and  o  are 
calculated  as  function  of  downwind  distance  and  atmospheric 
stability  have  been  modified  slightly  to  give  a  better  fit  to  the 
Pasquill-Gifford  curves. 

iii)  Liquid  Fraction  Entrained  in  Vapor  Cloud 

The  release  of  a  compressed  liquefied  gas  into  the  atmosphere 
results  in  the  flash  vaporization  of  a  part  of  the  released  mass. 
A  mass  fraation  f  manifests  as  saturated  vapor  and  the  remainder 
as  saturated  liquid  of  mass  fraction  fL  (fL  +  fv  «  1) .  However,  a 
part  of  this  liquid  mass  can  get  entrained  into  the  dispersing 
vapor  cloud/plume.  In  fact,  the  evidence  from  field  tests  (the 
Burro  series  of  pressurized  ammonia  releases,  the  Goldfish  series 
of  hydrogen  fluoride  tests)  seems  to  indicate  that  a  large  part  of 
the  liquid  gets  entrained  as  aerosol  into  the  dispersing  cloud. 

In  order  to  take  into  account  different  amounts  of  liquid 
entrainment  into  the  initially  formed  vapor,  a  parameter  <|>  was 
introduced  in  the  model  in  ADAM  (see  Section  2. 6. 2. 2,  Equation 
2.6.5,  p . 2-41 ,  of  the  report  by  Raj  and  Morris,  1987).  This  factor 
represents  the  (mass)  fraction  of  the  saturated  liquid  mass  formed 
after  flash  which  gets  entrained  into  the  dispersing  vapor.  The 
value  of  $  can  be  in  the  range  0  5  <)>  <  1. 

In  the  original  ADAM  code,  this  parameter  value  was  a  user  input. 
Unfortunately,  this  led  to  confusion.  Hence,  we  have  removed  this 
as  an  input  parameter  and  included  that  as  a  data  item  in  the 
ohemical  property  database.  The  default  value  is  net  to  ■  1 
(i.e.,  all  liquid  ie  entrained  into  the  vapor  phase).  We  do 
realize  that  the  value  of  4  is  not  dependant  solely  on  the  chemical 
property  but  also  on  the  conditions,  shape,  orientation,  etc.  of 
the  hole  and  whether  the  jet  impacts  the  ground  or  not.  However, 
by  including  the  value  as  a  data  item  in  the  chemical  property 
database,  it  provides  the  flexibility  for  an  informed  user  to 
change  the  value  if  deemed  necessary.  Other  users  may  not  want  to 
even  bother  with  the  parameter  (it  is  transparent  to  them), 


3.4.3  Panel/Menu  Driven  Inputs 

A  substantial  improvement  has  been  implemented  to  ADAM.  This 
involves  the  complete  revamping  of  the  data  Input  system.  The 
previous  system  was  based  on  inputting  data  through  an  editor 
routine.  This  was  complex  and  sometimes  confusing. 

The  revised  data  input  system  is  based  on  the  use  of  menu-driven 
panels.  Also  provided  are  "help"  menus  to  define  the  various  terms 
in  the  input.  More  detailB  of  these  input  panels,  and  other 
functions  to  run  ADAM  are  indicated  in  Volume  II  of  this  report. 
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CHAPTER  4 


COMPARISON  OF  EXPERIMENTAL  DATA  AND  MODEL  RESULTS 
FOR  HYDROGEN  FLUORIDE 


In  this  chapter,  we  compare  the  data  from  the  Goldfish  series  of  HF 
release  tests  with  the  calculated  results  obtained  by  exercising 
the  ADAM  code  simulating  HF  dispersion  under  the  test  conditions. 
First,  we  briefly  describe  the  test  conditions  and  other  pertinent 
data.  Than,  the  model  results  are  compared  with  the  test  results. 
A  discussion  on  the  results  is  also  given. 

4.1  Brief  Description  of  Field  Tests 

A  series  of  six  field  tests  (called  the  Goldfish  Series)  involving 
the  release  of  anhydrous  hydrofluoric  acid  is  reported  by  Blewitt, 
et  al.  (1987).  The  tests  are  also  described  in  a  recent  report  by 
Hanna,  et  al.  (1989).  Only  three  tests  in  the  series  were 
primarily  HF  vapor  dispersion  teats. 

The  tests  consisted  of  releasing  liquid  anhydrous  HF  from  a  tank 
pressurized  with  nitrogen.  The  temperature  of  the  liquid  HF  in  the 
spill  pipe  was  maintained  close  to  40  nC  with  an  electric  heater. 
The  spill  pipe  was  horizontal,  at  1  m  above  ground  and  pointed  in 
the  downwind  direction.  Soma  of  the  data  on  the  test  facility  and 
test  conditions  are  given  in  Table  4.1,  The  HF  release  rata  was 
measured  by  weighing  the  HF  trailer  continuously  on  a  load  cell  and 
calculating  the  slope  of  the  weight  vs.  time  curve.  The  mean  flow 
rate  thus  measured  is  indicated  in  Table  4.2. 

It  is  reported  by  Blewitt,  et  al.  (1987)  that  two  types  of  HF 
concentration  measuring  instruments  were  used;  only  one  type,  the 
Integrated  Filter  Sampler  (IFS),  worked  properly.  Hence,  only  the 
concentration  data  from  this  type  of  instrument  were  used  in  our 
analyeis.  Temperature  data  wero  also  measured  at  several  locations 
downwind.  No  instruments  were  ueod  to  measure  the  aerosol  content 
or  the  quality  of  HF  in  the  plume  (i.e.,  concentration  of  various 
oligomers) ,  Also,  the  plume  temperature  data  and  concentration 
data  were  not  measured  at  the  same  locations.  More  detailed 
information  on  the  instrumentation,  their  location,  accuracy  of 
measurements  can  be  obtained  from  the  paper  by  Blewitt,  et  al. 
(1987),  and  in  the  report  by  Hanna,  et  al.  (1969). 

In  the  following  section,  we  compare  the  experimentally  measured 
parameter  valu&s  with  ADAM  pradictiona. 

4.2  Rslssss  Rata  and  Flash 

Several  parameter  values  calculated  using  ADAM  for  the  three  test 
conditions  are  indicated  in  Table  4.2.  This  table  gives  the  valuos 
for  two  sets  of  parameters,  namely,  the  thermodynamic  parameters 
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and  the  flow  parameters  of  the  jet  just  after  the  exit  from  the 
orifice  plate.  There  are  only  limited  data  from  the  field  tests 
with  which  these  calculated  values  can  be  compared. 

Table  4.2  shows  the  HF  temperature  and  pressure  upstream  of  the 
orifice''*.  The  saturation  pressures  of  HF  corresponding  to  the 
temperatures  of  HF  upstream  of  orifice  are  indicated  in  this  table. 
The  saturation  pressure  corresponding  to  the  HF  temperature 
upstream  of  the  nozzle  is  lower  than  the  pressure  maintained  in  the 
tank  with  nitrogen.  Hence,  the  HF  liquid  is  in  a  compressed  liquid 
state.  It  can  be  shown,  however,  that  because  of  the  relatively 
low  enthalpy  of  liquid  compression  the  specific  enthalpy  of  the 
liquid  HF  upstream  of  the  orifice  is  almost  the  same  as  the 
specific  enthalpy  of  the  saturated  liquid  corresponding  to  the 
upstream  temperature.  The  details  of  this  analysis  are  indicated 
in  Appendix  A. 

When  this  compressed  liquid  is  released  to  the  atmosphere,  a  part 
of  the  liquid  flashes  to  form  vapor.  The  vapor  and  liquid  formed 
arc  in  the  respective  saturated  states  corresponding  to  the 
atmospheric  pressure.  The  calculated  mass  fraction  of  the  released 
HF  which  flashes  to  vapor  is  indicated  in  this  table.  Also  shown 
are  the  temperature  (saturation)  of  the  vapor  and  liquid 
immediately  downstream  of  the  orifice. 

The  vapor  formed  efter  the  flash  will  be  a  mixture  of  various 
polymern  of  HF.  Therefore,  the  molecular  weight  of  the  vapor  will 
be  considerably  different  from  that  of  the  monomer.  The  effective 
molecular  weight  of  the  saturated  vapor  is  calculated  by  the 
correlation  (see  Figure  2.2) 

H  vap  (T)  “  156 . 67  -  0.3  T  (4.1) 

where 

H  “  (T)  >=  apparent  molecular  weight  of  the  saturated 

Vap  vapor  (kg/kmole) 

T  »  saturation  temperature  (K) 

These  apparent  molecular  weights  of  the  vapor  at  a  position 
immediately  downstream  of  the  orifice  are  indicated  in  Table  4.2 
for  the  conditions  of  the  tests.  The  density  of  vapor,  calculated 
using  these  molecular  weight  values,  ambient  pressure,  saturation 
temperature  and  the  perfect  gas  assumption  is  also  indicated  in 
Table  4.2.  Unfortunately,  test  data  are  not  available  to  compare 
with  these  calculated  values. 

**  There  is  a  discrepancy  between  the  values  for  this  parameter 
reported  by  Blewitt,  et  al,  (1987)  and  by  Hanna,  et  al. 
(1989).  We  have  used  the  values  published  by  Hanna,  et  al. 
(1989)  . 


The  second  part  of  Table  4.2  contains  the  calculated  flow  rate 
values.  The  flow  rates  are  calculated  in  ADAM  usinq  the  Bernoulli 
equation  and  assuming  a  coefficient  of  discharge  (C0)  value  of  0.6. 
This  coefficient  multiplies  the  Bernoulli  exit  velocity  to  give  a 
reduced  velocity.  The  mass  flow  rate  i3  then  determined  by  the 
product  of  orifice  nominal  area,  reduced  velocity  and  liquid 
(saturated)  density  at  ambient  pressure.  These  calculated  values 
for  the  reduced  velocity,  liquid  density  and  the  mass  flow  rate  are 
indicated  in  the  table. 

The  flashing  process  results  in  the  formation  of  a  two-phase  jet. 
Because  of  the  reduction  in  the  mean  density  (compared  to  thu 
liquid  density)  of  the  flow  consequent  to  the  flashing  process,  the 
jet  size  expands  vsry  close  to  the  exit  section  of  the  orifica. 
The  calculated  values  of  ths  mean  density  of  the  two-phase  flow  are 
indicated  in  the  table.  These  values  are  determined  based  on  the 
assumption  that  all  of  ths  liquid  fraction  (in  the  post  flash 
stream)  is  entrained  as  aerosols  into  the  jet.  Also  indicated  in 
the  table  are  the  calculated  cross  section  of  the  two-phase  flow  in 
the  jet  very  close  to  the  orifice. 

Comparison  of  Source  Condition  Data  with  Calculated  Results 

No  sxpsrimental  data  seemed  to  have  been  measured  very  close  to  the 
exit  section  of  the  orifice.  Hence,  the  calculated  results  for  jet 
velocity,  quality  of  the  two  phase  mixture,  temperature,  etc.  at 
the  exit  section  cannot  be  compared  with  test  data.  The  calculated 
value  for  the  initial  diameter  of  the  two  phase  jet  can  be  compared 
to  the  jet  diameter  from  photographs  with  some  uncertainty  on  the 
photographic  scales  (due  to  perspective  and  mirage  effects) .  The 
only  source  parameter  that  can  be  compared  is  the  mass  flovr  rats, 
sines  this  can  be  inferred  from  the  experimental  data  on  the  HF 
trailer  weight  vs.  time. 

The  calculated  mass  flow  rate  values  and  the  experimental  values 
for  tests  2  and  3  agree  very  closely;  however,  for  test  1,  ADAM 
predicts  about  9%  higher  values.  It  is  noted  that  there  ace 
discrepancies  in  the  reported  values  for  mass  flow  rates  (test 
data)  by  different  researchers.  For  example,  for  test  3,  Hanna,  et 
al.  (1989)  report  a  value  of  10.27  kg/s.  The  value  determined  by 
us  using  the  graph  of  weight  vs,  time  for  this  test  published  by 
Blewitt,  et  al.  (1987)  Ib  10.9  kg/s.  Blewitt,  et  al.,  gives  the 
mass  flow  rate  values  in  a  table  in  units  of  gpm  without  indicating 
what  values  of  liquid  density  were  used  in  converting  the  raw  data 
(kg/s)  to  gpm.  Assuming  the  liquid  density  (for  test  3,  see  Table 
4.2)  to  be  964.44  kg/m3,  we  can  convert  the  gpm  data  of  Blewitt,  et 
al.  to  kg/s;  this  calculation  leads  to  12.54  kg/s l  However,  the 
ADAM  calculated  value  and  the  value  measured  off  the  graph  for  test 
3  agree  very  closely.  Therefore,  we  feel  that  the  algorithm  for 
source  condition  calculations  in  ADAM  is  quite  accurate. 

Comparison  of  other  calculated  source  condition  parameter  values 
with  test  data  are  not  possible  because  none  of  the  other 
parameters  were  measured. 
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4.3  Dispersion  Results 


We  consider  the  comparison  of  the  results  in  two  distinct 
f low/dispersion  regimes.  The  first  is  the  jet  flow  regime  where 
the  flow  velocities  in  the  plume  are  considerably  higher  than  the 
velocity  of  the  prevailing  wind.  The  second  regime,  combines  heavy 
gas  dispersion  and  passive  dispersion. 

4.3.1  Jet  Dispersion  Ragims 

It  is  seen  from  Table  4.2  that  the  calculated  jet  velocities  are 
about  23  m/s  whereas  the  wind  speeds  are  in  the  range  4.2  m/s  to 
5.6  m/s.  Because  of  the  relatively  high  speed  with  which  the  jet 
flows  out  of  the  orifice,  air  is  entrained  into  the  jet.  In  ADAM, 
the  jet  dilution  is  calculated.  The  calculations  are  terminated 
when  the  mean  velocity  in  the  jot  is  within  5%  of  the  wind  speed  at 
the  level  of  the  top  of  the  jet.  The  dispersion  is  then  modeled 
using  the  heavy  gas  model  until  the  density  of  the  cloud  is  very 
close  to  that  of  air.  At  that  time,  the  dispersion  modeling  is  done 
using  the  Gaussian  model. 

In  Figures  4.1o,  4.1b  and  4.1c,  the  calculated  values  for  variation 
of  jet  width,  peak  HF  concentration’,  mean  velocity  in  the  jet,  and 
the  centerline  plume  temperature  are  plotted  as  functions  of 
downwind  distance  from  the  orifice.  The  Figure  4.1a  results  are 
for  conditions  of  Goldfish  Test  #1.  similarly,  Figure  4.1b  and 
Figure  4.1c  represent  the  plots  for  Test  #2  and  Test  #3, 
respectively.  An  aerial  photograph  of  the  plume  was  available  for 
Test  #2.  The  data  on  plurr.o  width  (obtained  by  measuring  the 
visible  cloud  width)  are  plotted  in  Figure  4.1b.  The  concentration 
at  which  the  plume  becomes  invisible  is  not  known.  This  will 
depend  on  the  relative  humidity  and  (perhaps)  on  the  optical 
properties  of  HF  vapor.  The  calculated  plume  width  plotted  in 
Figure  4.1b  is  the  equivalent,  "box"  width  in  the  jet  region. 

Data  on  plume  temperature  measured  during  the  tests  are  also 
plotted  on  these  figures. 


Discussions  on  Jet  Regime  Results 

The  general  trend  for  the  variation  of  each  of  the  parameters  with 
distance  for  all  three  tests  are  similar.  The  temperature 
initially  decreases  from  the  saturation  temperature  corresponding 
to  the  ambient  pressure  up  to  about  15  m  distance  and  then 
increases.  The  initial  dip  in  the  temperature  is  due  to  the 


2  Concentration  is  given  in  units  of  kg  of  total  Hl;  per  unit  volume  (nr1)  of  space.  To 
convert  this  into  molar  or  ppm  units,  a  chemical  molecular  weight  is  needed.  We  use 
the  molecular  weight  value  corresponding  to  the  monomer  HF  (20  kg/kmole).  Hence, 
the  concentration  results  expressed  in  ppm  are  monomer  equivalent  values. 
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DOWN  WIND  DISTANCE  (meters) 

FIGURE  4.1a:  Variation  of  Jet  Parameters  vs  Oist 
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FIGURE  4.1b:  Variation  of  Jet  Parameters  vs  Dist 
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DOWN  WIM)  DISTANCE  [meters) 

FIGURE  4.1c:  Variation  of  Jet  Parameters  vs  Dist 


evaporation  of  liquid  aerosols.  At  later  stages,  the  exothermic 
reaction  with  moistur  and  the  enthalpy  of  air  together  contribute 
to  the  gradual  rise  in  plume  temperature.  The  lowest  temperatures 
recorded  by  instruments  at  different  distances  from  the  source  are 
also  plotted  in  the  above  figures.  More  detailed  discussion  on  the 
temperature  predictions  are  given  in  a  later  section. 

The  jet  oenterline  velocity  decreases  -  a  requirement  of 
conservation  of  momentum  (with  air  mass  entrained)  and  ground 
friction.  We  terminate  the  jet  calculations  when  the  mean  velocity 
in  the  jet  is  within  5%  of  the  wind  speed  at  the  height  of  the  top 
of  the  plume.  Based  on  this  criterion,  the  "jet  length"  (i.e.,  the 
distance  from  the  orifice  within  which  the  jet  velocity  is  greater 
than  wind  speed  at  least  by  51)  has  been  calculated  for  all  three 
tests.  It  is  seen  that  these  jet  regimes  sxtsnd  up  to  176  m  for 
Tsst  #1,  140  m  for  Tost  #2,  and  108  m  for  Test  #3,  respectively. 
These  are  indicated  in  Table  4.3. 

It  is  notiosd  from  ths  results  prssented  in  ths  above  three  figures 
that  the  HF  concentration  at  the  centerline  varies  rapidly  reaching 
1  to  2%  level  by  the  end  of  the  jet  regime.  Hence,  in  the  jet 
regime,  a  2  order  of  magnitude  reduction  in  concentration  is 
effected.  In  the  remainder  of  the  dispersion  regimes,  another  2  to 
3  orders  of  magnitude  reduction  in  the  concentration  occurs. 

Ths  variation  of  the  width  of  the  plume  with  distance  is  also  shown 
in  the  above  three  figures.  It  is  seen  from  the  ADAM  results  that 
the  width  seems  to  increase  linearly  with  distance  -  a 
characteristic  of  axi-aymmetria  turbulent  jets.  This  implies  that 
the  lateral  expansion  due  to  negative  buoyancy  effects  is  small,  at 
least  in  the  near  field  where  jet  velocity  is  considerably  larger 
than  the  cross  stream  velocity  induced  by  negative  buoyancy. 

Only  for  one  HF  dispersion  test  (Test  #2)  was  an  aerial  photograph 
of  the  plume  available.  A  copy  of  this  photograph  is  indicated  as 
Plate  4.1.  This  picture  was  taken  from  a  aamera  at  500  m  above 
ground  and  500  m  upwind  of  the  spill  point.  Hence,  the  picture 
suffers  from  the  perspective  distortion.  We  have  determined  the 
apparent  horizontal  (cross  wind)  and  vertical  (downwind)  scale  by 
noting  the  distance  to  the  300  m  arc  of  masts  and  the  separation 
distance  between  two  instrument  polea  at  300  m  (this  separation 
distance  is  15  m)  .  Using  this  photographic  uata,  we  have  measured 
the  visible  plume  width  at  various  downwind  distances.  These  are 
plotted  as  data  points  in  Figure  4.1b. 

It  is  seen  that  the  predicted  width  is  larger  than  the  measured 
width.  A  number  of  reasons  can  be  put  forth  for  this  discrepancy. 
The  first  one  is  that  the  jet  model  over-dilutes  the  jet  (perhaps 
due  to  the  use  of  higher  than  actual  entrainment  coefficient).  The 
second  reason  is  that  our  determination  of  the  scale  of  the 
photograph  is  not  entirely  correct  since  we  have  not  considered  the 
effect  of  perspective  view  on  the  change  in  the  horizontal  scale  as 
one  goes  away  from  the  camera.  The  third  possible  reason  is  that 
in  the  above  comparison,  we  are  comparing  the  width  of  an 
equivalent  "box"  typo  of  jet  with  a  real  plume  whose  visible  width 


4-10 


TAIL!  4 .3 1  IILICTID  MUITI  MON  JIT  RCOINI  CALCULATION*  (1) 


MOTH  i 

1)  Tht  v*lut  of  orodynomle  rougHnoai  Zo  ■  0.0002  m. 

2)  Nlttlvt  hunldlty  <•  itiuntd  to  bo  17. M. 


4-11 


PI. ATE  4.1s  Experimental  Set  up  ami  Vapor  I’luma 
in**  Holi'anc  Trsit  at  I )( iE  Tint  Facility 


depends  on  both  the  aerosoi  content  of  th*'  plume  and  the  relative 
humidity  in  the  atmosphere.  However,  a  strong  argument  can  be  made 
for  modifying  the  entrainment  coefficient  for  two  phase  jets  with 
its  values  substantially  different  from  the  ones  for  a  single  phase 
jet  with  the  same  initial  momentum. 

We  are  unable  to  compare  other  predicted  results  (HF  concentration, 
jet  velocity,  aerosol  content,  etc.)  with  test  data  because  none  of 
these  were  measured  within  the  jet  region  in  the  field  experiments. 


4.3.2  Plume  Temperature 

The  plume  temperatures  predicted  by  the  thermodynamic  model 
(discussed  in  Chapter  2)  are  indicated  in  Figure  4.2a  and  Figure 
4.2b  for  the  conditions  of  uhe  three  Goldfish  dispersion  tests.  In 
Figure  4.2a,  the  relationship  batwaen  cloud  tsmperature  and  total 
HF  concentration  (monomer  equivalent)  over  the  entire  range  of  0  to 
95%  is  presented,  m  figure  4.?b,  the  temperature  variation  over 
the  HF  concentration  range  of  0%  to  5%  is  shown.  It  is  assn  that 
the  plume  temperature  is  above  the  (initial)  saturation  temperature 
of  290  K  for  HF  concentration  values  bnlow  3  to  4.5%,  depending  on 
the  test  conditions.  Using  thooe  results  and  tho  ADAM  calculated 
results  for  downwind  plume  centerline  concentration  at  ground  level 
with  distance,  we  obtain  the  variation  of  temperature  at  ground 
level  with  distance.  These  calculated  results  are  plotted  in 
Figure  4.3a,  b  and  c,  and  respectively  for  Tests  #1,  #2,  and  # 3. 

It  is  indicated  by  Blewitt,  et  al  (1987)  that  plume  temperatures 
were  measured  at  20m,  60m,  100m,  ?.00m  and  1000m  from  the  spill 
point.  The  plume  temperature  data  are  provided  (ibid)  for  the  20m, 
60m,  100m,  and  200m  ae  function  of  time  for  all  three  tests.  We 
have  assumed  that  these  data  refer  to  tho  plume  centerline  at  1 
meter  height  above  ground  (Blewitt,  et  al,  do  not  specify  this 
important  information  in  their  paper) .  The  loweat  measured 
temperature  at  each  location  is  plotted  as  a  function  of  the 
measuring  location.  Theso  are  indicuted  in  Figure  4.3a,  Figure  4.3b 
and  Figure  4.3c.  Table  4.4  aloe  shows  these  "measured"  lowest 
temperature  data  for  all  three  experiments.  Particularly 
noteworthy  are  the  significant  differences  in  the  reported  data  for 
test  #  3. 

Discussions  on  the  Temperature  Results 

In  general,  the  temperatures  measured  in  the  three  Goldfish  tests 
and  those  predicted  by  ADAM  are  in  reasonable  agreement  as  can  bo 
seen  from  Figuru  4.3a  through  Figure  4.3c.  It  should  be  noted  that 
the  temperatures  predicted  by  the  model  are  for  a  steady  and 
uniform  release;  therefore,  the  plume  temperature  at  specified 
downwind  centerline  location  does  not  vary  with  time.  The  data 
published  by  Blewitt,  et  al.,  (1987)  show  the  temperature  at  the 
location  of  the  instruments  to  vary  with  time,  significantly,  even 
though  the  flow  rata  out  of  the  pipe  is  constant  over  a  significant 
period  of  the  test.  We  speculate  that  this  temporal  variation  in 
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TABLE  4.4:  GOLDFISH  TEST  DATA  FOB  PLUME  TEMPERATURE 


DISTANCE  FROM 

SOURCE  (m) 

TEST  #1 

TEST  02 
(Note  2) 

TEST  n 
(Note  2) 

TEST  n 

(Nott  3) 

20 

2S3.0 

240.0 

264.0 

271.0 

40 

N/A 

N/A 

N/A 

291.0 

60 

255,0 

275.0 

283.0 

299.0 

100 

290.8 

299.0 

298.6 

301.0 

200 

3C1.0 

mgm 

302.0 

304.0 

1000 

310.0 

309.0 

307.0 

307.0 

NOTES: 

1)  Awbltnt  air  temperature  valuta  art  Indicated  In  Tibia  4.1.  These  valuta  art  uatd  In  the 
abovt  tibia. 

2)  Thtat  rtftr  to  tha  lowest  matured  temperature  as  read  from  tht  graphs  pubtlihed  by  Blewett, 
at  al  (1087).  Thtat  graphs  provldt  tht  valuta  for  tawperatura  dtprtaalon  with  raipact 

to  tht  ambient  temperature. 

3)  Obtained  from  tht  data  prtatnttd  In  Blawitt,  ft  al  (1990). 

NA:  Not  Avallablt 
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the  experiments  may  have  been  caused  by  one  or  more  of  the 
following:  i)  meandering  of  the  plume,  ii)  accumulation  of  HF 
liquid  droplets  on  the  thermocouple  and  their  subsequent 
evaporation,  or  iii)  an  unknown  instrument  problem. 

It  can  be  seen  from  Figure  4.2a  that  the  thermodynamic  model  based 
minimum  temperature  that  can  be  attained  in  the  cloud  for  the 
conditions  of  any  of  the  three  tests  is  260  K.  In  Test  #1,  the 
measured  minimum  temperature  is  253  K  (see  Figure  4.3a).  It  is 
doubtful  that  this  minimum  temperature  measured  is  accurate  because 
of  the  very  noisy  temperature  data  output  from  the  recorders.  The 
ADAM  predicted  temperature  for  Test  #1  conditions  slightly 
overpredicts  the  plume  temperature  for  distances  less  than  100  m 
and  underpredicts  at  distances  greater  than  100m  when  compared  to 
the  lowest  temperatures  measured  at  the  different  points.  In  the 
case  of  Test  #2  and  Test  #3,  the  predicted  temperature  is  very 
close  to  the  measured  values.  In  Figure  4.3c  are  shown  the 
temperature  data  indicated  by  Blewitt,  et  al.,  in  their  two 
publications.  The  discrepancy  in  the  measured  temperature 
(reported)  can  be  clearly  seen.  In  view  of  the  uncertainty  in  the 
measured  values  of  the  plume  temperature  and  noting  that  in  Figures 
4.3a  through  4.3c,  we  are  comparing  the  minimum  temperature 
measured  at  any  point  to  the  ADAM  predicted  temperature  (which 
assumes  a  steady  flow  both  in  time  and  in  direction)  the  agreement 
between  the  measured  and  predicted  temperatures  is  very  good. 

The  predicted  temperature  curves  shown  in  Figure  4.3a  through 
Figure  4.3c  show  some  erratic  behavior  at  about  300  m.  This 
distance  is  close  to  the  point  at  which  all  of  the  liquid  in  the 
aerosol  form  evaporates  (see  Table  4.3).  The  kinks  in  the  curve 
may  be  a  result  of  some  difficulty  in  convergence  in  the 
thermodynamic  program.  The  possible  model  temperature  variation  at 
distances  beyond  about  300  -  400  m  is  shown  by  dotted  lines  in 
Figures  4.3a  to  Figure  4.3c. 

4.3.3  Heavy  Gas  Dispersion  Regime 

Concentration  Data  and  Predictions 

We  have  obtained  the  measured  HF  concentration  data  principally 
from  two  references,  namely,  Blewitt,  et  al.  (1987),  and  Blewitt, 
et  al.  (1990) . 

The  concentration  data  in  the  Goldfish  tests  were  measured  at  300m, 
1000m  and  3000m  distances  (Blewitt,  et  al.,  1987).  Based  on  the 
information  in  the»  paper  by  Blewitt,  et  al.  (1987),  we  infer  that 
the  concentrations  reported  represent  average  values  over  66.6 
seconds,  83.3  seconds  or  100  seconds  (no  additional  details  are 
available  as  to  which  sensors  were  set  to  what  averaging  times) . 
Blewitt,  et  al.,  have  provided  concentration  contours  for  300,  1000 
and  3000m  locations  (in  general)  at  two  different  times.  Also 
provided  in  the  referenced  paper  is  a  3-D  contour  map  of 
concentration  as  a  function  of  crosswind  position  and  time  at  a 
specified  sensor  elevation. 
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Table  4.5  shows  the  HF  concentration  data  obtained  from  the  above 
references.  The  numbers  indicated  in  the  table  were  measured  off 
the  figures  published  in  the  two  papers.  To  the  extent  that 
accurate  measurements  cannot  be  made  from  small  figures,  we 
anticipate  some  errors  (at  best,  a  15%  error  in  reading  off 
graphs) .  It  can  be  clearly  seen  that  the  concentration  values  at 
specified  locations  vary  quite  significantly,  not  only  between 
different  figures  in  the  same  publication  but  also  between  the  two 
referenced  publications. 

The  ADAM  predicted  values  for  the  variation  of  plume  centerline, 
ground  level  concentrations  with  downwind  distance  are  plotted  in 
Figure  4.4a,  Figure  4.4b,  and  Figure  4.4c  for  the  conditions  of 
Test  #1,  Test  #2,  and  Test  #3,  respectively.  Also  plotted  on  the 
same  figures  are  the  respective  measured  values  (indicated  in  Table 
4.5).  The  predicted  concentrations  indicated  are  the  "monomer 
equivalent"  values.  I.e.,  the  cumulative  concentration  of  all 
oligomers  of  HF  (both  vapor  and  liquid)  are  expressed  as  if  they 
had  a  molecular  weight  of  20. 


Discussions  on  the  Concentration  Results 

There  are  no  concentration  data  for  distances  less  than  300  m.  The 
jet  plume  region  is  within  this  distance.  Hence,  we  are  unable  to 
compare  the  "near  field"  concentration  predictions  with  test  data. 

The  predicted  concentration  in  Figure  4.4a  is  within  a  factor  of  2 
with  the  data  presented  by  Blewitt,  et  al.,  (1990).  However,  in 
view  of  the  significantly  large  data  scatter,  it  can  be  concluded 
that  the  predicted  curve  is  well  within  the  acceptable  accuracy. 

In  Figure  4.4b,  the  predicted  concentration  values  are 
significantly  lower  than  the  experimental  values.  In  fact,  in  this 
Test  #2,  the  data  scatter  do  not  seem  to  be  large.  Unfortunately 
for  this  test,  no  data  have  been  reported  for  the  3000  m 
instruments. 

The  predicted  values  for  Test  #3,  shown  in  Figure  4.4c,  seem  to 
agree  very  closely  with  the  data  presented  in  the  first  paper  by 
Blewitt,  et  al.  (1987).  Compared  to  the  data  indicated  in  the 
second  paper  (1990),  the  predicted  values  are  lower  in  some 
locations  by  as  much  as  a  factor  of  3.  We  are  unable  to  explain 
such  a  large  discrepancy. 

The  concentration  predictions  by  ADAM  are  based  on  a  concentration 
averaging  time  of  66  seconds.  This  value  was  chosen,  based  on  the 
information  indicated  by  Blewitt,  et  al.  (1987)  on  the  design  of 
the  HF  samplers  in  the  tests;  the  sampling  time  for  the  filter 
cassettes  was  66.6s,  or  83.3s  or  100s.  We  also  note  that  the  ADAM 
code  predictions  are  based  on  the  assumption  of  a  steady  state 
plume  (i.e.,  the  source  flow  was  steady).  The  data  provided 
Blewitt,  et  al.,  are  difficult  to  interpret.  As  an  example,  we 
consider  the  conditions  of  Test  #3.  The  duration  of  the  steady 
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TABLE  4. Si  MEASURED  DATA  ON  HE  CONCENTRATIONS 


HE  CONCENTRATION 

REFERENCE  0 

(ppn) 

AND  SUBREFERENCE  i 

3100 

Figure  4.1  of  Raf(A) 

21000 

Elgura  4.2  of  Raf(A) 

20150 

Figure  12  of  Ref(B) 

17920 

Elgura  4.3  of  Ref(A) 

1500 

Figure  4.6  of  Rof(A) 

2700 

Figure  4.7  of  Ref(A) 

2300 

Figure  4.6  of  Raf(A) 

3000 

Figure  12  of  Raf(B) 

400 

Figur#  4.10  of  R«f(A) 

ME  CONCENTRATION 

REFERENCE  » 

<PE») 

AND  SUBREFERENCE 

16300 

Figure  13  of  Ref(B) 

9300 

Figure  4.11  of  Ref (A) 

6300 

Figure  4.12  of  Rof (A) 

1900 

Figure  4.13  of  Rof(A) 

1950 

Figure  13  of  Ref(B) 

ME  CONCENTRATION 

REFERENCE  0 

(ppm) 

AND  SUBREFERENCE 

6300 

Figure  4.15  of  Ref(A) 

7600 

Figure  4.16  of  Ref(A) 

16420 

Figure  4.17  of  Ref(A) 

15700 

Figure  14  of  Ref(B) 

600 

Figure  4.18  of  Ref (A) 

700 

Figure  4.19  of  Ref(A> 

2323 

Figure  4.20  of  Ref(A) 

2340 

Figure  14  of  Ref(A) 

222 

figure  4.22  of  Ref(A) 

200 

Figure  14  of  Ref(B) 

NOTE It 

Ref(A)  ll  the  paper  by  Bleultt,  ft  ft.,  198T. 

•  Atf(B)  ll  the  paper  by  Bleultt,  ft  •(.,  1990. 

*  All  reneentratlon  valuai  refer  to  location*  1  m  abovi  the  ground  at  ptune  centerllr 
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figure  4.4b:  Ground  Level  Concentration  vs  Distance 
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figure  4.4c:  Ground  Level  Concentration  vs  Distance 


state  release  was  about  380  seconds.  Blewitt,  et  al.  (1987)  give 
two  different  concentration  contours  for  113s;  in  one  the  peak  1  m 
level  concentration  is  6300  ppm  and  for  the  other  it  is  7800  ppm. 
These  data  were  obtained  from  the  contour  plots  (Figure  4.15  and 
Figure  4.16)  indicated  by  Blewitt,  et  al.  (1987).  What  is  not 
clear  from  these  figures  is  the  averaging  time  used  for  the  plot  of 
the  concentration  contours.  other  questions  are:  Do  the  data 
presented  refer  to  only  one  filter  cassette?  What  was  the 
variation  of  concentrations  measured  at  the  same  location  by 
different  filter  cassettes?  What  is  th  basis  of  the  3-D 
concentration  contour  figures  presented  by  Blewitt,  et  al.  (1987)? 
Finally,  why  are  there  discrepancies  between  concentration  values 
published  by  the  same  author  in  two  different  publications  for  the 
identical  test  and  measurement  locations? 

Also  not  certain  is  the  algorithm  by  which  the  experimentally 
measured  concentration  values  were  converted  to  molar 
concentrations  which  are  reported  in  the  papers  by  Blewitt,  et  al. 
In  the  tests,  integrated  mass  dose  of  total  HF  (vapor  and  liquid 
aerosols,  if  any)  are  collected  for  66s,  88.3s  or  100s  on  filter 
paper.  The  aspiration  rate  is  given  to  be  3.5  liters  per  minute. 
Since  no  real  time  measurement  of  the  temperature  or  molecular 
weight  of  the  species  were  made  at  the  positions  where 
concentration  measurements  were  made,  it  is  our  opinion  that  it  is 
not  possible  to  convert  the  dose  data  (actually  obtained  in  the 
tests)  to  the  ppm  units  reported.  It  is  possible  that  a  monomer 
equivalent  is  assumed  (i.e.,  a  molecular  weight  of  20  kg/kmole) . 
Even  in  this  case,  the  temperature  of  the  cloud  must  be  assumed  in 
order  to  convert  the  raw  data  to  the  ppm  values  for  the 
concentration.  Blewitt,  et  al.,  have  not  indicated  how  this 
conversion  to  the  reported  units  are  made. 

In  view  of  the  many  uncertainties  in  the  reported  data,  both  in 
magnitude  and  the  way  in  which  the  raw  data  may  have  been  converted 
into  the  (reported)  molar  concentration  values,  we  conclude  that 
the  ADAM  predictions  represent  the  test  data  reasonably  well. 
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CHAPTER  8 


FLUORINE 

MODELS  FOR  AIR  MIXING  AMD  DISPERSION 


5.1  INTRODUCTION 

Under  normal  ambient  condition*,  fluorine  is  a  light  yellow 
"permanent"  gas  with  a  critical  temperature  of  144  K.  Due  to  its 
very  reactive  character,  the  chemical  is  shipped  commercially  in 
ganeous  state  in  small  steel  cylinders  (maximum  2.7  kg)  at 
relatively  low  pressures.  However,  when  larger  quantities  are  to 
be  shipped,  the  chemical  is  transported  in  a  liquid  state  near  its 
normal  boiling  point  of  85  K.  Accidents  involving  there  tanks 
could  lead  to  the  release  of  large  quantities  of  fluorine  and  the 
consequences  of  such  an  event  provide  the  focus  of  the  present 
discussion. 

5.1.1  Details  of  LF2  Transport 

Liquid  fluorine  (LF3)  is  transported  in  road  tanker  trailers  which 
consist  of  three  horizontal,  aoncentric  vessels.  LF3  is  contained 
within  the  inner  tank  made  of  monel  or  stainless  steel  304.  An 
outer  concentric  tank  made  of  stainless  steel  304  contains  liquid 
nitrogen  and  this  tank  is  enclosed  in  a  carbon  steel  cylinder.  The 
annular  space  between  the  nitrogen  tank  and  the  carbon  steel  tank 
is  filled  with  powdered  insulation  and  evacuated.  Figure  5.1  shows 
a  cross-sectional  view  of  the  LF3  and  other  tanks  on  the  road 
tanker.  Figure  5.2  shows  the  tank  flow  system  in  the  fluorine 
trailer. 

The  inner  product  vessel  has  a  capacity  of  1.77  m1  (468  gallons) 
with  an  ullage  volume  of  about  1.5  mJ.  That  in,  the  maximum  product 
mass  is  2260  kg  (5000  lbs) .  This  vessel  has  a  service  rating  of 
0.482  MPag  (70  psig).  Manual  operation  of  the  fluorine  fill  valves 
in  the  vapor  phase  is  the  only  form  of  pressure  relief  device 
provided.  In  case  of  an  emergency,  a  fifty-foot  coil  of  copper 
tubing  with  the  necessary  fittings  is  available  for  the  purpose  of 
venting  excess  fluorine  pressure.  Changes  in  the  temperature  of 
the  product  are  monitored  with  the  use  of  fluorine  product  pressure 
gauge.  Normal  indication  on  this  pressure  gauge  is  0.  The  high 
pressure  alarm  is  sat  at  0.124  MPag  (18  psig)  which  activates  the 
warning  lights  and  horn. 

The  temperature  of  the  liquid  fluorine  is  maintained  by  the  middle 
vessel  of  liquid  nitrogen  (LN2)  whose  boiling  point  is  77.4  K 
(-  360  8F) .  This  temperature  is  about  8  K  (14.5  °F)  less  than  the 
boiling  point  of  LF5.  The  annular  space  between  the  nitrogen  tank 
arid  the  outer  carbon  steel  tank  is  filled  with  powdered  insulation 
and  evacuated  to  less  than  one  millimeter  of  mercury,  absolute 
pressure.  The  small  heat  leak  from  the  atmosphere  causes  the  LN} 
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to  boil  and  the  vaporized  nitrogen  is  vented  to  the  atmosphere  thus 
maintaining  the  LN2  at  its  boiling  temperature. 

The  liquid  nitrogen  tank  has  a  volume  capacity  of  1.287  m5  (340 
gallons).  The  boil-off  rate  is  about  30.3  liters/day 

(8  gallons/day) .  The  vessel  has  a  service  pressure  rating  of  0.31 
MPag  (45  psig)  and  is  equipped  with  n  spring  loaded  pressure  relief 
valve  set  at  0.207  MPAg  (30  psig)  and  a  pressure  relief  rupture 
disk  designed  for  0.276  MPag  (40  psig)  rupture. 

Table  5.1  indicates  the  principal  specif icationa  and  transport 
volumes  of  the  fluorine  road  tanker. 


5.2  Properties  of  Fluorine 

5.2.1  Physical  Propsrtiss 

Fluorins  is  normally  a  gas  at  ambient  temperature  and  pressure. 
Its  aritical  tsmpsrature  is  144  X.  That  is,  in  order  to  liquefy 
the  element,  it  haa  to  be  cooled  to  temperatures  below  144  K. 
Fluorine  has  a  sharp,  penetrating  odor  detectable  by  human  beings 
at  as  low  a  concentration  as  0.1  ppm.  The  principal  thormodyriamir 
properties  of  Fluorine  of  interest  to  this  study  are  indicated  in 
Table  5.2. 

5.2.2  Chemical  Properties 

Fluorine  is  a  very  reactive  element  and  has  the  ability  to  combine 
with  most  organic  and  inorganic  materials  at  or  below  room 
temperature.  Reactions  with  organia  materials  are  strongly 
exothermic.  Organic  and  hydrogen  containing  compounds  especially 
can  burn  or  explode  when  exposed  to  fluorine.  Elemental  fluorine 
reacts  directly  with  the  noble  gases  xenon,  radon,  and  krypton  to 
form  fluorides.  Nitrogen  and  oxygen  form  fluorides  but  do  not 
react  directly  with  fluorine  except  in  the  presence  of  an  electric 
charge. 

Fluorine  is  the  first  member  of  the  halogen  family  and  has  the 
lowest  enthalpy  of  dissociation  relative  to  the  other  halogens, 
which  is  in  part  rssponsible  for  its  greater  reactivity. 
Furthermore,  the  strength  of  the  fluorine  bond  with  other  atoms  is 
greater  than  that  of  the  other  halogens.  The  reactivity  of 
fluorine  compounds  varies  from  extremely  stable  compounds  (such  as 
SF6 ,  NF3,  and  fluorocarbons)  to  extremely  reactive  compounds  such 
as  the  halogen  fluorides.  Volatile  metal  compounds  such  as  WF6 , 
MoF6  are  produced  by  the  reaction  of  the  metal  with  elemental 
fluorine. 

Most  reaction  chemistry  studies  involving  fluorine  have  been 
conducted  at  ambient  temperatures  and  above,  and  the  discussion 
below  relates  to  findings  in  this  temperature  range.  Leaks  of  vary 
cold  fluorine  gas  or  liquid  may  ahow  a  different  behavior,  i.e., 


TABLE  5.1 


LIQUID  FLUORINE  TRANSPORT  SPECIFICATIONS 


1.  Tank  Specifications;  Liquid  Fluorine 

Length 

Diaaeter 


Voluae 

1.770 

BA3 

Tank  Pressure  Rating 

590.000 

kPa  abs 

2.  Specifications  of  Liquid  Fluorine  (LF2) 

Voluae  transported 

1.500 

ba3 

Hass  of  LF2  in  tank 

2260.000 

kg 

Liquid  Teaperature 

77.400 

K 

Tank  Pressure  ** 

108,000 

kPa  (abs) 

3.  Tank  Specifications:  Liquid  Nitrogen 

Length 

Voluae 

1.287 

n'3 

Tank  Pressure  Rating 

310.000 

kPa  g 

4.  Specifications  of  Liquid  Nitrogen  (LN2) 

Voluae  transported 

1.287 

a*3 

Hass  of  LN2  in  tank 

1018.000 

kg 

Liquid  Teaperature 

77.400 

K 

Tank  Pressure  ** 

100.000 

kPa  (abs) 

**  The  ullage  voluae  of  the  LF2  tank  is  filled  with  gaseous 
Heliuo  which  aaintains  an  above  anbient  pressure  in  tank. 

LF2  partial  pressure  at  the  liquid  nitrogen  boiling 
teaperature  (77.4  K)  is  41.4  kPa  (abs)  or  0.41  ata  abs. 
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TABLE  5.2 


THERMODYNAMIC  PROPERTIES  OF 
FLUORINE  and  NITROGEN 


Property 

Fluorine 

Nitroqen 

Units 

Molecular  Weight 

38.00 

28.00 

kq/kaol 

Critical  Teaperature 

144.30 

126.20 

K 

critical  Pressure 

5.52  E+06 

3.39  E+06 

N/l‘2 

Boiling  Point 

84.45 

77.40 

K 

Freer inq  Point 

53.50 

63.30 

K 

Liquid  Density 
(  §  NBT  ) 

1559.30 

790.90 

kg/1'3 

Sat  Liq  Enthalpy 
(  §  NBT  ) 

0.00 
(«  NBT) 

0.00 

(S  273.16  K) 

J/ltq  K 

Sat  Vap  Enthalpy 
(  e  NBT  ) 

1.73  E+05 

NA 

J/kq  K 

Equations  for  Teoperature  dependent  paraaeters  (all  in  SI  units) 

FLUORINE: 

Saturated  Vap  Pr  p  ■  per  *  exp{ (-6.1S224  y  +  1.18  y'1.5  -  1.16555  y'3  •  1.50167  y*6]/(l-y|) 

with  y  «  (1-T/Tc) 

Liq  Heat  Cap  CL  *  1087  +  1.79  T 

Heat  of  Vaporization  XL  »  1.73  E+05  ((l-y)/(l-yr))'0.38  with  yr  ■  Tref/Tc 

Tref  *  84.45  K 

NITROGEN: 

Saturated  Vap  Pr  p  =  per  *  exp([-6. 09676  y  +  1.1367  y'1.5  -  1.04072  y*3  -  1.93306  y‘6]/(l-y)} 

with  y  =  (1-T/Tc) 

Liq  Heat  Cap  CL  -  2309  -  11.1  T  +  0.1  TA2 

Heat  of  Vaporization  XL  >  1.99  E+05  ((l-y)/(l-yr)('0.38  with  yr  «  Tref/Tc 

Tref  *  77.4  K 
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slower  reactions,  but  it  is  possible  that  there  could  be  local 
wanning  of  pockets  of  fluorine  so  that  it  is  only  prudent  to 
consider  reactions  occurring  under  ambient  conditions. 


a)  Reactions  with  Water:  There  seems  to  be  some  controversy 
regarding  the  reactivity  of  fluorine  with  water  or  water  vapor. 
Cady  and  Burger  (1950)  state  that  fluorine  reacts  with  water  to 
form  oxygen  fluoride  and  hydrogen  fluoride. 

2P2  +  H20  -  0F2  +  2 HP  (5.1) 

Rudge  (1962)  considers  OF,  to  be  unstable  and  to  decompose  to  02  and 
F2.  He  also  claims  that  h202,  fluorine  monoxide  (F,0) ,  and,  in  some 
cases,  ozone  are  formed.  Landau  and  Rosen  (1951)  state  that 
"fluorine  normally  reacts  vigorously  with  water  in  the  form  of 
either  a  vapor  or  a  liquid  to  form  hydrofluoric  acid  and  oxygen. 
Nevertheless,  inhibition  has  frequently  been  observed;  that  is, 
quantities  of  water  vapor  may  accumulate  in  the  presence  of 
fluorine  until  the  reaction  is  suddenly  initiated  with  explosive 
violence."  They  also  note  that  "When  fluorine  emerges  into  the  wet 
atmosphere,  e.g.,  through  a  leaking  pipe  flange,  it  reacts  with 
condensed  water  and  forms  sparks  and  flashes  that  are  readily 
visible  at  a  considerable  distance."  The  reaction 

F2(g)  +  H20(g)  -  (1/2)0.  (g)  +  2HF (g)  (5.2) 

is  very  exothermic  with  an  enthalpy  of  reaction  of  about  8  MJ/kg  P2 
(Stull  and  Prophet,  1971)  and  one  kg  of  fluorine  reacts  with  about 
0.47  kg  of  water. 

Thus,  while  there  is  some  uncertainty  concerning  the  rates  of 
reaction  between  fluorine  and  water,  most  authors  agree  that 
reaction  is  possible.  At  low  temperatures,  rates  of  reaction 
should  be  small.  Also,  as  air  dilutes  the  fluorine,  the  reactivity 
decreases  significantly  as  discussed  below. 


b)  Reactions  with  Organic  Materials:  All  investigators  warn  that 
fluorine  in  high  concentrations  will  "burn"  organic  materials  with 
the  presence  of  flames.  The  enthalpies  of  reaction  are  high  and 
much  larger  than  comparable  enthalpies  of  oxidation.  For  example, 
methane  reacts  (in  theory)  to  form  CF4  and  HF.  The  enthalpy  of 
this  reaction  is  »  -44  MJ/kg  CH4  whereas  if  CH4  reacts  with  02  to 
form  C02  and  water,  AH  »  -2lMJ/kg  CH4. 

The  probability  of  ignition  between  fluorine  and  organic  compounds 
appears  to  be  a  strong  (but  as  yet  unquantified)  function  of  the 
fluorine  concentration.  In  some  crude  but  illuminating  tests, 
Landau  and  Rosen  (1951)  placed  organic  materials  various  distances 
from  a  jet  of  fluorine  issuing  from  an  orifice  on  a  fluorine  gas 
tank  at  4.5  bar.  The  orifice  sizes  were  3  and  9  mm  in  diameter. 
Cotton,  wool,  and  wood  surfaces  were  used.  Ignition  occurred  when 
the  distances  were  less  than  about  15  cm.  In  a  very  few  tests 
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where  the  fluorine  was  diluted  to  20  volume  percent  with  nitrogen , 
no  ignition  took  place  even  when  the  material  (cotton  cloth  in  most 
cases)  was  in  very  close  proximity  to  the  jet.  Also,  with  20% 
fluorine/80%  nitrogen,  there  were  no  violent  reactions  with  butyl 
rubber  or  neoprene  gaskets,  but  the  physical  strength  of  the 
rubbers  deteriorated  rapidly. 

One  concludes  that  fluorine  reacts  violently  with  organic  materials 
when  pure  and  at  ambient  temperatures  and  above.  The  reactivity  of 
fluorine  decreases  to  lower  levels  when  diluted  to  10-20%  or 
below.  The  effect  of  temperature  is  unknown,  but  it  is  reasonable 
to  assume  the  reactivity  diminishes  as  the  temperature  is  lowered. 
Thus,  in  an  accident,  if  air  dilution  can  occur  before  the  fluorine 
warms,  there  is  a  lower  probability  of  ignition  with  organic 
materials  (and  water  vapor)  in  the  immediate  vicinity  of  the  spill. 


c)  Reactions  with  Metals  and  Other  Materials:  Due  to  the 
formation  of  a  passive  metal  fluoride  surface  film,  the  reaction  of 
fluorine  with  most  common  metals  (iron,  stainless  steel,  Monel, 
copper,  aluminum)  is  relatively  slow  at  room  temperature  and  below 
but  is  vigorous  at  elevated  temperatures.  Landau  and  Rosen  (1951) 
state  that  leaks  of  ambient  temperature  fluorine  gas  can  burn  iron 
piping  should  there  be  organic  materials  or  water  vapor  in  the 
vicinity  of  the  leak  to  initiate  the  reaction.  Experiments  with 
clean  copper,  brass,  and  stainless  steel  tubing  showed  no  reaction 
with  pure  fluorine  gas  flawing  within  the  tubing.  If,  however,  the 
tubing  were  contaminated  by  oil  or  grease,  tha  copper  tubing  burned 
while  the  brass  and  stainless  steel  tubing  were  heated  to  a  rad 
heat.  No  reactions  were  observed  in  any  of  the  tubing,  even  when 
oily,  when  the  gas  was  20%  fluorine  in  nitrogen  (Landau  &  Rosen, 
1951) . 

Many  inorganic  compounds  react  with  fluorine  with  the  fluorine 
replacing  oxygen  in  tha  material,  but  usually  only  at  temperatures 
well  above  ambient.  Rudge  (1962)  states  that  at  room  temperature, 
in  the  absence  of  hydrofluoric  acid,  "fluorine  is  without 
detectable  effect  on  glass." 


5.3  Release  Scenarios 

There  are  basically  two  fluorine  release  scenarios  that  nead  to  be 
considered  in  modeling  the  dispersion  hazard.  These  scenarios 
arise  because  of  the  nature  of  the  chemical  and  the  transport 
system  used.  These  two  scenarios  are  discussed  below. 


5.3.1  Explosive  Release  of  Fluorine  and  Flash  Vaporization 
(SCENARIO  A:) 

consider  an  accident  or  an  operational  malfunction  in  which  the 
nitrogen  coolant  is  drained  out  of  the  outer  tank  through  a  rupture 
hole  and  released  into  the  atmosphere.  It  is  assumed  that  the 
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inner  LF2  tank  is  not  damaged.  In  such  an  accident,  the  loss  of 
the  LN2  cooling  will  result  in  the  temperature  of  liquid  fluorine 
to  increase  with  a  corresponding  increase  in  the  fluorine  tank 
pressure.  The  fluorine  tank  has  a  service  rating  of  0.59  MPa 
(abs)  .  At  this  pressure,  the  temperature  of  fluorine  is  only  about 
20  K  higher  than  its  normal  boiling  point  of  85  K. 

Continued  heat  leak  beyond  the  point  at  which  the  tank  pressure  is 
0.59  MPa  may  lead  to  the  rupture  of  the  tank  and  the  consequent 
very  rapid  depressurization.  In  this  scenario,  the  compressed, 
supersaturated  fluorine  liquid  will  flash  instantly  into  vapor  and 
liquid.  The  liquid  may  be  entrained  in  the  vapor  cloud  in  the  form 
of  a  fine  mist  at  a  temperature  of  84.5  K.  Some  part  of  the  liquid 
may  spill  on  the  ground,  but  it  is  assumed  that  the  violence  of  the 
release  may  ensure  that  all  of  the  liquid  is  entrained  in  the  vapor 
cloud.  This  scenario  is  depicted  schematically  in  Figure  5.3. 

The  violence  of  release  will  also  entrain  the  ambient  air  into  the 
vapor  cloud.  The  density  of  this  vapor  cloud  containing  liquid  F2 
aerosols  will  be  higher  than  that  of  air;  consequently  the  vapor 
cloud  will  disperse  downwind  as  a  heavy  vapor  cloud.  During  this 
dispersion,  additional  air  will  be  entrained,  liquid  fluorine 
aerosols  will  evaporate,  water  vapor  from  the  atmosphere  will 
condense  into  the  cloud  and  "react"  with  F2  forming  HF.  All  of 
these  thermodynamic  phenomena  are  discussed  in  section  5.4. 


5.3.2  Liquid  Release  and  Boiling  Pool  Source  of  Vapor  (SCENARIO 
B!) 

In  an  accident,  instead  of  just  the  nitrogen  tank  being  punctured, 
both  nitrogen  and  fluorine  tanks  are  punctured,  then  two  cryogenic 
liquid  streams  will  issue  out  of  the  damaged  section  of  tho  tank. 
Because  both  liquids  will  be  at  about  the  saturation  temperature  of 
nitrogen  (i.e,,  77  K) ,  there  will  be  no  flash  vaporization  during 
release. 

When  the  two  liquids  form  a  puddle  on  the  ground,  rapid 
vaporization  is  expected  to  occur  because  of  heat  leak  from  the 
relatively  warm  ground  and  from  the  ambient  wind.  Only  vapors  of 
nitrogen  and  fluorine  will  be  generated.  The  vapor  evolved  will  be 
carried  downwind  and  dispersed.  This  scenario  of  release  is  Bhown 
in  Figure  5.4. 

In  section  5.5,  a  source  model  is  developed  to  analyze  the 
distillation  and  evaporation  from  a  pool  containing  two  volatile 
liquids.  More  details  of  this  model  are  indicated  in  Appendix  B. 


5.3.3  Flashing  of  the  Liquid  Released 

In  Scenario  A,  the  liquid  fluorine  is  released  from  an  elevated 
pressure  and  temperature  (compared  to  the  normal  boiling 
temperature) .  This  results  in  the  flashing  of  the  liquid  forming 
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TIPPED  OVER  TRAILER 


TIPPED  OVER  TRAILER 


BOTH  INNER  AND  OUTER  TANKS 
PUNCTURED  BY  ACCIDENT 


GROUND 


VAPOR  CLOUD 
CONTAINING  ONLY  VAPOR 
<L  OF  COLD  FLUORINE  AND 
NITROGEN 


T?  /  77 


V— ~  ...  BOILING  POC 


BOILING  POOL  OF  LN. 
AND  LF„  ' 


FIGURE  5.4:  Schematic  Representation  of  Scenario  "B"  Type  Release  of  Fluorine 
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ambient  pressure  saturated  vapor  and  saturated  liquid.  This 
phenomenon  is  discussed  below  and  the  fraction  of  liquid  which 
flashes  is  calculated. 

The  phenomenon  of  flashing  of  a  compressed  liquid  when  released 
into  ambient  pressure  has  been  discussed  in  an  earlier  companion 
report  (see  section  2.4  of  report  by  Raj  and  Morris,  1987).  We 
repeat  here  the  formula  for  the  mass  fraction  of  vapor  generated  by 
the  flashing  process. 

sat  sat 

C  h  (T  )  -  h  (T  )  ] 

Li  LB 

f  »  - - - —  (5.3) 

v  X(T  ) 

B 

where,  fy  represents  the  mass  fraction  of  ambient  pressure 
saturated  vapor  produced,  h's  represent  the  saturated  liquid 
enthalpies  at  the  respective  saturation  temperatures  indicated  and 
X  represents  the  heat  of  liquid  vaporization  at  ambient  pressure  or 
at  the  normal  boiling  temperature. 

Table  5.3  shows  the  conditions  of  the  LF.  at  various  temperatures 
and  the  values  of  the  enthalpies  at  these  temperatures .  Also  shown 
are  the  mass  fractions  of  vapor  flash  for  LFZ  releases  from  various 
saturated  conditions.  These  conditions  correspond  to  the  states 
that  LF2  will  be  in  during  its  heating  from  the  loss  of  LNZ  coolant. 
The  density  of  the  final  mixture  of  vapor  and  liquid  is  also  shown 
in  this  table.  (In  generating  these  mixture  density  values,  it  is 
assumed  that  all  of  the  liquid  in  the  post-flash  stream  is 
entrained  into  the  vapor  phase.)  Figure  5.5  shows  the  mass 
fraction  of  vapor  generated  by  the  flashing  of  LF2  from  different 
temperatures  (or  pressures) .  Also  shown  in  the  figure  is  the 
variation  of  the  saturated  vapor  pressure  of  Fluorine  with 
temperature . 

It  is  seen  that,  if  LFZ  is  released  as  a  con  aquencs  of  the  tank 
failure  at  the  rated  pressure  (590  kPa  abs) ,  then  about  13.6  %  of 
the  mass  of  fluorine  released  will  flash  to  vapor. 

The  density  of  the  mixture  of  vapor  and  the  post  flash  liquid 
entrained  within  the  vapor  cloud  as  fine  aerosols  is  calculated  as 
follows: 


Pmi*  "  - J - -  (5,4) 

[  v  +  ] 

Pv  Pi 

where,  pffllx  is  the  density  of  the  mixture  of  vapor  and  entrained 
liquid  aerosol,  py  and  are,  respectively,  the  saturated  vapor 
density  and  saturated  liquid  density  at  ambient  pressure  and  f  and 
f,  are  respectively  the  mass  fractions  of  vapor  and  liquid  produced 
after  the  flashing. 
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TABLE  5.3 

FLUORINE  FLASH  CALCULATIONS 


Liquid  Sat  Pressure  Liquid 

Temp  Enthalpy  ! 

(note  1) 

(K)  (kPa  abs)(atn  abs)  (J/kg  K) 

< —  Post  Flash  — > 
conditions 

Latent  Mass  Density 

Heat  of  fraction  of 
boiling  of  Vapor  Final 
forned  Vapor-liq 

Mixture 

REMARKS 

(kJ/kg)  (8)  (kg/nA3) 

77.4 

41.4 

0.41 

-1018.1 

180.5 

Liq  N2  saturation  tenperature 

80.0 

58.7 

0.58 

-655.8 

177.8 

83.0 

85.4 

0.84 

•222.6 

174.6 

84.5 

101.9 

1.01 

0.0 

172.9 

0.00 

1559.3  <—  Normal  Boiling  Point 

85.0 

107.9 

1.06 

619.4 

172.4 

0.36 

771.6 

90.0 

184.3 

1.82 

6837.5 

166.7 

3.95 

127.8 

95.0 

296.4 

2.92 

13100.4 

160.7 

7.57 

69.4 

100.0 

453.3 

4.47 

19408.0 

154.3 

11.22 

47.5 

101.0 

490.9 

4.85 

20674.9 

153.0 

11.95 

44.7 

102.0 

530.8 

5.24 

21943.6 

151.6 

12.69 

42.2 

103.0 

573.0 

5.66 

23214.0 

150.3 

13.42 

39.9 

103.3 

584.0 

5.76 

23531.9 

149.9 

13.61 

39.4  <-•  Pressure  at  which  the  tank  is  rated 

104.0 

617.7 

6.10 

24486.3 

148.9 

14.16 

3V.8 

105.0 

664.8 

6.56 

25760.4 

147.4 

14.90 

36.1 

Notes! 

1.  Liquid  enthalpy  reference  tenperature  is  the  norsal  boiling  point. 

That  is  at  84.5  K  the  saturated  liquid  enthalpy  is  equal  to  0. 

2.  Liquid  enthalpy  is  calculated  by  integrating  the  tenperature  dependent 
liquid  specific  heat  (see  Tablo  5.2)  x  dT  froo  the  reference 
tenperature  to  the  liquid  tenperature . 

3.  Heat  of  vaporization  is  obtained  fron  the  foraula  indicated  in  Table  5.2 

4.  The  saturated  vapor  density  at  ambient  pressure  *  5.48  kg/n'3 
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Fluorine  Saturation  Temperature  (K) 

FRACTION  FLASH  vs  FLUORINE  TEMPERATURE 

FIGURE  5.5 


The  final  mixture  density  of  the  vapor-aerosol  mixture  is  shown 
plotted  in  Figure  5.6  as  a  function  of  the  initial  temperature  of 
fluorine  in  the  tank.  This  density  is  calculated  on  the  basis  that 
all  of  the  post-flash  liquid  is  entrained  in  the  vapor  cloud 
formed.  It  is  seen  that  the  combined  density  of  the  vapor  cloud 
when  fluorine  is  released  from  a  pressure  corresponding  to  the 
rated  failure  pressure  of  the  tank  (590  kPa  abs)  is  39.4  kg/m*3. 
This  density  is  about  33  times  the  density  of  ambient  air  at  293  K1 


5.4  Thermodynamic  Model  for  Fluorine  Liquid  Aerosol,  Vapor  and 
Humid  Air  Mixing 

5.4.1  Description  of  the  Mixing  Phenomenon 

The  final  atate  of  the  fluorine-air-water  vapor  cloud  depends  upon 
the  relative  quantities  of  these  components  as  well  as  upon  any 
energy  input  from  the  ground,  vegetation,  etc.  Also,  the 
thermodynamic  state  of  the  cloud  varies  with  time  (or  downwind 
distance)  as  air  is  continually  entrained. 

We  consider  the  oaae  where  a  specified  quantity  of  fluorine  ia 
mixed  with  a  definite  amount  of  air  with  a  given  relative  humidity 
-  and  with  an  energy  input  that  ia  also  specified.  From  an  overall 
net  enthalpy  balance  (since  the  mixing  step  ia  ieobaric) ,  one  can 
determine  the  final  state.  In  addition,  equilibrium  relationships 
are  employed  to  ascertain  if  condensed  phases  (liquid  fluorine, 
water  ice,  or  water  liquid)  are  present. 

Before  mixing  with  any  air,  the  liquid  (aerosol)  and  vapor  fluorine 
are  in  equilibrium  at  the  given  ambient  pressure.  If  thie  pressure 
were  about  one  bar,  the  equilibrium  fluorine  temperature  is  »  85  K. 

When  mixed  with  air,  initially,  the  partial  pressure  of  the 
fluorine  decreases  and,  if  liquid  fluorine  is  prsssnt,  evaporation 
oaoura.  Thus,  while  the  warmer  air  tends  to  increase  the 
temperature,  as  does  the  condensation  (and  freezing)  of  water,  the 
fluorine  evaporation  step  ik.  the  more  important  and  the  cloud 
temperature  decreases.  This  trend  continues  with  the  admixture  of 
additional  air  until  the  liquid  fluorine  has  completely 
evaporated.  The  cloud  temperature  is  then  below  the  usual  boiling 
temperature  of  fluorine  (85  K)  and  essentially  all  water  is  in  the 
form  of  ice. 

Following  this  temperature  drop  caused  by  fluorine  evaporation, 
mixing  with  more  air  results  in  a  rise  in  temperature.  As  this 
continues,  some  water  ioe  sublimes  so  as  to  maintain  its 
equilibrium  vapor  pressure.  When  the  mixture  temperature  reaches 

273.2  K,  the  ice  melts  to  liquid  water.  Above  this  temperature, 
only  liquid  and  vapor  water  are  present. 

The  scenario  given  above  is  true  evan  if  there  is  a  significant 
energy  input. 
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Fluorine  Saturation  Temperature  [K] 

POST  FLASH  VAP-AEROSOL  MIX  DENSITY  vs  TANK  TEMP 


5.4.2 


calculation  Procedure 


Thus  in  any  calculational  scheme,  one  must  allow  for  various  final 
"domains"  depending  upon  the  circumstances.  We  define  the  inlet 
conditions  by  specifying: 

o  mass  released,  fraction  flash,  and  temperature  of  the 
escaping  fluorine. 

o  amount  of  air  mixed  with  the  chemical,  temperature,  and 
relative  humidity  of  the  ambient  air. 

o  energy  input  from  the  ambient  into  the  fluorine-air 
mixture. 

We  employ  the  information  opacified  above,  with  appropriate 
equations  representing  the  enthalpies  of  fluorine  (vapor  and 
liquid),  water  (vapor,  liquid,  and  solid),  air  (vapor)  as  functions 
of  temperature,  to  calculate  the  "initial  stream  enthalpy."  This 
term  is  called  H,  and  has  the  dimensions  of  joules.  No  mixing  has 
as  yet  occurred  ao  the  property  values  of  pure  fluorine  are  used. 

We  add  to  Hj  the  given  energy  input,  Q,  also  in  joules.  The  sum 
H,  +  Q  represents  the  total  stream  enthalpy  at  any  time.  Because 
of  conservation  of  energy,  the  final  enthalpy  after  mixing,  Hr,  is 
set  equal  to  the  initial  stream  enthalpy  plus  any  external  heat 
added  Q. 

When  carrying  out  a  calculation,  one  does  not  know  at  the  start 
which  condensed  phases  are  present  -  and  of  those  which  are,  the 
quantities  of  eaoh.  To  overcome  this  difficulty,  we  first  bracket 
the  domains  and  compute  Hf  values  at  the  extreme  conditions.  This 
technique  readily  allows  one  to  determine  the  temperature  range 
within  which  the  final  mixture  temperature  will  occur.  Thu  final 
value  can  then  be  obtained  quickly  by  an  iterative  procedure 
locating  the  final  state  which  satisfies  the  equation  H,  *  Hi  +  Q, 
and  the  condensed  phase  equilibrium  relationships  (if  condensed 
phases  are  found  to  exist) . 

We  have  identified  five  principal  domains  of  the  final  mixture 


conditions, 

as  shown  below? 

Cfim&ln 

URL 

Phases* 

i 

>273.2  to  373** 

F2(gas);  H}0(gas,  liquid) 

ii 

«  273.2 

F](gas);  H30(gas,  liquid,  and 
solid) 

hi 

84.5  to  273.2 

Fj(gaa)  ;  H:0(gas,  solid) 

IV 

■  84.5 

Fa(gas,  liquid);  HjO(solid) 

V 

<  84.5 

Fj(gas,  liquid);  HjO(eolid) 

*  Air  is  always  present  as  a  gas 

**  Chosen  arbitrarily  as  the  maximum  possible  temperature 
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Thus  we  begin  with  Domain  I  and  calculate  Hf  at  both  extremes.  If 
II,  +  Q  lies  between  these  limits,  we  know  the  solution  lies  in 
Domain  I.  If  Hj  +  Q  <  Hf  at  the  lower  temperature  end  of  Domain  I, 
then  we  repeat  the  procedure  for  Domain  II,  etc.  In  a  few  steps, 
we  can  locate  the  appropriate  domain  and  then  direct  our  attention 
to  the  computation  of  the  final  mixed  state. 


As  noted  earlier,  if  there  are  fluorine  liquid  aerosols  in  the 
release,  the  temperature  of  the  cloud  will  drop  upon  mixing  with 
air  as  the  liquid  fluorine  evaporates.  In  developing  and  testing 
the  model,  one  must  make  certain  that  the  decrease  in  temperature 
docs  not  cause  a  condensation  of  oxygen-enriched  air.  Pure  air 
begins  to  condense  at  about  79  K.  In  the  early  phases  of  mixing, 
the  partial  pressure  of  air  is  below  one  bar  so  the  condensation 
temperature  is  below  79  X.  The  following  table  indicates  the 
condensation  temperature  of  air  as  a  function  of  the  partial 
pressure  of  air. 


Air  Condensation 


Air  Partial 


Temperature 

(K) 


Pressure 

(bar) 


78 

75 

70 

65 

64 


0.91 

0.64 

0.32 

0.15 

0.12 


Only  by  carrying  out  sets  of  example  calculations  can  one  determine 
if  it  is  possible  to  condense  air.  If  so,  this  step  would  need  to 
be  included  in  the  enthalpy  balance.  Preliminary  calculations  do 
not  suggest  air  condensation  is  a  problem  unless  the  fraction 
liquid  fluorine  is  high. 

Hence,  the  important  step  in  determining  the  state  of  the  final 
mixture  (and  hence,  the  mixture  thermodynamic  parameter  values)  is 
the  evaluation  of  the  final  mixture  enthalpy  Hf.  There  are  two 
approaches  to  performing  this  calculation.  In  the  first  method, 
the  possible  chemical  reaction  between  fluorine  and  water  vapor  is 
neglected.  In  the  second  method,  the  complexity  of  fluorine  - 
water  vapor  reaction  is  considered.  This  reaction  is  indicated  in 
equation  5.1  and  results  in  the  formation  of  Hydrogen  fluoride  (HF) 
and  oxygen  (02)  . 

In  considering  the  fluorine  -  water  vapor  reaction,  all  the  water 
in  the  air  is  allowed  to  react  with  fluorine  upon  mixing  with  air. 
Here,  the  energy  input  is  8  MJ/kg  of  F2  reacting  or  17  MJ/kg  of 
water  consumed.  Then,  one  would  have  no  further  water,  but  the 
enthalpies  of  the  HF  and  02  products  would  have  to  be  included.  As 
an  approximation,  the  masses  of  02  and  HF  formed  are  added  to  the 
air  present  and  not  considered  as  individual  species.  For  each 
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mole  of  water  reacted,  2  moles  of  HF  and  one-half  mole  of  02  are 
produced  -  or  per  kg  of  water  reacted,  0.89  kg  02  and  2.22  kg  HF. 

The  rates  of  reaction  of  fluorine  and  water  vapor  are,  however, 
expected  to  be  negligible  at  the  low  temperatures  prevailing  in  a 
liquid  fluorine  release.  There  are,  in  fact,  no  reaction  kinetic 
data  for  fluorine  and  water  vapor  reaction  at  temperatures  close  to 
the  saturation  temperature  of  fluorine.  When  the  temperature  has 
risen  to  anywhere  near  ambient,  the  fraction  fluorine  in  the  vapor 
cloud  is  quite  low  so,  again,  reactions  would  not  be  important. 
Hence,  the  neglect  of  the  heats  of  fluorine  -  water  vapor  reaction 
in  determining  the  final  thermodynamic  state  of  the  mixture  would 
not  in  any  way  reduce  the  correctness  of  the  calculations. 


5.4.3  Results  from  the  Fluorine  -  Humid  Air  Nixing  Model 

The  results  obtained  from  the  above  model  are  indicated  in 
Figure  5.7  and  Figure  5.8.  The  variation  of  the  final  mixture 
density  and  the  mixture  temperature  with  the  amount  of  air  added 
are  indicated  in  Figure  5.7.  The  temperature  of  air  is  set  at 
293  K  and  the  relative  humidity  value  is  a  variable.  The  mass  of 
liquid  present  in  the  mixture  in  the  form  of  aerosol  is  plotted  in 
Figure  5.8  for  various  ambient  relative  humidities  and  dilution 
ratios.  These  plots  are  for  a  given  initial  composition  of  the 
saturated  fluorine  stream  released.  In  the  example  shown,  this  is 
assumed  to  be  the  condition  after  flash  occurring  as  a  result  of 
tank  burst  (the  stream  contains  13.61  %  by  mass  of  saturated 
vapor) . 

It  is  seen  from  Figure  5.7  that  the  mixture  density  is  very  high 
(35  kg/m3)  at  low  dilution  ratios.  This  is  because  we  have  assumed 
that  all  of  the  liquid  released  after  flash  is  entrained  into  the 
vapor  cloud  in  the  form  of  liquid  aerosols.  The  density  decreases 
continuously  with  increase  in  dilution.  Also,  the  variation  in  the 
relative  humidity  of  the  ambient  has  very  little  effect  on  the 
final  mixture  density  at  any  specified  dilution. 

The  temperature  variation  is  also  shown  in  Figure  5.7.  The 
temperature  decreases,  albeit  very  slightly  (about  4  or  5  K)  for 
dilution  ratio  below  0.6.  This  temperature  decrease  is  due  to  the 
evaporation  of  the  liquid  F,  aerosols  at  the  expense  of  the 
sensible  heat  of  the  system.  The  enthalpy  provided  by  the  mixing 
in  air  is  insufficient  to  raise  the  temperature.  The  effect  of 
partial  pressure  decrease  (of  F.)  is  dominating  this  early  process. 
However,  as  more  and  more  air  is  mixed,  the  liquid  aerosols 
evaporate  at  a  dilution  ratio  of  about  0.6.  Any  additional  air 
added  increases  the  temperature.  The  effect  of  water  is  not 
noticeable  in  small  dilutions.  However,  at  larger  dilutions,  the 
heat  of  condensation  of  water  provides  tho  thermal  energy  to  raise 
the  temperature  of  the  mixture  (compared  to  the  dry  air  case  at  the 
same  dilution) . 
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atr/ltas  of  F2) 


The  variation  of  the  liquid  mass  in  the  final  mixture  is 
illustrated  in  Figure  5.8.  As  can  be  seen,  the  initial  0.864  kg/kg 
of  F,  liquid  rapidly  evaporates  and  the  liquid  is  essentially 
completely  vaporised  by  the  time  the  dilution  ratio  is  0.6  of  0.7. 
However,  the  water  at  this  stage  is  in  the  form  of  ice  in  the 
mixture.  As  the  dilution  increases  beyond  the  above  0.6  or  0.7 
value,  there  is  no  liquid  in  the  system  until  the  water  ice  melts 
to  liquid  water.  This  occurs  at  a  dilution  ratio  of  about  17  in 
the  case  of  100  %  relative  humidity  air  and  at  about  150  in  the 
case  of  50  %  relative  humidity.  Mote  also  the  re-evaporation  of 
the  water  liquid  from  the  mixture  as  the  dilution  is  greater  than 
200.  This  is  because  at  a  dilution  of  200  the  temperature  of 
mixture  is  just  above  the  wet  bulb  temperature  in  air  at  the  50% 
relative  humidity.  Table  5.4  gives  the  detailed  numerical  values 
for  the  various  thermodynamic  parameters  of  the  final  mixture  of  F2 
and  humid  air  for  different  dilution  ratios  and  ambient  relative 
humidities. 


5.5  LN2  and  LF2  Pool  Boiling  Model  (SCENARIO  B) 

5.5.1  Source  Rate  Calculations 

In  section  5.3.2,  we  described  the  second  type  of  release  accident 
which  can  result  in  the  formation  of  a  liquid  pool  on  the  ground. 
The  pool  will  contain  both  liquid  nitrogen  and  liquid  fluorine. 
These  two  liquids  being  very  cold  will  rapidly  boil  on  the  ground 
surface,  evaporate  and  form  a  vapor  cloud  consisting  of  cold 
fluorine  and  nitrogen  vapors.  The  evolved  vapor  is  carried 
downwind  and  mixed  with  ambient  air.  No  liquid  fluorine  or  liquid 
nitrogen  aerosols  are  present  in  the  source  vapor.  This  boiling 
phenomenon  is  discussed  in  this  section. 

The  detailed  mathematical  representation  of  the  phenomenon  of 
boiling  and  distillation  of  a  pool  containing  two  liquid  components 
is  discussed  in  Appendix  B.  When  a  release  of  both  LN2  and  LF2 
occurs  from  the  road  transport  as  a  result  of  an  accident, 
initially  a  pool  of  liquid  is  formed  on  the  ground.  The  extent  of 
spread  of  this  pool  on  the  ground  depends  on  the  local  topography, 
nature  of  the  soil,  slope  in  the  ground  and  whether  there  are  any 
containment  depressions  on  the  ground.  In  the  model  discussed  in 
Appendix  B,  we  assume  that  the  pool  spreads  to  an  area 
corresponding  to  an  initial  depth  of  l  cm. 

This  liquid  pool  evaporates  due  to  the  transfer  of  heat  from  the 
ground.  The  heat  transfer  from  the  ground  to  the  liquids  is 
modeled  using  the  one  dimensional  heat  transfer  model  (discussed  in 
a  previous  report  by  Raj  and  Morris,  1987).  Because  the  pool 
contains  two  liquids  with  differing  boiling  temperatures  (77  K  for 
N2  and  84.5  K  for  F2) ,  the  pool  boiling  process  is  similar  to  a 
process  of  distillation  of  a  two  component  mixture.  In  the  case  of 
a  LN2  -  LF2  pool,  nitrogen  evaporates  preferentially  during  the 
initial  stages.  Therefore,  in  the  beginning  the  vapor  evolved  is 
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TABLE  5.4:  RESULTS  OF  CALCULATIONS 
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richer  in  nitrogen  than  in  fluorine.  However,  at  later  stages  the 
vapor  consists  essentially  of  pure  fluorine. 

The  important  parameters  of  interest  to  the  dispersion  process  are 
the  rate  of  generation  of  vapor  and  its  quality  (or  concentration 
of  F2 ) .  Also  of  interest  is  the  total  duration  over  which  the  cold 
liquids  evaporate  on  the  ground. 

Detailed  results  from  this  vapor  source  rate  model  are  presented  in 
Appendix  B.  Figure  5.9  shows  the  variation  with  time  of  the 
concentration  of  fluorine  and  nitrogen  in  the  vapor  evolved.  It  is 
clearly  seen  that  the  vapor  progressively  becomes  richer  as  time 
goes  on.  At  first  glance,  this  may  seem  to  imply  that  the  downwind 
dispersion  distance  for  fluorine  (for  specified  hazard 
concentration)  will  increase  with  time.  However,  the  rate  of 
evaporation  decreases  with  time  resulting  in  lesser  mass  rates  of 
fluorine  injection  into  the  atmosphere.  In  Figure  5.10,  the  mass 
rate  of  evaporation  of  fluorine  and  nitrogen  are  indicated.  As  can 
be  seen,  the  rates  of  evaporation  of  the  two  components  are  close 
and  have  the  same  trend  with  time. 

One  important  result  that  is  noticed  from  Figure  5.9  is  that  all 
nitrogen  and  fluorine  are  evaporated  over  a  relatively  short  period 
of  time,  namely,  about  14  seconds  1  This  is  due  to  the  limited 
volume  of  liquids  spilled  (several  hundred  gallons)  and  the 
extremely  high  heat  transfer  rates  from  the  ground. 

5.5.2  Consideration  of  source  Hate  in  Dispersion  Model 

It  is  anticipated  that  the  above  duration  for  pool  evaporation 
(i.e.,  the  duration  over  which  the  vapors  are  released)  will  be 
very  small  compared  to  the  potential  dispersion  duration  to  dilute 
to  concentrations  below  the  F2  hazard  concentrations.  Therefore, 
the  dispersion  calculations  can  be  performed  in  one  of  two  ways, 
namely,  by  assuming  that: 

1.  all  of  the  F2  vapor  is  released  instantaneously.  in  this 
case,  an  initial  cloud  of  cold  F2  and  N2  vapors  is  assumed  to 
be  released.  The  dispersion  i3  then  modeled  as  an  initially 
diluted  puff. 

2.  F2  vapor  is  released  at  a  constant  rate  (corresponding  to  the 
mean  rate)  along  with  nitrogen  vapor.  The  concentration  of  F2 
vapor  in  this  steady  stream  of  vapor  is  assumed  to  be  at  the 
mean  value. 

It  is  our  contention  that  the  former  assumption  makes  more  sense 
and  provides  a  more  conservative  estimate  of  the  extent  of  hazard 
posed  by  a  F2  release  from  a  transport.  Hence,  in  ADAM,  we  have 
coded  the  Scenario  B  discussed  earlier  (see  section  5.5.1)  as  an 
instantaneous  releaFe  of  a  f2  -  N*  cloud.  The  rationale  for  this 
assumption  is  that  any  accident  which  is  severe  enough  to  puncture 
both  the  outer  nitrogen  tank  and  the  inner  fluorine  tank  will 
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both  the  outer  nitrogan  tank  and  the  inner  fluorine  tank  will 
oreate  large  holes  in  both  tanks  leading  to  a  rapid  draining  of  the 
tank  contents.  The  evaporation  of  the  liquids  on  the  ground  will 
be  very  rapid  as  has  been  shown  in  the  previous  section. 

Therefore,  the  fluorine  vapor  source  for  the  second  scenario  of 
release  (Scenario  B)  is  specified  as  a  cylindrical  vapor  cloud 
consisting  only  of  cold  vapors  of  fluorine  and  nitrogen,  completely 
mixed,  and  containing  no  liquid  aerosols.  This  heavy  gas  cloud  is 
than  allowed  to  disperse  in  the  prevailing  atmospheric  condition. 
The  dispersion  routines  in  ADAM  are  then  utilized. 


S.f  Discussion  of  Results  from  the  fluorine  Dispersion 
Calculations 

Calculations  were  performed  for  the  case  of  fluorine  release  from 
a  road  transport  vehicle  carrying  the  chemical.  The  details  of  the 
volumes  of  liquid  fluorine  and  liquid  nitrogen  carried  were 
indicated  in  Table  5.1.  These  values  are  used  in  the  dispersion 
calculations  discussed  in  this  section.  Both  soenario  A  type 
release  and  scenario  B  type  of  release  have  been  ooneidered. 

In  the  case  of  scenario  A  type  release  (occurring  as  a  result  of 
nitrogen  loss  and  inner  fluorine  tank  failure  due  to  overpressure)  , 
fluorine  released  flashes.  It  is  calculated  that  13.6%  of  the 
liquid  flashes  to  vapor  and  the  remaining  86.4%  of  the  mass 
manifests  itself  as  liquid.  We  have  aeaumed  that  all  of  this 
liquid  will  be  entrained  into  the  vapor  as  aerosols.  We  also 
assume  that  a  mass  of  air  equal  to  10  times  the  mass  of  air 
contained  within  a  volume  occupied  by  the  vapor  released  from  the 
explosive  failure  of  the  tank.  Two  atmospheric  stability 
conditions,  namely,  D  and  F  stabilities  were  used  in  the 
calculations. 

Figure  5.11  shows  the  results  for  the  case  of  fluorine  release  from 
the  bursting  of  the  tank  due  to  overpressure  under  neutral 
stability  (D  stability)  conditions  with  wind  speed  of  5  m/s.  The 
contour  plotted  is  the  l  ppm  contour  which  is  the  short  term 
exposure  (toxic)  concentration  limi "■  for  fluorine.  As  can  be  seen 
that  because  of  the  very  low  hazard  concentration  for  F3i  the  hazard 
distance  extends  to  significant  distances,  namely,  about  21  kml 

t 

The  dispersion  of  the  same  amount  of  fluorine  released  in  a  stable 
F  type  atmosphere  is  indicated  in  Figure  5.12.  The  maximum  hazard 
in  this  case  extends  to  an  even  longer  distance  of  about  64  kml 

Because  of  the  low  concentration  to  which  the  cloud  has  to  be 
diluted,  it  is  estimated  that  a  substantial  part  of  the  dispersion 
takes  place  in  the  near  neutral  density  dispersion  regime.  As  a 
matter  of  fact,  the  hazard  distance  can  be  calculated  with  neutral 
density,  point  source  release  Gaussian  model  and  come  with 
substantially  the  same  result. 
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SCENARIO  Ai  TANK  BURST  STABILITY  D  (4.0) 


Grid  radii  ■  6000,00  m 
1  Inch  >  7921.85  * 

miss  <wih>  to  cqmtinui 


Figure  5.11i  Dispersion  under  Neutral  Stability  Fluorine 
Cloud  Released  by  a  Transport  Tank  Burst 

SCENARIO  Al  TANK  BURST  STABILITY  F  (5.5) 


HWH 

Uapor  Diapers Ion  Hods I 


CHMICftL  RELEASED  I 
FLUORINE 

Source  nasi  2228.81  kg 
Cone  contour  1,0  ppM 
Wind  speed  (10  n)  5.0  M/s 
Diaper  duration  3513,0  a 
Hazard  area  129631100,  ■' 
flux  haetrd  diet  61117.9  n 


Uapor  Dispersion  Node  I 

CIIEMICftL  RELEASED: 
FLUORINE 

Source  ness  2225,81  ky 
Cono  contour  1,0  ppn 
Wind  speed  (10  ■)  5,0  m/i 
Dliper  duration  3517.7  i 


Ha sard  area  333S0S90,  «« 
Max  hastrd  diet  20890. 9 


Grid  radii  ■  ooeo.oo  m 
1  Inch  ■  7013,27  m 

FR1SS  <EMTUI>  TO  CONTINUE 


Figure  5 . 12 1  Dispernion  under  Stable  Weather  Conditions 
of  a  Fluorine  Vapor  Cloud  Released  by  a 
Transport  Tank  Burst 
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The  results  for  the  case  of  release  of  fluorina  and  nitrogen 
liquids  and  the  dispersion  of  the  vapors  thus  produced  are  shown 
respectively*  for  0  stability  condition  and  F  stability  condition* 
in  Figure  5.13  and  Figure  5.14.  Again*  it  is  noticed  that  the 
distances  are  comparable  to  the  previous  numbers  obtained  in  the 
case  of  tank  burst.  This  is  because  of  two  reasons.  The  first  is 
that  the  total  mass  of  fluorine  released  is  the  same  in  both 
scenarios.  Secondly*  the  dilution  that  is  to  be  achieved  to  reach 
1  ppm  concentration  is  it  io4.  This  is  a  very  large  dilution  and 
the  contribution  from  the  initial  dilution  from  the  cold  nitrogen 
vapors  is  relatively  small. 
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SCENARIO  B:  TANK  BREACH i  STABILITY  D  (4.0) 


Figure  5.13:  Dispersion  under  Neutral  Stability  Fluorine 
Cloud  ReieRHed  by  a  Transport  Tank  Breach 


SCENARIO  Bi  TANK  BREACH t  STABILITY  P  (5.5) 


Figure  5.14;  Dispersion  under  Stable  Weather  Conditions 
of  a  Fluorine  Vapor  Cloud  Released  by  a 
Transport  Tank  Breach 
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CHAPTER  6 


CONCLUOXONS 


The  scope  of  Air  Force  Dispersion  Assessment  Model  (ADAM)  has  been 
expanded  to  include  the  simulation  of  accidental  releasee  of  two 
additional  chemicals,  namely,  anhydrous  hydrogen  fluoride  (HF)  and 
fluorine.  The  analyses  presented  in  this  report  include  modeling 
both  the  release  scenario  ("source  modeling")  and  dispersion  of  the 
chemical  plume  or  cloud  formed  after  release. 

Inclusion  of  HF  in  the  list  of  chemicals  in  ADAM  required  the 
evaluation  of  available  HF-humid  air  reaction  thermodynamic  models, 
modifications  and  improvements.  The  dispersion  algorithm  for  HF  is 
the  same  as  the  one  for  other  chemicals,  namely  a  hybrid  of  a  heavy 
gas  (modified)  box  model  and  a  modified  finite  size  source 
Gaussian. 

Release  of  fluorine  from  a  transport  vehicle  was  analyzed.  Two 
different  scenarios  of  fluorine  release  were  identified.  Detailed 
source  models  have  been  developed  for  the  evaporation  of  a  two 
component  liquid  pool  (consisting  of  liquid  nitrogen  and  liquid 
fluorine)  and  the  flashing  release  of  pressurized  liquid  fluorine. 
Also,  a  thermodynamic  model  was  developed  to  calculate  the  effect 
of  mixing  of  fluorine  vapor  (or  vapor  with  aerosols)  and  humid  air. 

Based  on  the  work  accomplished  in  this  project,  we  conclude  the 
following: 


Hydrogen  Fluoride 

1.  The  HF-humid  air  model  is  a  modified  Schotts  model.  It  has 
been  modified  to  expand  its  range  of  applicability  and  improve 
its  acauraoy.  The  model  is  robust  and  calculates  the  final 
HF-air  mixture  conditions  for  a  variety  of  initial  conditions 
of  air  and  HF. 

2.  The  source  and  dispersion  calculations  give  results  close  to 
the  values  observed  in  one  series  of  field  tests  with  HF. 

3.  There  are  significant  uncertainties  in  the  field  test  data. 
Many  questions  regarding  instrument  calibration,  conversion  of 
raw  field  data  to  concentration  values,  discrepancies  in 
reported  values  for  the  same  data  items  in  two  publications 
remain  unresolved. 
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4.  The  concentration  values  predicted  by  ADAM  seem  to  be  lower 
than  observed  values.  Also  at  the  transition  regime  (between 
heavy  gas  to  passive  dispersion) ,  the  concentration  vs. 
distance  curve  seems  to  have  a  different  trend  than  what  one 
expects.  This  aspect  needs  a  significantly  more  detailed 
model  parameter  analysis. 

5.  The  HF  model  in  ADAM  predicts  results  that  have  accuracy 
acceptable  for  a  rapid  assessment  of  a  potentially  hazardous 
situation  in  an  accident  environment. 


Fluorine 

6.  Two  different  fluorine  release  conditions  can  occur  depending 
on  the  accident  scenario.  These  are,  (i)  the  simultaneous 
release  of  liquid  nitrogen  (from  the  jacket)  and  liquid 
fluorine  at  the  saturation  temperature  of  nitrogen  and  (ii) 
the  flashing  release  of  fluorine  (only)  due  to  inner  tank 
rupture  caused  by  elevated  pressure. 

7.  Two  component  (LN2  and  LF2)  liquid  pool  evaporation  model 
developed  indicates  that  nitrogen  evaporates  preferentially. 
However,  because  of  the  very  low  boiling  temperature  of  both 
nitrogen  and  fluorine,  and  the  consequent  high  heat  transfer 
rates  from  ground  to  the  liquid  pool,  the  evaporation  of 
fluorine  is  very  rapid  (within  a  matter  of  tens  of  seconds) . 

8.  In  the  flashing  release  of  fluorine  about  13.5%  of  the  mass 
released  flashes  to  vapor  instantaneously.  The  cloud  formed  by 
the  mixture  of  the  liquid  aerosol  and  the  saturated  fluorine 
vapor  has  very  high  density  relative  to  that  of  air 
(32.8  to  1) . 

9.  Mixing  of  cold  fluorine  vapor  or  a  mixture  of  vapor  and  liquid 
aerosol  with  humid  air  results  in  five  different  final 
conditions  of  the  mixture  depending  upon  the  amount  of  air 
mixed  and  the  relative  humidity  of  air.  The  final  conditions 
refer  to  the  presence  of  solids  (water  ice) ,  liquid  aerosols 
(water  and  fluorine)  and  all  components  in  the  vapor  phase. 

10.  The  results  of  the  fluorine  -  air  mixing  model  indicate  that 
the  final  mixture  density  is  relatively  insensitive  to  the 
relative  humidity  of  air  but  dependents  on  the  amount  of  air 
mixed  (for  specified  initial  conditions  of  fluorine). 

11.  The  temperature  of  the  mixture  is  relatively  constant  until 
about  equal  mass  of  air  is  mixed  with  fluorine.  Beyond  this 
dilution  the  mixture  temperature  rises  and  the  rate  of  rise  is 
dependent  on  the  air  relative  humidity. 
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12.  The  dispersion  distances  calculated  for  the  case  of  flashing 
fluorine  release  from  a  road  transport  for  1  ppm  level  of 
hazard  are  21  km  and  64  km,  respectively  for  neutral  and 
stable  atmospheric  conditions.  These  large  distances  arise  due 
to  the  very  low  concentration  levels  to  which  fluorine  has  to 
be  diluted. 

13.  Even  in  the  case  of  liquid  pool  formation  and  evaporation  the 
calculated  hazard  distances  are  comparable  to  the  above.  This 
is  because  the  total  mass  of  fluorine  released  is  the  same  and 
the  pool  evaporates  in  relatively  short  time.  Therefore, 
during  most  of  the  dispersion  regime  the  cloud  is  dispersing 
as  if  it  were  released  instantaneously. 

14.  Because  of  the  very  low  hazard  concentration  of  fluorine  the 
extent  of  dilution  to  safety  is  of  the  order  of  10*.  By  100  to 
1  dilution  the  fluorine  air  mixture  h^s  a  density  that  is 
within  0.5  t  of  that  of  ambient  air.  That  is,  in  the  next  104 
dilution  the  cloud  disperses  essentially  as  a  neutral  density 
vapor . 

Hence,  the  overall  effects  of  the  higher  than  air  density  of 
the  vapor-aerosol  mixture  affect  the  dispersion  process  only 
during  the  initial  stages;  they  have  very  little  effect  on  the 
final  outcome  of  the  dispersion  model.  Ordinary  Gaussian 
models  can  be  used  therfore  for  evaluating  the  hazard 
distances  for  fluorine  spills.  However,  if  near  field 
concentration  values  or  cloud  widths  are  needed  a  heavy  gas 
model  such  as  ADAM  has  to  be  used. 

15.  There  are  no  experimental  data  on  fluorine  releases  with  which 
to  compare  the  ADAM  dispersion  predictions. 
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APPENDIX  A 


FLASHING  OF  ANHYDROUS  HYDROGEN  FLUORIDE 

RELEASED  FROM  A  COMPRESSED  LIQUID  STATE 

The  normal  boiling  point  of  anhydrous  hydrogen  fluoride  (HF)  at 
atmospheric  pressure  is  20  °C.  In  the  case  of  HF  stored  at  ambient 
temperature  but  at  an  elevated  pressure,  the  release  of  liquid 
results  in  the  formation  of  vapor  due  to  flash  vaporization  of  the 
liquid  released.  The  pressure  Inside  the  tank  may  be  higher  than 
the  ambient  pressure  either  due  to  the  ambient  temperature  being 
higher  than  the  saturation  temperature  (20  °C)  corresponding  to  l 
atmospheric  pressure  or  due  to  the  pressurization  of  the  tank  by 
another  inert  gas  (such  as  nitrogen). 

In  this  appendix  we  discuss  the  calculation  of  the  percent  flash 
into  vapor  for  HF  releases  from  under  higher  than  saturation 
pressure. 


A.  It  FLASH  FRACTION 

The  mass  fraction  of  released  liquid  HF  which  flash  to  vapor  is 
given  by, 


fV  “  [  ^Uq(T-f)  -  hH,(PJ  ]/MPJ  (A.i) 

where, 

hUq(TT)  "  Enthalpy  of  liquid  at  temperature  of  tank  TT 

hltq(Pltm)  -  Enthalpy  of  saturated  liquid  at  ambient  pressure 

X (P-tM)  *  Latent  heat  of  vaporization  at  ambient  pressure 

Tt  -  Liquid  storage  temperature  in  the  tank 

The  key  to  determining  the  flash  fraction  is  to  determine  the 
liquid  enthalpy  at  the  tank  temperature.  In  order  to  calculate  this 
it  is  necessary  to  calculate  the  saturation  pressure  of  HF 
corresponding  to  the  tank  storage  temperature. 

Let, 


PT  ■  Tank  pressure 

P#,t(TT)  ■  Saturation  pressure  corresponding  to  Tt 
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Depending  on  the  relative  magnitudes  of  the  above  two  pressures 
there  are  three  cases  to  consider  in  evaluating  the  liquid 
enthalpy. 


Case  is  PT  <-  1.01  *  P"*(Tt) 

That  is  the  tank  pressure  is  within  1  %  of  the  saturation  pressure 
corresponding  to  the  tank  temperature.  In  this  case  we  oan 
approximate  the  liquid  enthalpy  as  follows! 


hUq(TT)  -  hltq(sat,  T7)  -  hUq(pMt(Tr))  (A. 2) 

Case  2:  PT  >  PMt(TT) 

The  liquid  in  the  tank  under  this  condition  is  in  a  compressed 
liquid  state.  Hence,  the  enthalpy  is  calculated  by  the  formula 

hl(p(TT)  -  hUq(sat,TT)  +  (fih/fiP)T  dP  (A. 3) 


where, 


(fih/AP)T  -  [  -  1  +  M  ]  /  Puq  <A-4) 


(1  -  -  (1/p)  (  ftp  )  (A. 5) 


Coefficient  of  volumetric  expansion  of  liquid 
with  temperature  under  constant  pressure. 


The  second  term  on  the  right  hand  side  of  equation  A. 3  (i.e.,  the 
integral)  represents  the  contribution  to  the  liquid  enthalpy  due  to 
liquid  compression  over  the  saturation  pressure.  This  term  is 
represented  by  the  symbol  "Ah". 


The  calculation  of  0  value  for  HF  liquid  is  indicated  in  a  section 
A. 2. 


Case  3:  PT  <  PMt(TT) 

This  case  implies  that  there  is  vacuum  or  that  the  liquid  is  not  at 
saturated  condition.  That  is,  the  liquid  is  superheated  with 
respect  to  the  tank  pressure.  This  condition  is  not  of  interest  to 
us. 


A  -2 


A.  2:  NUMERICAL  VALUES 


The  values  for  the  HF  liquid  enthalpies  for  non  saturated 
conditions  are  not  available  in  the  literature.  Also  not  available 
is  the  value  for  the  coefficient  of  liquid  expansion.  Therefore,  we 
use  the  indirect  and  approximate  method  to  determine  the  values  for 
the  non  saturated  liquid  enthalpy. 

The  following  values  are  indicated  in  handbooks  for  coefficient  of 
liquid  expansion,  p. 

P (hydrochloric  acid)  -  0.49  x  10'3  K'1,  Refs  Handbook  of 

Mech  Engineering,  p 
4-9,  Table  12. 

p (organic  liquids)  -  0.043 14/ (Tc  -  T)0-6*1  K1,  Refj  Chemioal 

Engineers 
Handbook,  p  3- 
227. 

We  assume  that  the  p  for  HF  is  close  to  that  of  hydrochloric  acid 
and  use  that  value  given  above. 

Density  of  HF  liquid  at  300  K  (p)  -  940  kg/m3 

Henoe,  using  equation  A. 4  and  the  above  values  for  p  and  p  we  can 
show  that 


(Ah/6P)T  -  -  9  x  10'*  (J/kg)  /  (N/m*)  at  300  K 

Also,  we  have  the  following  thermodynamic  values  for  the  HF 
conditions  at  313  K  the  tank  temperature. 


PMt  (313) 

-  2.02  X  10* 

N/m* 

hllq(sat,313) 

-  7.11  X  10* 

J/kg, 

with  the  enthalpy 
base  -  0  at  292.7  K 

Pt 

-  8.66  X  10* 

N/m* 

-  Tank  pressure 

-  3.744  X  103 

J/kg 

Hence,  we  get  for  Ah,  the  contribution  to  liquid  enthalpy  due  to 
liquid  compression, 

Ah  *■  -  598  J/kg 


It  is  noticed  by  comparing  the  saturated  liquid  enthalpy  at  the 
tank  temperature  and  the  value  of  Ah  above  that  the  latter  forms  a 
very  small  fraction  of  the  liquid  enthalpy.  Hence,  no  great  loss  of 
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accuracy  rasults  if  the  saturated  liquid  enthalpy  value  is  at  the 
tank  temperature  is  used  for  the  actual  liquid  enthalpy. 


A. 3!  FLASH  VALUES 

The  saturated  liquid  enthalpy  value  for  anhydrous  HF  can  be 
represented  by  the  equation, 

hUq(sat,T)  -  3500  (T  -  292.7)  J/kg  (A. 6) 

where  T  is  the  temperature  of  the  liquid  in  degrees  Kelvin. 

Using  equation  A. 6  and  equation  A.l  we  can  show  that  the  value  of 
vapor  mass  fraction  resulting  from  the  flash  is, 

fv  (in  I)  -  0.935  (T  -  292.7)  (A. 7) 

The  fraction  that  flashes  to  vapor  for  liquid  release  from  313  K  is 
therefore  calculated  to  be  19.12  I. 
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APPENDIX  B 


DISTILLATION  OP  A  TWO  COMPONENT  SYSTEM 


1.0  INTRODUCTION 

In  this  appendix,  wa  consider  the  distillation  of  a  pool  of 
cryogenic  liquid  containing  two  volatile  components,  namely, 
fluorine  and  nitrogen.  The  pool  of  this  cold  liquid  is  on  .the 
ground  and  evaporating  due  to  heat  transfer  from  the  ground. 

The  important  objective  of  the  calculations  presented  in  this 
appendix  is  to  determine  the  total  time  for  evaporation  of  the 
pool.  It  ie  also  desired  to  calculate  the  rate  of  evaporation  of 
the  liquids,  the  molar  concentrations  of  the  two  components  in  the 
liquid  pool  and  the  vapor  generated,  as  functions  of  time. 

The  details  of  the  model  formulation,  equations,  the  computer  oode 
written  in  FORTRAN  to  solve  the  ooupled  differential  equations  and 
the  results  calculated  are  indicated  in  this  appendix. 


2.0  THE  EVAPORATION  MODEL 

The  two  component  liquid  mixture  evaporating  from  a  pool  is  assumed 
to  be  located  on  the  ground.  The  energy  for  vaporisation  is 
assumed  to  come  from  the  ground.  The  components  are  represented  by 
subscripts  1  and  2.  In  our  nomenclature,  component  1  refers  to 
Fluorine  and  component  2  refers  to  Nitrogen.  (The  definition  of 
the  symbols  used  in  this  model  formulation  are  given  at  the  and  of 
this  appendix.) 

We  now  define  the  molar  contents  of  the  liquid  pool  as  follows! 

L,  -  X,L  (la) 

Lj  «  XaL  (lb) 

The  heat  from  the  ground  evaporates  both  components  of  the  liquid. 
For  an  incremental  input  of  heat,  the  incremental  change  in  the 
molar  oontents  of  the  pool  are  given  by, 


dQ  -  (h1dL1  +  hjdLj)  (2) 

Heat  input  to  pool 


The  negative  sign  indicates  that  as  heat  is  added  to  the  pool  the 
moles  in  the  liquid  decrease.  The  molar  heats  of  vaporization  are 
indicated  by  h,  and  h2. 
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Using  Equations  Is  and  lb,  Equation  2  is  raduoad  tot 


do  -  -  [  L  i  hxdXx  +  ha  dXx)  +  { Ja^  ♦  haXa I  dL  ] 


Wa  also  nota  that 


X,  +  X2 


Hanoa,  aubstituting  Equation  4  in  Equation  3,  wa  gat 
do  -  -  [L  Oji-h,!  dXx  ♦  +  Jja<i  -  ^))  dL] 


-  -  lx,thx-/aa)-g^-  +  DijJfj  hjd-Xj) 


If  a  ralationship  is  astablishad  batwaan  L  and  Xir  than,  for 
spaoif iad  haat  input  rata  (LHS  of  Equation  Sb) ,  wa  oan  aolva  for 
tha  tamporal  variation  of  X,  (and  hanoa,  L  and  X.) .  Tha  aaotion 
balow  illuatratas  tha  davalopmant  of  tha  ralationship  batwaan  X. 
and  L. 


3.0  RELATIONSHIP  BETWEEN  X,  AND  L 

Tha  avaporation  of  a  liquid  oomponant  rsaults  in  tha  inoraasa  in 
tha  mola  oontant  of  tha  vapor.  Using  a  mass  balanoa, 


dL  +  dV 


-  0 


d(X,L)  +  Y,  dV  -  0  (6b) 

whara,  V1  raprasants  tha  mola  fraotion  of  oomponant  l  in  tha  vapor 
ganaratad  and  Y,  dV  is  tha  total  molaa  of  oomponant  1  ganaratad  in 
an  infinitasimal  tima. 


d  (XxL) 

— ar“  " + 


i.s 


•  i 


dX±  dL 

Iyx-X;)  “  T 


(7  b) 


As  V,  snd  X.  ars  rslatsd  by  phasa  squilibria,  Equation  7b  oan  ba 
intagratad. 

Wo  now  dalina 


P  -  ^Tv"  "Relative  Volatility  of  l  wrt  2  (8) 


a  -  Y^X, 


(») 


substituting  Equation  9  in  Equation  8  and  noting  tho  rolationship 
in  Equation  4,  ws  obtain 


a 


(10) 


Equation  7b  oan  bs  writtsn  (substituting  Equation  9  and  Equation 
10) 


«k _ 

L  Xx  (a-1) 


_ 

JL _ _ 


i 


(11a) 


i.s 


•  * 


dL  _  dX±  tpjfr  +  (1-^)1 
L  Xx  (P-1)  (1-Xj) 


tfi r 


(lib) 
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integrating  the  above  equation,  we  get 


InL  -  [ln^  -  pind-jrj]  +  const 


(12a) 


i.e 


•  i 


ln< 


f  1-^ 
1  1-jfj 


(I2jb) 


or 


Jl.  . 

(M 

'  l-Jtf^ 

A' 

L  i 

■71 

1 1-Jfl  J 

iAr 


(12c) 


where  X1,,  and  L*  represent  some  known  values  (say,  the  initial 
values  in  the  pool) . 


4.0  MOLAR  EVAPORATION  RATE  OF  A  COMPONENT  OF  INTEREST 
4.1  Heat  Balance  Equation 

We  repeat  Equation  S  here  for  the  sake  of  continuity. 


(5b) 


-  Lht 
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substituting  for  dL/L  from  Equation  11,  wa  gat 


- 

dt 


Li <-£>  * 


jsa 

dt 


(1  -  + 


0f1+ A  (1-^)} 
(FT) 


{x/iT^r} 


(13) 


Wa  now  dafina  a  function  F  as  follows: 


(l-*i) 


-Lt P— p 

Xx  (l-xj 


(14) 


Substituting  Equations  12o  and  14  in  Equation  13  and  raarranging, 
wa  gat: 


dX, 


_ -..ALfcJ _ _ 

TJJFixJTUL1) 


with  xx  -  Xx  at  e  -  o 


(15) 


whara, 


4(t)  -  dQ/dt 


Nota  that  L/L*  is  a  function  of  X1  only. 

Tha  abova  aquation  oan  ba  solvad  by  Runga  Xutta  mathod  for 
apaoifiad  q(t). 


4.2  Maas  Rata  of  Evaporation 

Tha  molar  avaporation  of  epaoia  1  from  tha  pool  is  givan  by: 

^i,(t)  *  “  ^ dt'"'"  (1®*) 


whara,  n,(t)  is  tha  molar  avaporation  rata  of  componant  l. 
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i- . e  •  i 


M^t) 


dX.  r  dx, 

— -A  +  L  — pi 

dt  dt 


(16Jb) 


Substituting  for  Hj/L  and  dL/dX1  from  Equation  11  wa  get, 


-  L 


dxx 

dX, 

dt 

7FU  i-3  1 

i  *  i 

l^c7|T 

JL 


(17a) 


1  •  •  •  / 


^(t) 


rif  M  P  dX i 
L  1  i'J  7Mui-*i)  dt 


(i7Jb) 


4.3  Non-Dimensional izat ion  of  Equations 
Let 

tch  -  characteristic  time  in  the  problem  (18a) 

t  “  t/teh  ■  non-dimensional  time  (18b) 

Wa  can  determine  the  value  for  the  characteristic  time  by  imposing 
the  following  relation! 

tf<eoA>  tch  -  <1») 


The  above  equation  implies  that  teh  is  the  total  time  for  the 
evaporation  of  the  entire  pool  if  it  contained  only  component  1  and 
received  heat  at  a  rate  equal  to  that  at  time  teh. 


u-6 


<7(t)  -  <Ht )/dt(teh) 


( 20 ) 


Let, 


Hence,  Equation  15  reduces  to: 


•— T “  “  “  •»  fit  t  "  0  (21) 

ax  F(Xx)  (L/L1)  1 


Equation  17b  ie  now  written  aa: 


d,(C>  c«„ 

Li 


f-Ll 

UjJ 


J _  ja 

TFiTTi^T  dx 


(22) 


We  oan  write  an  equation  similar  to  equation  22  for  the  component  2 
also.  The  mass  rate  of  evaporation  of  component  2  la  given  by, 


^ch  ,t  /t  h  1  dXj 

“3 - (L/L>  TRTHraTtT  “37 


(23) 


where  0  is  defined  in  equation  6.  We  also  note  that, 

X2  -  1  -  X, 

Therefore,  equation  23  oan  be  written  as, 

dx. 


(24) 


Jklsh 

r.  i 


1 

irqjTX  “3T 


(25) 
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5.0  GROUND  NEXT  TRANSFER  MODEL 

The  ground  heat  transfer  rate  is  assumed  to  be  given  by  the  one 
dimensional  model  described  in  an  earlier  report  (see  Raj  and 
Morris,  1987).  This  model  assumes  the  following  formulation  for 
the  heat  transfer  rate  into  the  pool  of  liquid. 


<*"(t)-^  (25) 


We  also  choose  the  following  as  the  characteristic  time  for  the 
problem: 


'Cft 


h  r  1  J 

f  JhJL .] 

28  A  T 


(26) 


where,  AT  is  the  temperature  difference  between  the  ground  and  the 
evaporating  pool. 


6.0  NUMERICAL  EXAMPLE 


6.1  Values  of  Parameters 


We  choose  the  following  values  for  the  various  parameters.  The: 
volumes  of  the  liquids  spilled  are  equal  to  the  volumes  carried  on 
a  road  transport.  The  depth  of  the  pool  is  assumed  to  be  1  cm. 
The  area  of  the  pool  is  then  calculated  with  the  known  volume  of 
the  liquids  spilled. 


Volume  of  Fluorine  spilled  - 
Volume  of  Nitrogen  spilled  » 
Assumed  depth  of  pool  » 

Liquid  pool  temperature  • 

Ground  Temperature  » 

Ground  heat  transfer  parameter 
Relative  Volatility  (  p)  - 


1.771  mA3 
1.287  mA3 
0.01  m 

80.75  K  (avg  of  boiling  temp 
Of  F2  6  N2) 

293.0  K 

1400. 0  J/(mA2  eAl/2) 

0.412 
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The  values  of  some  of  the  calculated  parameter*  for  the  initial 
conditions  are  as  follows: 

Area  of  the  pool  ■  305.8  mA2 

Initial  molar  contents 

of  F2  ■  L,  (i)  -  72.66  kmoles 

of  N2  -  Lj(i)  -  36.36  kmoles 

Mole  Fraction  of  F2  =>  0.656 

in  the  liquid  [Xl(i)] 

Characteristic  time  (t„h)  -  15.51  s 


6.2  Results 

To  solve  the  equations  indicated  in  this  appendix,  a  FORTRAN 
program  was  developed.  Using  this  program,  results  were  obtained. 

The  result-,  ire  presented  in  the  attached  figures.  Figure  B-l 
shows  the  u.iation  of  the  molar  fractions  of  Fluorine  and  Nitrogen 
in  the  liquid  phase  as  a  function  of  time.  It  is  seen  that  the 
molar  r  itent  of  F2  increases  and  towards  the  end  of  the 
evapora  .on  the  liquid  is  rich  in  Fluorine.  This  is  because 
nitrogen  is  much  more  volatile  and  therefore  evaporates 
preferentially. 

Figure  B-2  shows  the  variation  of  the  mole  fractions  of  F2  and  N2 
in  the  vapor  phase  with  time.  Again  due  to  the  volatility  of 
nitrogen,  the  concentration  of  nitrogen  in  vapor  is  high  initially 
but  decreases  continuously  with  time. 

Figure  B-3  is  a  plot  of  the  molar  contents  of  the  liquid  pool  with 
time.  As  expected,  the  total  molar  content  decreases  with  time  due 
to  evaporation.  The  rate  of  decrease  of  nitrogen  in  the  liquid 
pool  is  higher  than  that  of  fluorine.  By  about  0.87  nondimensional 
time  units,  almost  all  of  the  nitrogen  has  evaporated. 

It  is  noticed  that,  at  about  0.95  nondimensional  time  (dimensional 
time  -  14.75  seconds),  all  of  the  liquid  in  the  pool  has 
evaporated.  That  is,  the  evaporation  of  the  F2  -  N2  pool  ie 
relatively  rapid. 


I  1  I  T  I  '  1  ’  I  '  I  '  I  '  l  '  I  1  I 

§§5§S2335§ 

9S0t)d  p.inbji  uj  suojpojj  9|o^ 
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dsoqcj  JodoA  ui  suojpojj  9|Ofl 
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Non  Dimensional  Time  [Tow) 

Figure  B.2:  Vapor  Mole  fraction  with  Time 


Figure  B.3:  Liq 
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APPENDIX  C 

THERMODYNAMIC  PROPERTIES  OF  HYDROGEN  FLUORIDE 

■10  3  0  HFXPROP.DAT 
VERSION  2.11 


112  0  MOLECULAR  WEIGHT  (XMWT) 

REF:  HACS 

20.006 

2  12  0  CRITICAL  POINT  TEMPERATURE  (TCRI)  K 
REF:  HACS 

461.00 

3  12  0  CRITICAL  POINT  PRESSURE  (PCRI)  N/m2 
REF:  HACS 

0.6480GE+07 

4  12  0  NORMAL  BOILING  POINT  (XNBP)  K 
REF:  HACS 

292.70 

5  12  0  NORMAL  FREEZING  POINT  (XNFP)  K 
REF:  HACS 

190.00 

6  8  2  0  VAPOR  HEAT  CAPACITY  (CPV)  PARAMETERS  J/kg  K 

REF*  HACS 

I.45300E  +  03  3.30500E-02  -1.0160E-04  1.25200E-07  600.00 

250.00  O.OOOOOE  +  00  O.OOOOOE+OO 

7  8  2  0  LIQ  HEAT  CAPACITY  (CPL)  PARAMETERS  J/kg  K 

REF:  HACS;  Upper  temp  set  to  TCRI  1 1/20/89  (was  293.00) 

3500.0  O.OOOOOE + 00  O.OOOOOE  +  OO  O.OOOOOE+OO  461.00 

265.00  3500.0  270.00 

8  7  2  0  LIQUID  DENSITY  (RHOL)  PARAMETERS  kg/m3 

REF:  HACS 

1.85300E+03  -3.9710  3.12500E-03  300.00  210.00 

O.OOOOOE+OO  O.OOOOOE+OO 

10  8  2  0  VAPOR  PRESSURE  (PSAT)  PARAMETERS  N/m2 

REF:  HACS 

-9.74369  4.68496  -2.98358  9,65825  0.0000 

461.00  273.16  1.00000 

11  2  2  0  ENTHALPY  OF  SAT  LIQUID  (HLIQS)  PARAMETERS  J/kg 

REF: 

O.OOOOOE+OO  292.7 

12  2  2  0  ENTHALPY  OF  SAT  VAPOR  (HVAPS)  PARAMETERS  J/kg 

REF: 

3.74400E+05  292.7 
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13  6  2  0  ENTHALPY  OF  VAPORIZATION  (XLAMDA)  J/kg 

REF:HACS  :XLAMDA  value  from  Allied  Chem  Co  publication  on  HF  properties. 

461.00  292.70  3.74400E+05  0.38000  461.00 

181.00 

14  6  2  0  SOLUBILITY  (SOL)  PARAMETERS  kg/ 100  kg 

REF:  MACS 

0.00000E  +  00  0.00000E+00  O.OOOOOE-t-OO  O.OOOOOE  +  OO  O.OOOOOE+O 

O.OOOOOE+OO 

15  1  2  0  ENTHALPY  OF  FUSION  (DHF)  i/kg 

REF:  HACS 

2.29000E+05 

16  1  2  0  ENTHALPY  OF  COMBUSTION  (DHC)  J/kg 

REF:  HACS 

O.OOOOOE  +  OO 

17  1  2  0  ENTHALPY  OF  DECOMPOSITION  (DHDC)  J/kg 

REF:  IIACS 

O.OOOOOE  +  OO 

18  1  2  0  ENTHALPY  OF  SOLUTION  (DHS)  J/kg 

REF:  HACS 

•0.30756E  +  07 

19  1  2  0  ENTHALPY  OF  REACTION  WITH  WATER  (DHWR)  J/kg 

REF:  HACS 

-0.30756E  +  07 

20  1  2  0  ENTHALPY  OF  POLYMERIZATION  (DHPY)  J/kg 

REF:  HACS 

O.OOOOOE  +  OO 

2112  0  AEROSOL  ENTRAINMENT  FRACTION  (XPHI) 

REF: 

1.0 

101  7  2  0  LIQ  THERMAL  COND  (XKL)  PARAMETERS  W/m  K 

REF:  HACS 

0.710  -8.6220E-04  -6.4400E-07  438.0  190.0 

0.00000E+00  O.OOOOOE+OO 

102  8  2  0  VAP  THERMAL  COND  (XKV)  PARAMETERS  W/m  K 

REF: 

3.85700E-03  5.27600E.05  2.26100E-08  O.OOOOOE+OO  1000.0 

175.000  O.OOOOOE  +  OO  0.00000E  +  00 

103  8  2  0  LIQ  VISCOSITY  (XMUL)  PARAMETERS  N  s/m2 

REF:  HACS 

•1.4040E  +  01  1.87900E  +  03  2.97500E-02  -3.0600E-05  453.00 

193.00  0.00000E  +  00  0.00000E  +  00 
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104  7  2  0  VAP  VISCOSITY  (XMUV)  PARAMETERS  N  s/m2 

REF:  Allied  Chemical  Corp;  "Hydrofluoric  Acid",  1978 

-1.0000E.06  4.50000E-08  0.00000E  +  00  373.000G+00  233.Q00E  +  0 

1.2000E-05  288.700E+00 

105  8  2  0  SURFACE  TENSION  (STEN)  N/m 

REF:  Allied  Chemical  Corp;  "Hydrofluoric  Acid",  1978 

8.6000E-03  461.0  273.15  1.555  297.2 

270.0  O.OOOOOE+OO  0.00000E+00 

106  2  2  0  INTERFACIAL  SURFACE  TENSION  (XITEN)  N/m 

REF:  HACS 

O.OOOOOE+OO  O.OOOOOE+OO 

201  l  2  0  LOWER  FLAM  ABILITY  LIMIT  (XLOFLM)  % 

REF:  HACS 

O.OOOOOE  +  OO 

202  1  2  0  UPPER  FLAMABILITY  LIMIT  (UPFLM)  % 

REF:  HACS 

O.OOOOOE  +  OO 

203  1  2  0  BURN  RATE  (BRAT)  m/s 

REF:  HACS 

O.OOOOOE+OO 

204  1  2  0  ADIABATIC  FLAME  TEMPERATURE  (ADFT)  K 

REF:  HACS 

0.00000E  +  00 

205  1  2  0  MOLAR  RATIO  REACTANTS/PRODUCTS  (MRAT) 

REF:  HACS 

O.OOOOOE+OO 

206  1  2  0  AIR  FUEL  RATIO  (AFRT) 

REF:  HACS 

O.OOOOOE  +  OO 

207  1  2  0  FLAME  TEMPERATURE  (FLTM)  K 

REF:  HACS 

0.00000E+00 

211  1  2  0  LIMITING  VALUE  OF  MOLECULAR  FN  CONC  (XMFRC) 

REF:  HACS 

O.OOOOOE+OO 

212  0  2  0  GREY  BODY  FIRE  EMISSIVE  POWER  (kW/nT2) 

REF: 

213  0  2  0  EMISSIVITY  (fraction) 

REF: 

214  0  2  0  EFFECTIVE  FIRE  TEMPERATURE  (K) 

REFs 

301  1  2  0  TOXIC  INHALATION  LIM1T/TLV  (TOX1)  ppm 

REF:  HACS 

3.0000 
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302  1  2  0  SHORT  TERM  INHALATION  LIMIT  (TOX2)  ppm 
REF*  HACS 

3.0000 

303  1  2  0  SHORT  TERM  INHALATION  TIME  LIMIT  (TOX3)  s 
REF:  HACS 

900.00 

304  1  2  0  LOWER  TOXICITY  INGESTION  LIMIT  (TOX4)  kg/kg 
REF:  HACS 

O.MOOOOE-04 

305  1  2  0  UPPER  TOXICITY  INGESTION  LIMIT  (TOX5)  kg/kg 
REF:  HACS 

O.OOOOOE  +  OO 

306  0  2  1  LATE  TOXICITY  LEVEL  (TOX6) 

REF:  HACS 

800  1  2  1  NAME  OF  CHEM  (CNAM) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL. DAT 
HYDROGEN  FLUORIDE 

801  l  2  1  SHIPPING  STATE  CODE 

REF:  HACS 

L 

802  1  2  1  HAZARD  CLASS  (HACL) 

REF:  COMPUTERIZED  DATA  IN  CHRIS_ALL.DAT 
Corrosive 

803  1  2  1  HAZARD  CLASS  NUMBER  (HACN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 

8 

804  1  2  1  DOT.  NUMBER  (DOTN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 

1052 

805  1  2  1  PATH  CODE  (PATC) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 
A-C-K-M-N-0 

806  0  2  1  ENVIRONMENT  BEHAVIOR  CODE 


CM 
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APPENDIX  C  (contd) 

THERMODYNAMIC  PROPERTIES  QE, FLUORINE 

-10  3  0  FXXPROP.DAT 
VERSION  2.11 


112  0  MOLECULAR  WEIGHT  (XMWT) 

REF:  HACS 
37.997 

2  12  0  CRITICAL  POINT  TEMPERATURE  (TCRI)K 
REF:  HACS 

144.30 

3  12  0  CRITICAL  POINT  PRESSURE  (PCRI)  N/m2 
REF:  HACS 

0.55200E  +  07 

4  12  0  NORMAL  BOILING  POINT  (XNBP)  K 
REF:  HACS 

84.45 

5  12  0  NORMAL  FREEZING  POINT  (XNFP)  K 
REF:  HACS 

53.500 

6  8  2  0  VAPOR  HEAT  CAPACITY  (CPV)  PARAMETERS  J/kg  K 

REF:  HACS 

611.00  0.96240  -9.50900E-04  3.16900E-07  1000.0 

84.45  O.OOOOOE+OO  O.OOOOOE+OO 

7  8  2  0  LIQ  HEAT  CAPACITY  (CPL)  PARAMETERS  J/kg  K 

REF:  HACS 

1087.0  1.79  O.OOOOOE+OO  O.OOOOOE  +  OO  110.0 

55.0  1238.26  84.5 

8  7  2  0  LIQUID  DENSITY  (RHOL)  PARAMETERS  kg/m3 

REF:  HACS 

1559.3  -0.73000  O.OOOOOE  +  OO  173.16  84.450 

1500.0  84.450 

10  8  2  0  VAPOR  PRESSURE  (PSAT)  PARAMETERS  N/m2 

REF:  HACS 

-6.18224  1.18062  -1.16555  -1.50167  O.OOOOOE+O 

144.30  64.000  1.0 

11  2  2  0  ENTHALPY  OF  SAT  LIQUID  (HLIQS)  PARAMETERS  J/kg 

REF: 

O.OOOOOE+OO  84.45 

12  2  2  0  ENTHALPY  OF  SAT  VAPOR  (HVAPS)  PARAMETERS  J/kg 

REF: 

0.17300E+06  84.45 


C-3 


FXXPROP.DAT:  page  2  of 

13  6  2  0  ENTHALPY  OF  VAPORIZATION  (XLAMDA)  J/kg 

REF:  HACS 

144.30  84.450  0.17300Ci  06  0.38000  144.30 

54.000 

14  6  2  0  SOLUBILITY  (SOL)  PARAMETERS  kg/ 100  kg 

REF:  HACS 

0.00000E+00  O.OOOOOE+OO  G.000Q0E+00  O.OOOOOE+OO  O.OOOOOE+O 

O.OOOOOE+OO 

15  1  2  0  ENTHALPY  OF  FUSION  (DHF)  J/kg 

REF:  HACS 

4.1000CE  +  04 

16  l  2  0  ENTHALPY  OF  COMBUSTION  (DHC)  J/kg 

REF:  HACS 

0.O0000E+00 

17  1  2  0  ENTHALPY  OF  DECOMPOSITION  (DHDC)  J/kg 

REF;  HACS 

O.OOOOOE  +  OO 

18  1  2  0  ENTHALPY  OF  SOLUTION  (DHS)  J/kg 

REF:  HACS 

0.O0000E  +  00 

19  1  2  0  ENTHALPY  OF  REACTION  WITH  WATER  (DMWR)  J/kg 

REF:  HACS 

O.OOOOOE  +  OO 

20  1  2  I)  ENTHALPY  OF  POLYMERIZATION  (DHPY)  J/kg 

REF:  HACS 

O.OOOOOE  +  OO 

21  1  2  0  AEROSOL  ENTRAINMENT  FRACTION  (PHI) 

REF: 

1.0 

101  7  2  0  LIQ  THERMAL  COND  (XKL)  PARAMETERS  W/m  K 

REF:  HACS 

0.2565  •6.7950E-04  •4.9580E-06  133.0  54.0 

().t)0(K)0E  +  0()  O.OOOOOE  +  OO 

102  8  2  0  VAP  THERMAL  COND  (XKV)  PARAMETERS  W/m  K 

REF: 

7.8 1200E-04  8.28700E-05  5.19300E-08  O.OOOOOE  +  OO  600.0 

145.0  0.0000QE+00  0.00000E+00 

103  8  2  0  LIO  VISCOSITY  (XMUL)  PARAMETERS  N  »/m2 

REF:  HACS 

•3.6290  1 .97200 E+ 02  -9.3780E-04  -6.2750E06  88.0 

54.0  O.OOOOOE  +  OO  O.OOOOOE  +  OO 
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104  7  2  0  VAP  VISCOSITY  (XMUV)  PARAMETERS  N  s/m2 
REF: 

2.09000E-04  O.OOOOOE+OO  O.OOOOOE+OO  1000.0  0.0 

,  O.OOOOOE+OO  O.OOOOOE  +  OO 

105  8  2  0  SURFACE  TENSION  (STEN)  N/m 

REF:  MACS 

0.01460E+00  144.300E+00  81.0000E+00  1.00000E+00  100.000E+0 

55.0000E+00  O.OOOOOE+OO  O.OOOOOE+OO 

106  2  2  0  INTERFACIAL  SURFACE  TENSION  (XITEN)  N/m 

REF:  HACS 

O.OOOOOE+OO  O.OOOOOE+OO 

201  1  2  0  LOWER  FLAMABILITY  LIMIT  (XLOFLM)  % 

REF:  HACS 

0.00000E+00 

202  1  2  0  UPPER  FLAMABILITY  LIMIT  (UPFLM)  % 

REF:  HACS 

O.OOOOOE+OO 

203  1  2  0  BURN  RATE  (BRAT)  m/s 

REF:  HACS 

O.OOOOOE+OO 

204  1  2  0  ADIABATIC  FLAME  TEMPERATURE  (ADFT)  K 

REF:  HACS 

0.00000E+00 

205  1  2  0  MOLAR  RATIO  REACTANTS/PRODUCTS  (MRAT) 

REF:  HACS 

0.00000E+00 

206  1  2  0  AIR  FUEL  RATIO  (AFRT) 

REF:  HACS 

O.OOOOOE  +  OO 

207  1  2  0  FLAME  TEMPERATURE  (FLTM)  K 

REF:  HACS 

O.OOOOOE  +  OO 

211  1  2  0  LIMITING  VALUE  OF  MOLECULAR  FN  CONC  (XMFRC) 

REF:  HACS 

O.OOOOOE+OO 

301  1  2  0  TOXIC  INHALATION  LIMIT/TLV  (TOXt)  ppm 

REF:  HACS 

1.0000 

.  302  1  2  0  SHORT  TERM  INHALATION  LIMIT  (TOX2)  ppm 

REF:  HACS 

0.50000 

-  303  1  2  0  SHORT  TERM  INHALATION  TIME  LIMIT  (TOX3)  s 

REF:  HACS 

300.00 
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304  1  2  0  LOWER  TOXICITY  INGESTION  LIMIT  (TOX4)  kg/kg 
REF:  HACS 

O.OOOOOE+OO 

305  1  2  0  UPPER  TOXICITY  INGESTION  LIMIT  (TOX5)  kg/kg 
REF:  HACS 

O.OOOOOE  +  OO 

306  0  2  1  LATE  TOXICITY  LEVEL  (TOX6) 

REF:  HACS 

800  1  2  1  NAME  OF  CHEM  (CNAM) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 
FLUORINE 

801  0  2  1  SHIPPING  STATE  CODE 

REF:  HACS 

802  1  2  1  HAZARD  CLASS  <HACL) 

REF:  COMPUTERIZED  DATA  IN  CHRIS_ALL.DAT 
Nonflammable  gus 

803  1  2  1  HAZARD  CLASS  NUMBER  (HACN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 

2 

804  1  2  1  DOT.  NUMBER  (DOTN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 

1045 

805  1  2  1  PATH  CODE  (PATC) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 
A-C 

806  0  2  1  ENVIRONMENT  BEHAVIOR  CODE 

REF: 
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APPENDIX  C  (contd) 

THERMODYNAMIC  PROPERTIES  OF  NITROGEN 

,  -1  0  3  0  N2.DAT  NITROGEN  ACIJ 

VERSION  2.11 


112  0  MOLECULAR  WEIGHT  (XMWT) 

REF:  MACS 
28.013 

2  12  0  CRITICAL  POINT  TEMPERATURE  (TCRI)  K 
REF:  MACS 

126.20 

3  12  0  CRITICAL  POINT  PRESSURE  (PCRI)  N/m2 
REF:  MACS 

33.9E+05 

4  12  0  NORMAL  BOILING  POINT  (XNBP)  K 
REF:  MACS 

77.40 

5  12  0  NORMAL  FREEZING  POINT  (XNFP)  K 
REF:  MACS 

63.30 

6  8  2  0  VAPOR  MEAT  CAPACITY  (CPV)  PARAMETERS  J/kuK 

REF:  REID 

1112.00  -0.8440  9.5670E-04  -4.1690E-07  600.00 

298.00  1.0420E+03  298,00 

7  8  2  0  LIO  MEAT  CAPACITY  (CPL)  PARAMETERS  J/kg  K 

REF:  YAWS, C.L., "PHYSICAL  PROPERTIES", MCGRAW-HILL, NY, P.220. 

2309.00  -11.100  0.100  O.OOOOE+OO  90.01) 

70.00  2.0490E+03  77.4 

8  7  2  0  LIQUID  DENSITY  (RHOL)  PARAMETERS  kg/m3 

REF:  MACS  (THI  from  CEH  p.3-167  (wus  293)) 

790.9  4.7700  -0.05818  120.00  70.00 

812.0  126.2 

10  8  2  0  VAPOR  PRESSURE  (PSAT)  PARAMETERS  N/m2 

REF:  REID  , P.658, N0.22 

*6.09676  1.13670  -1.04072  -1.93306  O.OOOOE  +  OO 

126.20  63.0  1.0 

11  2  2  0  ENTHALPY  OF  SATURATED  LIQUID  (I ILIQS)  J/kg 

REF:  BASIS 

.00000  273.16 

12  2  2  0  ENTHALPY  OF  SATURATED  VAPOR  (HVAPS)  J/kg 

REF:  MACS 

O.OOOOE+OO  O.OOOOE+OO 
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13  0  2  0  ENTHALPY  OF  VAPORIZATION  (XLAMDA)  J/kg 
REF:  HACS 

126.20  77.40  1.9900E  +  05  .38000  1 10.0 

75.00 

14  6  3  0  SOLUBILITY  (SOL)  PARAMETERS  kg/100  kg 

REF:  SCONCE, J.S., "CHLORINE", KREIGER  PUBLISI  IINC, HUNTINGTON, NY, t972.P.33. 
REF:  FLDS  1-4  HACS  FLDS  5,6 
5.3927  -.01569  323.16  282.77  .65000 

298.16 

15  1  2  0  ENTHALPY  OF  FUSION  (DKF)  J/kg 

REF: 

.90330E+05 

16  12  0  ENTHALPY  OF  COMBUSTION  (DHC)  J/kg 
REF:  HACS 

.00000 

17  1  2  0  ENTHALPY  OF  DECOMPOSITION  (DHDC)  J/kg 

REF:  HACS 

.00000 

18  1  2  0  ENTHALPY  OF  SOLUTION  (DHS)  J/kg 

REF:  HACS 

.00000 

19  1  2  0  ENTHALPY  OF  REACTION  WITH  WATER  (DHWR)  J/kg 

REF:  HACS 

.00000 

20  1  2  0  ENTHALPY  OF  POLYMERIZATION  (DHPY)  J/kg 

REF:  HACS 

.00000 

2112  0  AEROSOL  ENTRAINMENT  FRACTION  (PHI) 

REF: 

1.0 

101  7  2  0  LIQ  THERMAL  COND  (XKL)  PARAMETERS  W/m  K 

REF:  YAWS,  P.224. 

2.629E-01  -1.5450E-03  -9.450E-07  121,0  64.00 

1.1500E-01  91.0 

102  8  2  0  VAP  THERMAL  COND  (XKV)  PARAMETERS  W/m  K 

REF:  YAWS  P.208. 

3.919E-04  9.8160E-05  -5.0670E-08  0.U000E  +  00  1470.00 

115.0  2.5100E-02  298.0 

103  8  2  0  LIQ  VISCOSITY  (XMUL)  PARAMETERS  N  s/m2 

REF'  YAWS  P2P 

-2.795E  +  01  ’  8.660E+02  2.7630E-01  -1.0840E-03  78.00 

68.00  1. 80000-04  73.0 
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104  7  2  0  VAP  VISCOSITY  (XMUV)  PARAMETERS  N  s/m2 

REF:  YAWS,  P.210. 

4.2150E-06  O.OOOOE+OO  4.5000E-08  373.0  173.00 

«  1.2000E-05  173.00 

105  8  2  0  SURFACE  TENSION  (STEN)  N/m 
REF:  YAWS,  P.218.(PARS  1-6)  IIACS(PARS  7,8) 

8.750E-03  126.2  77.80  1.235  100.00 

78.00  O.OOOOE+OO  O.OOOOE+OO 

106  2  2  0  INTERFACIAL  SURFACE  TENSION  (XITEN)  N/m 

REF:  HACS 

.00000  .00000 

201  1  2  0  LOWER  FLAMABILITY  LIMIT  (XLOFLM)  % 

REF:  HACS 

.00000 

202  1  2  0  UPPER  FLAMABILITY  LIMIT  (UPFLM)  % 

REF:  HACS 

.00000 

203  1  2  0  BURN  RATE  (BRAT)  m/s 

REF:  HACS 

.00000 

204  1  2  0  ADIABATIC  FLAME  TEMPERATURE  (ADFT)  K 

REF:  HACS 

.00000 

205  1  2  0  MOLAR  RATIO  REACTANTS/ PRODUCTS  (MRAT) 

REF:  HACS 

.00000 

206  1  2  0  AIR  FUEL  RATIO  (AFRT) 

REF:  HACS 

.00000 

207  1  2  0  FLAME  TEMPERATURE  (FLTM)  K 

REF:  HACS 

.00000 

301  1  2  0  TOXIC  INHALATION  LIMIT/TLV  (TOX1)  ppm 

REF;  HACS 

1.0000 

302  I  2  0  SHORT  TERM  INHALATION  LIMIT  (TOX2)  ppm 

REF:  HACS 

3.0000 

*  303  1  2  0  SHORT  TERM  INHALATION  TIME  LIMIT  (TOX3)  s 

REF:  HACS 

300.00 

.  304  1  2  0  LOWER  TOXICITY  INGESTION  LIMIT  (TOX4)  kg/kg 

REF:  HACS 

.00000 
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305  1  2  0  UPPER  TOXICITY  INGESTION  LIMIT  (TOX5)  kg/kg 

REF:  HACS 

.00000 

306  0  2  1  LATE  TOXICITY  LEVEL  (TOX6) 

REF:  MACS 

701  1  2  0  LIMITING  VALUE  OF  MOLECULAR  FN  CONC  (XMFRC) 

REF:  HACS 

.00000 

800  1  2  1  NAME  OF  CHEM  (CNAM) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 
NITROGEN 

801  1  2  1  SHIPPING  STATE  CODE 

REF:  HACS 

L 

802  1  2  1  HAZARD  CLASS  (HACL) 

REF:  COMPUTERIZED  DATA  IN  CHRIS.  ALL.DAT 
Nonflammable  gas;  Poison 

803  1  2  1  HAZARD  CLASS  NUMBER  (IlACN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS.  ALL.DAT 

2  and  6 

804  1  2  1  DOT.  NUMBER  (DOTN) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 

1017 

H05  1  2  1  PATH  CODE  (PATC) 

REF:  COMPUTERIZED  DATA  IN  CHRIS  ALL.DAT 
A-C-I-J 

806  1  2  1  ENVIRONMENT  BEHAVIOR  (ENVB) 

REF: 

V\V 
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